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Chapitre I 


Introduction aux concepts fondamentaux 
de la mécanique quantique 


I- Rappel succinct des théories de la physique classique 

Avant la naissance de la physique quantique vers le début du 20 ème siècle, la 
physique classique distinguait deux catégories d'objets, les corpuscules d'une part et les 
ondes d'autre part : 

- les corpuscules sont des entités matérielles discrètes et localisables dans l'espace à 
tout instant; ils sont régis par les lois classiques de la mécanique newtonienne (1687). 

- les ondes sont des phénomènes continus, non localisables et occupent tout l'espace à 
tout instant ; elles sont décrites par la théorie des ondes électromagnétiques de 
Maxwell (1867). 

1-1° Mécanique newtonienne 

En mécanique classique, le mouvement d'un point matériel de masse m se déplaçant 
avec une vitesse v sous l'action d'une force F est régi par la seconde loi de Newton F = my . 
En terme de la quantité de mouvement (ou impulsion) p = mv ou bien en terme de l'énergie 
potentielle V(x,y,z), la loi fondamentale de la mécanique classique peut s'écrire sous les 
deux formes suivantes : 


Dans le cas où F = -VV , on dit que la force dérive d'un potentiel et la relation n'étant 
valable que dans le cas d'un système conservatif, c'est-à-dire un système pour lequel son 
énergie totale E = T + V (somme des énergies cinétique T et potentielle V) reste constante 
dans le temps, (somme des énergies cinétique T et potentielle V) reste constante dans le 
temps. 

Lors des applications de la mécanique classique, on suppose toujours que la masse m 
du point matériel demeure constante quelque soit la vitesse v avec laquelle elle se déplace. 
Lors des applications de cette mécanique, les succès ne sont obtenus que dans le cas des 
corps macroscopiques ayant des grandes masses et des vitesses relativement petites. Or, 
pour des corps ayant des petites masses se déplaçant avec des vitesses relativement 
grandes, cette mécanique cesse d'être applicable et en alternative on fait appel soit à la 
mécanique relativiste d'Einstein soit à la mécanique quantique. 

1-2° Mécanique relativiste d'Einstein 

En 1905, Einstein pose les bases de la théorie de la relativité restreinte qui a été 
développée indépendamment de la mécanique quantique. Il montre que la masse d'un 
corps en mouvement varie avec sa vitesse de déplacement. Si mo est la masse de la particule 
au repos, alors sa masse m lorsqu'il acquiert une vitesse v relativement grande est : 


( dV/dx' 



F = -VV = dV/dy 
K dV/dZ; 


dt 
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La masse m du corps augmente avec la vitesse v et tend vers l'infini lorsque v tend vers la 
vitesse c de la lumière (3x1 0 8 m/s). On doit donc admettre que c est la limite supérieure 
assignée à la vitesse v de tout corps (v<c), c'est-à-dire qu'aucun objet matériel en 
mouvement ne peut dépasser la célérité c de la lumière. Dans le cas d'un photon qui se 
déplace dans le vide à la vitesse v = c , sa masse au repos m 0 = m (1 - v 2 / c 2 ) 1/2 est nulle. 

Selon la théorie de la relativité, il existe une équivalence entre la masse et l'énergie 
qui se traduit par la célèbre relation : 

E = me 2 


qui veut dire que toute masse m se déplaçant à la vitesse de la lumière se transforme en 
énergie et vice versa. Remplaçons m par sa relation : 


Vl-v 2 /c 2 

et portons l'expression de l'énergie au carré, on obtient : 

m 2 ^4 


Pour une particule matérielle se déplaçant à une vitesse v assez limitée, (v/c) 2 (( 1, on peut 
effectuer le développement approché suivant : 

E 2 = m?c 4 fl + 4 


Sachant que p=mov, l'équation précédente conduit à la célèbre relation de l'énergie 
relativiste : 


E 2 = hIqC 4 +p 2 c 2 


ou encore : 


/ 9 2 9 4 

y p c + m 0 c = m 0 c 


v m o c 


E = m 0 C 1 + -^-2 
2ïï1qC 2 


Pour trouver l'origine du terme d'énergie relativiste, effectuons un développement limité 
de E en puissance de (p/moc). On obtient : 

P 4 

8moC 4 

ou encore : 

E = m 0 c 2 + L- 

2m 0 8moC 

On voit que l'énergie d'une particule relativiste se compose, en plus de l'énergie au repos 
(m Q c 2 ) et de l'énergie cinétique non-relativiste (p 2 /2m Q ), d'un terme supplémentaire 
(- p 4 /8m ( jC 2 ) due à la variation relativiste de la masse avec la vitesse. 

Selon la théorie de la relativité, toute force F à laquelle est soumis un corps doit 
également varier avec la masse. En effet : 

- dp d(mv) dv dm _ dm 

F = — = = m — + v = my + v 

dt dt dt dt dt 
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Aux faibles vitesses (v«c), lorsque les variations de la masse sont négligeables 
(dm/dt-O), on retrouve la formulation habituelle F = my. Ainsi, la relation F = my ne 
s'applique qu'aux objets non relativistes. 


1-3° Théorie des ondes électromagnétiques de Maxwell 


En toute généralité, les ondes électromagnétiques peuvent être émises par diverses 
sources: soleil, corps chauffés, corps incandescents, tubes à décharge, ... Ces ondes 
comprennent en toute généralité les rayons y , les rayons X, l'ultraviolet, l'infrarouge et les 
ondes radios. La lumière n'est autre qu'une onde électromagnétique qui lorsqu'elle se 
propage dans le vide, elle acquiert une vitesse limite c = 3xl0 8 m/s. 

Selon la théorie de Maxwell, une onde se propageant dans l'espace (x,y,z) et dans le 
temps (t) obéit à l'équation d'onde suivante : 


1 d 2x ¥ 


- A = /ôx^ + ô ^ / ôy^ + ô ^ 18 : le Laplacien, 

- : l'amplitude ou fonction d'onde, 

- v : la vitesse de propagation de l'onde (v=c dans le vide). 


Les caractéristiques de l'onde sont bien définies, si l'on connaît : 

- sa longueur d'onde À, = v/ v ou sa pulsation co = 2nv , 
son intensité 1= | T 7 1 2 = V P* V P, avec 'F* le conjugué de T. 

Si l'onde est supposée monochromatique, c'est-à-dire se propageant dans une seule 
direction (par exemple x), l'équation d'onde devient simplement : 

d 2v F(x,t) _ 1 d 2 ¥(x,t) 

Sx 2 v 2 St 2 


Dans le cas d'une onde stationnaire, la fonction T^xd) peut se mettre sous la forme de 
produit de deux fonctions indépendantes : 

T(x,t) = O(x) T(t) 

Le remplacement de ce produit dans l'équation et la division par 0T conduit à : 

V 2 d 2 P _ 1 ô 2 T 

®ÔÔdx 2 


Le membre de gauche de cette relation ne contient que la variable d'espace (x) et celui de 
droite ne contient que la variable temps (t). La seule et unique façon d'égaler ces deux 
expressions différentielles est de les considérer chacune comme constante. Si cette constante 
commune est désignée par -co 2 , nous obtenons deux équations différentielles 
indépendantes : 


d 2 $ 

CO 2 


d 2 © 

(7.1 if 


+ — r-0 = 0 

ou 


+ 0=0 

dx 2 

v 2 


dx 2 

U J 

d 2 T 

d# 

+ 

s 

H 

o 

ou 

d 2 T 

dt 2 

+ (2ttv) 2 T = 0 
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Leurs solutions générales sont de la forme : 

j 2nx ,j 2sx 

o(x) = Aje ^ + A 2 e ^ 

T(t) = B,e i2,tvt + B 2 c~ i2nvt 


où Ai, A 2 , Bi et B 2 sont des constantes arbitraires qui se déterminent à partir des conditions 
initiales et des conditions aux limites pour les variables x et t. Si l'onde est supposée se 
propager dans la direction x positive (coefficient A 2 = 0) et que son amplitude s'annule à 
t oo (coefficient Bi = 0), alors les fonctions spatiale et temporelle s'écrivent : 
j2rot 

O(x) = Aje ^ T(t) = B 2 e _l27rvt 


Dans ces conditions, l'amplitude de l'onde ^(xA) devient : 

T(x,t) = Ce ^ ou ^(xa) = O(x) e _i27t vt 


La fonction O(x) est solution de l'équation spatiale : 
d 2 0 47T 2 

dx 2 X 2 


-0 = 0 


qui décrit le phénomène de propagation de l'onde stationnaire dans l'unique direction x. Si 
l'on suppose maintenant que l'onde se propage dans toutes les directions (x,y,z), le même 
raisonnement précédent conduit à l'équation stationnaire suivante : 


AO + 


471 2 


0=0 


avec A étant le Laplacien et 0(x,y,z) l'amplitude de l'onde. Dans la théorie de Maxwell, 
l'énergie E d'une onde électromagnétique est supposée seulement proportionnelle à son 
intensité I : 

E = k I 


: 2nx 
Aje x 


^ 2kx 


Aje 


Aj Aj = |Aj| 


L'énergie E ne dépend donc que de l'intensité I de l'onde, elle est indépendante de sa 
fréquence v ou de sa longueur d'onde X , par suite elle est continue. 


1-4° Insuffisances des théories de la physique classique 

La physique classique ne trouve son succès que lorsque l'on applique aux objets 
macroscopiques. Cependant, dans le cas des objets microscopiques (atomes, molécules, ou 
encore des objets plus petits appelés particules élémentaires) doués de petites masses ou de 
grandes vitesses, l'application de cette physique se heurte à des difficultés insurmontables. 
La contestation la plus remarquable à relever est que cette physique prévoit que l'énergie 
de la matière, l'énergie du rayonnement et l'énergie d'échange rayonnement-matière sont 
continues : 

E classique = continue 

Or, en réalité ces énergies sont discontinues, elles ne prennent que certaines valeurs 
permises. 
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Le cas le plus illustratif de l'insuffisance de la physique classique est l'étude de la 
stabilité des édifices atomiques ou moléculaires. Si les lois de cette physique s'appliquaient 
à un électron gravitant autour d'un noyau, celui-ci devrait émettre continuellement des 
ondes électromagnétiques en perdant de l'énergie. Ainsi, l'électron rayonnant de l'énergie 
tomberait sur le noyau et l'atome serait instable. Parmi les célèbres expériences qui ont mis 
la physique classique en échec, on cite : 

- le rayonnement du corps noir, 

- l'effet photoélectrique, 

- le spectre de raies atomiques 

Ces trois expériences ont permis une révision radicale de la physique classique de l'époque 
pour l'établissement d'une nouvelle théorie de la physique dénommée "physique 
quantique" ou "mécanique quantique". La mécanique quantique est la théorie qui décrit le 
monde microscopique à l'échelle des atomes, des molécules et des particules élémentaires. 
Comme toute nouvelle théorie, la mécanique quantique repose sur un certain nombre 
d'hypothèses et principes. 


II- Hypothèses et principes fondamentaux de la mécanique quantique 


11-1° Hypothèse de Planck - Rayonnement du corps noir (1900) 

Par définition, un corps noir est un corps capable d'absorber tout rayonnement reçu 
quelque soit sa fréquence v , et inversement il restitue le même rayonnement par émission. 
Pour interpréter l'expérience du rayonnement de ce corps, impossible à résoudre par les 
théories classiques, Planck stipule l'hypothèse suivante : 

Les échanges d'énergie entre la matière et le rayonnement ne peuvent 
s'effectuer que par quantités discrètes appelées quanta d'énergie : 
E = hv 

où h est une constante de proportionnalité universelle appelée constante de Planck, elle 
vaut 6.6262 x 10-34 J.s. 


Expérience de rayonnement du corps noir 

Il est bien connu qu'un corps solide porté à une certaine température émet de 
l'énergie sous forme de rayonnement électromagnétique de fréquence v spécifique. Le 
modèle idéal d'un corps noir consiste à prendre une enceinte sphérique couverte à 
l'intérieur d'une couche assez épaisse de noir et menée d'un petit orifice laissant échapper 
la radiation à mesurer : 
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Les caractéristiques du spectre émis par un corps solide dépendent de la température. 
Lorsque la température T s'élève, la fréquence v de la radiation mesurée augmente et l'on 
passe d'une radiation rouge à une radiation jaune puis blanche. A une température donnée, 
l'intensité de la radiation émise croît au départ, elle passe par un maximum puis elle décroît 
en fonction de la fréquence. 

Les lois classiques connues de l' électromagnétisme et de la thermodynamique n'ont 
pas pu expliquer les propriétés du spectre du corps noir. Selon ces lois, les échanges 
d'énergie à l'intérieur d'un corps s'effectuent de manière continue et proportionnelle à la 
température T : 

E f - Ei = kT 

Ei : énergie de l'état initial avant absorption de l'énergie thermique, 

Ef : énergie de l'état final après absorption de cette énergie, 
k : constante de Boltzmann valant 1.3807 x 10‘ 23 J/ K. 

Ainsi, le corps dans un état final (Ef) revient à l'état initial (Ei) en rayonnant de l'énergie. 
Cette énergie rayonnée selon les lois classiques augmente indéfiniment avec la température 
et que le spectre émis est continu (toutes les fréquences devraient être représentées dans le 
spectre). Cette prévision classique se trouve alors en complet désaccord avec les résultats 
expérimentaux puisque le spectre d'émission d'un corps solide était formé de raies 
distinctes (spectre discontinu). Ces contradictions ont incité à une révision en profondeur 
les théories classiques de l'époque. 

Grâce aux travaux de Planck, ces contradictions ont été résolues et une nouvelle 
théorie est née en 1900. Il montre que l'état final d'un corps ne peut avoir n'importe quelle 
énergie, mais que cette énergie doit avoir des valeurs déterminées non pas en fonction de la 
température mais en fonction de la fréquence v émise selon la relation : 

Ef - Ei = hv 

En d'autre terme, un corps ne peut absorber ou émettre de l'énergie que par quantités 
discrètes hv appelées quanta. Ainsi, Planck introduit pour la première fois la notion de 
quantification de l'énergie qui a ouvert la voie vers l'établissement d'une nouvelle 
mécanique. 


11-2° Hypothèse d'Einstein - Effet photoélectrique (1905) 

L'effet photoélectrique a été découvert en premier par Hertz en 1887 et aucune 
théorie de la physique classique de l'époque n'a permis de l'expliquer. Le phénomène de 
l'effet photoélectrique consiste en l'aptitude d'une lumière incidente à arracher des 
électrons d'une surface métallique. Pour explique ce phénomène, Einstein émis l'hypothèse 
suivante : 


La lumière est constituée d'un flux de particules d'énergie hv appelées photons. 


Par cette hypothèse, Einstein montre clairement l'aspect corpusculaire de la lumière (ou de 
toute onde électromagnétique en générale) et que l'énergie de la lumière est quantifiée. En 
percutant les électrons de la surface métallique, sous l'effet des chocs, les photons ayant une 
énergie hv suffisante arrachent les électrons de cette surface. 
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Expérience de l'effet photoélectrique : 


Il est bien connu que les métaux sont conducteurs et présentent à la surface des 
électrons de nature très mobiles. Un électron ne peut être éjecté du métal que s'il reçoit 
d'une façon ou d'une autre une certaine énergie suffisante puisque cet électron est lié aux 
noyaux constituants le solide. 

Le principe de l'expérience de l'effet 
photoélectrique est illustré par le dispositif simple ci- 
contre. La surface métallique est considérée comme une 
cathode recevant une lumière et une anode est disposée 
pour recueillir les électrons. Pour mesurer l'aptitude du 
métal à injecter des électrons, un courant électrique et une 
différence de potentiel entre les deux électrodes sont 
mesurés. L'expérience a montré que l'effet photoélectrique 
ne se produit que lorsque la fréquence de la lumière 
incidente v est supérieure à une certaine valeur limite dite fréquence-seuil vs : 
si v > vs l'effet se produisant spontanément 
si v < vs l'effet ne se produisant pas 

L'extraction des électrons ne dépend donc que de la fréquence de la lumière et non de son 
intensité. La théorie classique de la lumière se trouve incapable d'expliquer l'effet 
photoélectrique puisque l'énergie d'une onde lumineuse est supposée seulement 
proportionnelle à son intensité indépendamment de sa fréquence. 

Devant ces difficultés de la physique classique, Einstein a pu expliquer l'effet 
photoélectrique à partir des collisions entre les photons d'énergie hv et les électrons du 
métal. Les électrons sont liés au métal par une certaine énergie Ws, appelée énergie de 
sortie ou énergie de liaison. Cette énergie est nécessaire pour rompre la liaison qui 
maintient l'électron au métal. Après extraction, l'électron se trouve libre avec une énergie 
cinétique V2 mv 2 . 




1 surface du méta 


A partir de la loi de conservation de l'énergie, on obtient l'équation dite équation d'Einstein 
pour l'effet photoélectrique : 

hv = W s +|mv 2 

On peut donc identifier la fréquence-seuil observée expérimentalement par v s = W s / h . La 
fréquence vs dépend bien de la nature du métal, chaque métal possède une valeur de Ws 
caractéristique. D'où : 

h(v-v s ) = |mv 2 

Comme l'énergie cinétique est un terme positif, il est clair que le phénomène de l'effet 
photoélectrique ne peut se produire que si le métal reçoit assez d'énergie lumineuse 
(hv > W s ), c'est-à-dire une lumière incidente de fréquence assez suffisante (v > vs). 
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En résumé l'hypothèse d'Einstein, montre que la lumière ne manifeste pas 
uniquement l'aspect ondulatoire de propagation de son onde, mais également l'aspect 
corpusculaire photonique qui emmagasine son énergie hv. Il existe donc pour la lumière 
une sorte de dualité onde-corpuscule : le photon est à la fois onde et corpuscule. 


11-3° Hypothèse de Louis de Broglie - Dualité onde-corpuscule (1924) 


D'après Louis de Broglie, la dualité onde-corpuscule observée pour les photons n'est 
pas propre seulement à la lumière, mais également à toutes les particules de la physique. Il 
stipule que : 


A toute particule de masse m en mouvement avec une vitesse v est associée 
une onde de longueur d'onde : 



mv 


Cette relation n'est qu'une généralisation de la relation Planck-Einstein, elle peut être 
retrouvée, dans le cas d'une onde lumineuse, à partir des relations suivantes : 

E = hv relation de Planck 

E = me 2 relation d'Einstein 


L'égalité de ces deux relations, avec v = c/ A, conduit à : 



me 


La relation de Louis de Broglie est générale à toutes les particules (photons, électrons, 
atomes, molécules, etc.), c'est-à-dire à tout objet microscopique ayant de faibles masses et 
de grandes vitesses. Toute particule en mouvement manifeste les deux aspects : l'aspect 
ondulatoire et l'aspect corpusculaire, il y a donc une sorte de dualité onde-corpuscule. 

L'énergie d'une onde de longueur d'onde A donnée par : 



montre que pour une particule donnée cette énergie est d'autant plus grande que la 
longueur d'onde associée est courte et inversement. Pour une énergie donnée, la longueur 
d'onde (A = h/mv ) est d'autant plus grande que la masse ou la vitesse de la particule est 
petite. 


Expérience des fentes de Young 
L'expérience des fentes de Young constitue une 
vérification expérimentale de l'hypothèse de Louis de 
Broglie dans le cas des électrons. Un faisceau d'électrons 
émis par une source est dirigé vers un écran en passant 
par deux fentes. On observe un motif d'interférence 
semblable aux franges d'interférence dans le cas de la 
lumière. L'observation de ces franges montre clairement l'aspect ondulatoire des électrons. 

Pour résumer, les électrons se comportent comme des corpuscules dans le 
phénomène photoélectrique et comme des ondes dans l'expérience des fentes de Young, ce 
qui montre l'aspect dual de l'électron. 
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11-4° Hypothèse de Heisenberg - Principes d'incertitude (1925) 

La mécanique classique permet de connaître à tout instant la position et la vitesse 
d'un objet mobile, c'est-à-dire connaître sa trajectoire. La mesure de l'une des deux 
grandeurs (position, vitesse) ne perturbe pas la mesure de l'autre. Ce principe en physique 
classique n'est valable que lorsqu'on a affaire à des objets macroscopiques. Cependant pour 
des objets microscopiques, ce principe cesse d'être valable. En effet, à l'échelle 
microscopique des expériences ont montré que la mesure d'une grandeur perturbe le 
système et entraîne à une incertitude à la mesure d'une autre grandeur. Ces imprécisions 
obéissent au principe d'incertitude de Heisenberg qui s'énonce de façon générale : 

On ne peut connaître simultanément et avec précision deux grandeurs 
physiques A et B incompatibles : 

AA AB > h avec h = h / 2k 

où AA et AB désignent respectivement les incertitudes sur les grandeurs A et B. Deux 
grandeurs physiques sont dites incompatibles lorsque la mesure de l'une perturbe la 
mesure de l'autre. D'après cette inégalité, on remarque que la mesure de A est d'autant 
plus précise (erreur AA petit) que B sera moins précise (erreur AB grand) et inversement. 
Une relation équivalente à la précédente est également proposée : 

8A 8B > - 
2 

Dans ce cas, les déviations 5A et 8B sont les écarts moyennes définis par : 

JôA = A-< A> 

[ôB = B-<B> 

ou les écarts quadratiques moyennes : 

JsA = V< A 2 >-< A> 2 
|ôA = V<B 2 > - < B > 2 


avec : 

< A > : valeur moyenne de la mesure de A 

< A > 2 : carré de < A > 

< A 2 > : valeur moyenne de la mesure de A 2 

Ces valeurs moyennes sont définies plus explicitement plus loin (voir Chapitre II). 

Les grandeurs physiques incompatibles qui obéissent aux inégalités de Heisenberg, 

sont : 

(position, impulsion) — » Ax Ap >ti 

(temps, énergie) — > At AE > h 

(angle, moment cinétique) — > Acp AL Z > h 

Devant l’impossibilité de connaître avec précision la position et l’impulsion d’une particule 
en mouvement, le principe d'incertitude interdit toute notion de trajectoire ordinaire. Par 
contre, une particule en mouvement devrait être vue comme un nuage et ce nuage est décrit 
par la notion de probabilité de présence de la particule. Cette probabilité est définie à partir 
d'un être mathématique T, appelé "fonction d’onde", solution de l’équation de 
Schrôdinger. 
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11-5° Hypothèse d'état quantique - Equation de Schrôdinger (1926) 

La notion d'onde associée à une particule a permis à Schrôdinger de développer en 
1926 sa célèbre équation dite équation de Schrôdinger. Pour déduire cette équation, on part 
de l'équation classique d'une onde se propageant dans l'espace avec une longueur d'onde X 
et une amplitude ®(x,y,z) : 


Une telle équation est une équation purement ondulatoire. Pour y introduire l'aspect 
corpusculaire, on fait appel à la relation de Louis de Broglie X = h/p : 


La quantité de mouvement p = mv peut s'exprimer en fonction de l'énergie totale E (somme 
des énergies cinétique et potentielle) : 

E = jmv 2 + V = ^— + V -> p = V2m(E-V) 

2 2m 

La substitution de p dans l'équation de propagation de l'onde conduit à l'équation : 

A® + — (E-V)® = 0 
h 

dite équation de Schrôdinger. Cette équation décrit le mouvement d'une particule de masse 
m soumise à un champ de potentiel V. Les solutions de cette équations à rechercher sont la 
fonction ®(x,y,z) et l'énergie E. La fonction ®(x,y,z) est appelée fonction d'onde qui décrit 
l'état quantique de la particule. 

L'équation de Schrôdinger peut se mettre sous la forme suivante : 


La quantité entre crochet représente une suite d'opérations mathématiques qu'on désigne 
par le symbole H appelé opérateur hamiltonien : 

n 2 

H — A+ V 

2m 

Par suite, l'équation de Schrôdinger prend la forme équivalente suivante : 

H® =E® 

appelée équation aux valeurs propres dont O est la fonction propre de l'opérateur H et E 
la valeur propre associée. 

11-6° Principe de correspondance 

La mécanique quantique se base sur la connaissance des opérateurs mathématiques 
qui gèrent les grandeurs physiques. Pour déduire la relation de correspondance entre 
grandeurs de la mécanique classique et opérateurs de la mécanique quantique, on part de la 
fonction d'onde classique associée à une onde monochromatique : 




A® + ^-® =0 
H 2 
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Cette fonction classique est de nature purement ondulatoire, on y introduit l'aspect 
corpusculaire via la relation de Louis de Broglie (À = h/p x ) et l'aspect quantique via la 
relation de Planck (v = h/E ) : 


T(x,t) = Ae^ 


(xpx-Et) 


Il est clair que si l'on veut obtenir les deux grandeurs classiques E (énergie) et p x 
(impulsion), il suffit de dériver T successivement par rapport à t et par rapport à x : 


= E¥(x,t), = p x T(x,t) 


D'où les deux relations de correspondance entre grandeurs en mécanique classique et leur 


correspondant en mécanique quantique 


Grandeur 

— >• opérateur 

E 

- ft =-M 

Px 

- p.=M 

Ainsi, à l'énergie classique E corresponde l'opérateur hamiltonien H et à la quantité de 

mouvement p x correspond l'opérateur 

impulsion p x . Les équations aux valeurs propres 

qui correspondent aux deux opérateurs s'écrivent : 

H * = 

1 ôt 



Dans les deux cas, la fonction d'onde doit dépendre des deux variables position x et temps t 
c'est-à-dire Tfxù). 


III- Les postulats de la mécanique quantique 

La mécanique quantique, comme toute nouvelle théorie, repose sur un certain 
nombre de postulats qui servent de principes généraux applicables à tout système 
quantique. 

111-1° Notion d'état quantique 

Postulat I : L'état quantique d'un système physique est complètement défini par une 
fonction d'onde 'F(q,t), fonction des coordonnées de position (q) et du temps (t) : 
T(q,t) = Etat 

La variable q est prise ici au sens large permettant de repérer la position spatiale de la 
particule, elle peut être des coordonnées cartésiennes (x,y,z), coordonnées sphériques (r,0,cp) 
ou tout autre jeu de variables. La connaissance de T permet d'accéder à toute information 
sur le système étudié tel que son énergie, son moment cinétique, et toute autre propriété. 

111-2° Principe de superposition des états 

Postulat II : Si un système physique peut se trouver dans des états décrits par les 
fonctions d'onde ®i, ®2, ..., ®n, il peut également se trouver dans un nouvel état 
décrit par une combinaison linéaire de celles-ci : 

¥= JcjOr 

i=l 

Ci étant les coefficients de participation des états individuels. Ce postulat est introduit afin 
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de satisfaire le phénomène d'interférence ou de recouvrement des ondes. 

111-3° Principe de correspondance généralisé 

Postulat III : A toute grandeur physique (A) en mécanique classique correspond un 
opérateur linéaire et hermitique ( A ) en mécanique quantique : 

A -> Â 

L'opérateur doit être linéaire pour satisfaire le principe de superposition des états (postulat 
II) et hermitique pour qu'à l'opérateur A correspond des valeurs propres réelles, c'est-à- 
dire mesurables ( voir Chapitre II). Par suite tout opérateur hermitique et linéaire est appelé 
opérateur "observable". Tout opérateur de la mécanique quantique s'obtient en appliquant 
les deux correspondances suivantes : 

q -> q=q p, -> P,=fJj 


En mécanique classique, l'impulsion p = mv est un vecteur à trois composantes (p x ,p y ,p z ) 
et les opérateurs quantiques correspondant s'écrivent : 


Px 

_ h ô 
i dx ' 

Py 

-ÏlJL f, -Tl d 

i Sy ' i dz 

Connaissant ces relations de base, tout autre opérateur de la mécanique quantique peut 
alors être construit : 

• impulsion 

P 

-> 

p = |-V p 2 = -h 2 A 

• énergie potentielle 

V(x,y,z) 


V = V(x,y, z) 

• énergie cinétique 

T = p 2 /2m 

— » 

T = p 2 / 2m = -h 2 A / 2m 

• énergie totale 

E = T + V 

— » 

H = T + V 

• moment cinétique (*) 

L = r Ap 

—> 

L = r Ap = yï aV, L 2 =-fc 2 (r 2 A-?V.Vr) 


(*) (a a b) • (c a d) = (ac)(bd) - (ad)(bc) 


111-4° Mesure des grandeurs physiques 


Postulat IV : Le résultat de la mesure d'une grandeur physique A dans un état 
quantique 4L ne peut être que l'une des valeurs propres a n de l'opérateur Â 
associé: 

Â4 / n =a n v P n 


En d'autres termes, les seules valeurs observables de 
la grandeur physique A sont les valeurs propres a n 
de l'opérateur Â associé. L'ensemble des valeurs a n 
forme ce que l'on appelle le spectre discret de Â ; la 
valeur ao la plus basse correspond à l'état 
fondamental de fonction d'onde 4L et les valeurs 
supérieures correspondent aux états excités. 


A 



•* états excités 


) L) ■* état fondamental 
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111-5° Evolution temporelle d'un système physique 

Postulat Y : L'évolution dans le temps d'un système quantique se trouvant dans un 
état H / (q,t) est régie par l'équation de Schrôdinger généralisée : 

HT 

i ôt 

où l'hamiltonien H est l'observable associé à l'énergie totale du système qui dépend 
explicitement du temps via la partie énergie potentielle : 

H=T + V(q,t) 

Les systèmes pour lesquels V dépend du temps, on cite par exemple les systèmes 
atomiques et moléculaires non isolés soumis à un champ extérieur électrique ou 
magnétique. De tels systèmes évolutifs au cours du temps sont dits systèmes 
non-conservatifs ou systèmes dynamiques. Au contraire, dans le cas des systèmes 
conservatifs, dits systèmes stationnaires, l'hamiltonien H est indépendant du temps : 

H = T + V(q) 


c'est l'exemple des atomes et molécules isolés. Dans ce dernier cas, H n'agit que par sa 
variable d'espace q, son action sur une fonction H / (q, t) est indifférente par rapport à la 
variable t; on peut donc séparer les variables dans l'équation de Schrôdinger. La fonction 
d'onde est alors représentée par le produit de deux fonctions indépendantes: fonction 
d'espace ®(q) et fonction de temps T(t): 

T(q,t) = 0(q)T(t) 


En substituant dans l'équation de Schrôdinger, on obtient: 

M>( q) dT(t) 
i dt~ 

ou bien 


T(t)H v F(q) = 


HO(q) _ i h dT(t) 
O(q) “ T(t) dt 


On retrouve deux équations différentielles à variables indépendantes; par conséquence 
chacun des deux membres de l'équation doit être égal à une constante et cette constante 
n'est autre que l'énergie E du système. L'équation dépendante du temps : 

1m . _æ dt 

T(t) t, 

a pour solution : 

T(t) = exp(-^t) 
n 


de telle manière que la fonction d'onde totale s'écrit : 

_æ 

^(q,t) = Ÿ(q)e 

La fonction 'L(q) indépendante du temps satisfait l'équation aux valeurs propres : 

H'F(q) = E'F(q) 

dite équation de Schrôdinger indépendante du temps dont ses solutions particulières 'F(q) 
décrivent les états stationnaires du système. 
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IV- Equation de Schrôdinger de quelques systèmes physiques 

En toute généralité, pour établir l'équation de Schrôdinger d'un système physique, 
on fait appel au principe de correspondance (postulat III). Selon ce principe, à toute 
grandeur physique en mécanique classique correspond un opérateur hermitique et linéaire 
en mécanique quantique. 


IV-1° Equation de Schrôdinger pour une particule 

Considérons une particule de masse m soumise à un potentiel scalaire V(x,y,z) et 
mobile dans un référentiel trirectangle Oxyz. A l'énergie totale E, somme de des énergies 
cinétique T et potentielle V, on fait correspondre l'opérateur hamiltonien H qui s'obtient de 
la manière suivante : 


E = T + V 


— (Px+PÇ+Pz) +V 
2m J 


— (Px +P^+Pz) +V 
2m 1 


Sachant que P p = y ^ ( q = x, y ou z ), l'hamiltonien devient : 


ou encore : 


h = 82 

2m \^ôx 2 ôy 2 ôz 2 


H = — —A + V 
2m 


V(x,y,z) 


L'équation de Schrôdinger s'obtient en faisant agir H sur une fonction d'onde 'P : 
h 2 

A + V l 'F = E'P 


ou encore : 


A'F + ^fE-VJ'F = 0 
h 2 


Une telle équation décrit les états stationnaires d'une particule de masse m mobile dans 
l'espace tridimensionnel Oxyz et soumise à un champ de potentiel V(x,y,z). Lorsque celle-ci 
est astreinte à se déplacer sur un axe, par exemple l'axe x, son équation de Schrôdinger 
devient : 


d 2v P 

dx 2 


2m 

h 2 


(E-V) V P = 0 


où l'énergie potentielle ne dépend dans ce cas que de la coordonnée axiale V(x). Dans tous 
les cas, l'équation de Schrôdinger est une équation différentielle du second ordre, son 
intégration nécessite la connaissance de la forme analytique de l'énergie potentielle V. Deux 
cas de figure peuvent être distingués : 
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(a) V = 0, c'est le cas d'une particule libre, son énergie totale est purement cinétique : 

% 2 % 2 d 2 T 

- — AT = ET ou _ = ^ 

2m 2m dx 2 

De telles équations décrivent les mouvements libres de translation ou de rotation de la 
particule. 

(b) VAO, c'est le cas d'une particule liée par un champ de potentiel issu de forces qui la 
retiennent: champ de force électrostatique ou champ de force de rappel. 

En toute généralité, un potentiel V traduit des 
interactions issues des forces F auxquelles la particule 
est soumise. Lorsqu'une force F est constamment 
dirigée vers un point fixe de l'espace, on dit que le 
champ de potentiel est de type central : 

V(x,y,z) = V(r) 


V" 


^ 0 

Le potentiel ne dépend que de la distance r qui sépare x ^ 

la particule à l'origine O. Si la force F = | F | qui décrit l'interaction dérive d'un potentiel, 
on écrit : 

F = -™ 
dr 

Ainsi, la connaissance de l'expression explicite de F permet la déduction de celle du 
potentiel V(r) par simple intégration entre les deux bornes r 0 et r : 


V(r) = - J^F dr 


la borne inférieure r 0 correspond à l'origine du potentiel et la borne extérieur r définit la 
position instantanée de la particule. En toute généralité, cette origine est choisie de telle 
manière que la force d'interaction F ou l'énergie d'interaction V soit nulle : 

r=r 0 -> F(r 0 )=0 V(r 0 ) = 0 

Cas de champs de forces électrostatiques 

Un champ de force électrostatique est un champ dans lequel deux charges 
électriques ponctuelles q et q 1 se repoussent ou s'attirent selon qu'elles sont de même signes 
ou de signes opposés. En accord avec la loi de Coulomb, la force F qui régit l'interaction des 
deux charges est inversement proportionnelle au carré de la distance r qui les séparent : 

F = lÆ q #q' 

r 2 r 

K étant un facteur de proportionnalité dont sa valeur dépend du système des unités 
employées : 

K = — en MKS 
47is 0 

K = 1 en CGS 
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En remplaçant cette force dans l'expression de l'énergie potentielle, on obtient : 


V(r) = -f r K^-dr = K^_ 
J0 ° r 2 r 


L'origine du potentiel, dans ce cas, correspond aux deux charges séparées à l'infini (F=0 si 
r— > go), c'est-à-dire une origine pour laquelle les deux charges sont sans interaction V(r)=0. 

Dans ces conditions, l'équation de Schrôdinger d'une particule de masse m se 
mouvant dans un champ de force électrostatique s'écrit : 


L'atome d'hydrogène en est un exemple où Tunique électron de charge négative (q = -e) 
gravite autour du noyau de charge positive (q'= +e), son équation de Schrôdinger dans le 
système des unités CGS s'écrit : 


La fonction d'onde à rechercher dans ce cas est tridimensionnelle s'exprimant en fonction 
des trois variables d'espace T^y/Z). Dans le cas d'un atome hydrogénoïde (système à un 
seul électron et un noyau de charge +Ze), l'équation de Schrôdinger s'écrit : 


Cas de champs de force de rappel 

Un champ de force de rappel est un champ dans lequel une masse m est 
constamment attirée vers un point fixe O par une force de rappel F. De tel système est à 
l'image d'une masse m reliée à une origine fixe O par un 
ressort de constante de raideur ou constante de force k, a a a £ 
dit aussi oscillateur. Sous l'effet de la force de rappel, le WWW 

système effectue des mouvements oscillatoires dits _ ^ J 

mouvements de vibration. La masse m subit dans ce cas 

des déplacements de quantité x = r - r 0 à partir de sa position d'équilibre r 0 avec r étant la 
position instantanée après allongement (r>r 0 ) ou rétrécissement (r<r 0 ) du ressort. Si 
l'oscillateur est harmonique, la force de rappel F est simplement proportionnelle au 
déplacement x en suivant la loi de Hook : 


Le signe négatif indique que la force F et le déplacement x s'opèrent dans des sens opposés. 
Comme la force de rappel dérive d'un potentiel, on a : 


En choisissant l'origine du potentiel au repos du système (F = 0 si x = 0), c'est-à-dire 
correspondant à la position d'équilibre (r = r 0 ), on aura : 





F = -kx 


V(x) = J*kx dx 
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Ainsi, l'équation de Schrôdinger de l'oscillateur harmonique s'écrit : 



W = 0 


La fonction d'onde T impliquée dans ce cas est une fonction monodimensionnelle 'L(x). Il 
existe également des cas où l'oscillateur est à plusieurs dimensions. Dans le cas par exemple 
d'un oscillateur qui se déplace dans le plan (x,y), la fonction d'onde est bidimensionnelle 
v F(x,y) solution de l'équation : 


Le terme k xy xy traduit l'interaction des deux modes d'oscillation. Parmi les systèmes 
physiques qui utilise le modèle de l'oscillateur harmonique à une ou deux dimensions, on 
peut citer le cas de molécules diatomiques AB (oscillateur à une dimension) et triatomiques 
ABC (oscillateur à deux dimensions). Dans les deux cas, les liaisons chimiques sont 
assimilées à des ressorts de constante de raideur k : 


IV-2° Equation de Schrôdinger pour un système à plusieurs particules 

Considérons un système physique formé de n particules de masses m 1 ,m 2 ,...,m n 
soumis à un champ de force dérivant d'un potentiel V. L'énergie potentielle doit s'exprimer 
en fonction des coordonnées de l'ensemble des particules VfxjyjZj, x 2 y 2 z 2 , ...,x n y n z n ) • 
L'hamiltonien d'un tel système s'écrit : 


Î-— A,+V'P = E'P 

b 2m i J 

Si toutes les particules ont la même masse = m (particules dites identiques), on aura : 


4 +V v F = E'P 

Li 2m J 

La fonction d'onde 'P dans ce cas doit être exprimée en fonction des 3n coordonnées du 
système : 


Pour simplifier l'écriture de la fonction d'onde, la notation v P(l,2,...,n) peut être introduite 
sachant que i désigne les trois coordonnées cartésiennes (xj, y^Zj) de la particule i. 







1 Sx 2 ôy 2 ôz 2 

Schrôdinger est : 


où A- =-@— + -^— + -^— est le Laplacien associé à la particule i de masse mi. L'équation de 

1 c^2 cs.,2 ^2 


v P(x 1 y 1 z 1 , x 2 y 2 z 2 , ...,x n y n z n ) s v P(l,2,...,n) 
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Parmi les systèmes physiques à particules identiques, on peut citer l'exemple des 
atomes. Un atome est un système formé d'un ensemble d'électrons de charges -e gravitant 
autour d'un noyau de charge +Ze. On est donc en présence d'un système électrostatique 
dont l'énergie potentielle est de type coulombien : 




Par exemple, l'atome d'hélium (Z=2) est un système à deux électrons en interaction avec un 
noyau de charge +2e. Dans ce cas, on définit deux types d'énergie potentielle dus à deux 
types d'interaction: énergie d'attraction électron-noyau (VNe) et énergie de répulsion 
électron-électron (V ee ) : 



Par suite, l'équation de Schrôdinger de l'atome s'écrit : 


h 2 

2m 


(Ai+A 2 )T 


2e 2 2e 2 

h r 2 


'¥=E'¥ 


La fonction d'onde du système est une fonction des deux variables électroniques : 

T(l,2) = 'P(x 1 y 1 z 1 ,x 2 y 2 z 2 ) 

En général, pour un atome comprenant n électrons gravitant autour d'un noyau de charge 
+Ze, les termes des énergies cinétiques et potentielles s'écrivent : 


T = Z-— Ai 

i 2m 


fc*I- 


V e e=I X — 

i<j r ü 



et l'équation de Schrôdinger correspondant s'écrit : 


2m r= 


+z Z— 

i < i r ü . 


Dans ce cas, la fonction d'onde est une fonction des coordonnées des n électrons : 
'P(b2,...,n)s'F(x 1 y 1 z 1 , x 2 y 2 z 2 , ...,x n y n z n ) 

Parmi les systèmes physiques à particules non-identiques, on peut citer le cas des 
molécules qui possèdent deux types de particules en mouvement: électrons qui se 
déplacent dans le champ des noyaux et noyaux qui se déplacent dans le champs des 
électrons. L'exemple de la molécule d'hydrogène peut bien illustrer les différents termes 
impliqués dans l'hamiltonien de Schrôdinger : 

H = T N + T e + V Ne + V ee + V NN 
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où T N et T e définissent les énergies cinétiques des noyaux et des électrons, V Ne , V ee et 
V NN les énergies potentielles d'attraction électrons-noyaux, répulsion électron-électron et 
répulsion noyau-noyau. Si la molécule comporte n électrons et N noyaux, les différents 
termes se généralisent comme suit : 


Tn ~Z-— — Aa./ 

a 2Ma 
n N 7 . P 2 

v Ne = S E-— 

i A TLA 


T e = S- — Ai 
i 2m e 

n n2 

Vee = Z I- 
i < j rij 



et la fonction d'onde dépendant des 3n variables électroniques et 3N variables nucléaires se 
note : 


T = 'F(x 1 y 1 z 1 , ...,x n y n z n , XjY^ , . . . , X N Y N Z N ) 
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CHapitre II 


La mécanique quantique, dans son contexte le plus général, est une théorie 
nécessitant une manipulation mathématique à la fois des fonctions d'ondes et des 
opérateurs. Elle incite donc à une revue d'ensemble d'outils mathématiques de base qui la 
constituent. 

I- Propriétés fondamentales des fonctions d'onde 
1-1° Notion de fonction d'onde 


La fonction d'ondeT en mécanique quantique est la représentation de l'état 
quantique d'un système physique dans la base de dimension infinie des positions (x,y,z) et 
du temps (t) : 

T = T (x,y,z, t) 

Si le système est dans un état stationnaire, la fonction d'onde s'écrit : 

_iE 

'Ffx, y, z, t) = 'Pfx, y, z) e n 


Sous cette forme v F(x,y,z,t) est une fonction complexe dépourvue de tout sens physique 
direct à cause de l'existence du nombre imaginaire i = V— T - Pour en attribuer un sens 
physique, on doit prendre le carré de son module : 

|q /| 2 = x ¥ *x ¥ 


afin d'éliminer la partie exponentielle portant le nombre imaginaire : 


|'F(x,y,z,t| 


T (x,y,z) e n v F(x,y,z) e 


'P (x,y,z)'F(x,y,z) 


Une telle quantité est appelée densité de probabilité de 
trouver la particule en un point M(x,y,z) de l'espace. La 
probabilité de trouver la particule dans un élément de 
volume dx=dxdydz entourant le point M est donnée par : 

dP = | Y | 2 dx 



dite probabilité élémentaire. Dans tout le domaine d'espace délimitant la position de la 
particule, cette probabilité est donnée par l'intégrale : 

P = J | v P| 2 dx = JJJ | l F| 2 dxdydz 


où P peut prendre des valeurs comprises entre 0 et 1 : 


0<P<1 


Ainsi, la condition : 


|P = 0 particule est nul part dans l' espace (W = 0) 

[P = 1 particule est quelque part dans l' espace (VP * 0) 

J |<P| 2 dx = l 
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forme ce qu'on appelle la condition de normalisation de la fonction d'onde et T dans ce cas 
est dite normée. Si l'intégrale donne une valeur différente de 1 : 

jl'FjdT = N 

on peut toujours la normaliser en cherchant une nouvelle fonction T' telle que : 

¥' = OF et £j v F'| 2 dT = 1 

En remplaçant T' dans la condition de normalisation, on trouve le facteur de normalisation 

C recherché 'F'= — . 

VN 


1-2° Conditions de bonne de fonction d'onde 


Pour qu'une fonction d'onde soit une solution acceptable de l'équation de 
Schrôdinger et représentant réellement un état quantique d'un système donné, elle doive 
satisfaire certaines conditions mathématiques dites conditions de bonne fonction d'onde : 


- continue : 

T ne doit strictement pas posséder de 
discontinuité en aucun point de l'espace, 
pour qu'il n'y ait pas de solutions 
singulières de l'équation de Schrôdinger. 

continue 

- uniforme : 

T ne possède qu'une seule valeur en un 
seul point de l'espace pour qu'il y ait une 



seule probabilité en un seul point. 

uniforme 

-finie : 

T ne doit posséder aucune valeur infinie 
en aucun point de l'espace pour qu'il y ait 

V 


toujours une probabilité finie. 



v p 



discontinue 


'ÿ 



non uniforme 
¥ 



finie (x— »oo) 


non finie (x— »co) 


- carré sommable , l'intégrale j| v F| di doit être finie et convergente afin qu'il y ait une 


- dérivable : 


probabilité certaine. 

¥ doit être infiniment dérivable : (d n 'P/dx I1 ) n _ >Q0 finies. 


Dans l'espace de définition de la particule, il existe des points particuliers pour lesquels la 
fonction d'onde T doit satisfaire certaines conditions appelées conditions aux bords ou 
conditions aux limites. Soit x e [a, b] est le domaine délimitant la position de la particule, 
dans ce cas la fonction d'onde doit avoir des valeurs T (a) et T (b) finies. 


1-3° Espace d'états et notation de Dirac 

Les fonctions d'onde en mécanique quantique ont pour propriété mathématique 
fondamentale d'appartenir à l'espace d'Hilbert direct 6 appelé espace d'états. Il est formé 
de fonctions d'onde { v F 1 , v F 2 ,..., v P n ,...} de dimension infinie (n— »oo) car les variables de 
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position dans chaque fonctions 'Pi(x,y,z) prennent une infinité de valeurs 
(-00 < x,y,z < +oo ). A l'espace direct S correspond son dual S* caractérisant les fonctions 
complexes conjuguées { 'h, , v F 2 ,..., v F n ,...} : 

'Pi e S -> 'F* e g’ 

L'espace d'états est un espace vectoriel hilbertien semblable à l'espace vectoriel euclidien. 
Une fonction d'onde 'Pj de l'espace direct S est symbolisée par le vecteur ket | ^P; ) , alors 

que sa conjuguée 'Pj dans l'espace dual S* est représentée par le vecteur bra | : 

|'P i )eS ('Pii e S’ 

Par définition, le ket est l'adjoint du bra et inversement le bra est l'adjoint du ket : 

oj'H'ir 

L'opération adjoint ( symbolisée par le signe +) correspond simplement à l'opération 
transposée-conjuguée. 

Etant donné que l'espace des états est un espace vectoriel, on peut toujours se servir 
d'une base de fonctions d'onde connue {cp 1 ,cp 2 ,...,(p n } pour développer une fonction d'onde 
ou vecteur du même espace : 

¥ = 1 ^ -> 'P*=Sc*cp* 

ou : 


|'P) = £c i |(p i ) ('P| = £c*(cp i | 


Ainsi, un vecteur ket et son bra peuvent être représentés, dans la base {cpi} considérée, sous 
la forme matricielle suivante : 


M = (ki>.|<p„ 


(*1 = (c*. 4 


(cp)(c) 


(c) + (cp) + 


Propriétés : 

• multiplication par une constante X : 

• dérivation : 

• linéarité : 

• antilinéarité : 


|W) = ?,|'P), (W|=X*(Y| 

|4'p 1 +4'p 2 )=4|'p 1 )+4|>P2> 
(x.i'p, + 4^ | = 4(>Pi 1+4(41 
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1-4° Produit scalaire 

Le produit scalaire de deux fonctions distinctes est défini par : 

= J^'Pjdx 

où x Pi(x,y,z) et v Fj(x,y / z) s'expriment en fonction des mêmes variables, l'intégration s'étend à 
tout l'espace de définition des variables et dx=dxdydz est un élément de volume entourant 
le point de coordonnées (x,y,z). Le produit scalaire en mécanique quantique n'est pas 
commutatif contrairement au cas du produit scalaire des vecteurs euclidiens ( VjVj = VjVj ) : 

(WW * WW 

Dans le premier membre, le produit est un scalaire alors que dans le second c'est une 
matrice ( voir plus bas). 

Condition d'orthonormalisation : 

Si le produit scalaire de deux fonctions d'onde distinctes est nul, celles-ci sont dites 
orthogonales : 

= 0 i*j 

Si le produit scalaire d'une fonction d'onde par lui-même est égal à l'unité, celle-ci est dite 
normée : 

d'ih) = i i-î 

Ces deux situations peuvent être rassemblées en une seule, appelée condition 
d'orthonormalisation : 

WW = ** 

Ôij est le symbole de Kronecker valant 0 si i ^ j et 1 si i = j. 

Décomposition d'une fonction d'onde : 

Soit {cp 1 ,cp 2 ,...,(p n } une base orthonormée c'est-à-dire ^(p i |(Pj) = 5 i j. Toute fonction T 
appartenant au même espace peut être définie dans cette base sous la forme d'un 
développement en une série suivant : 

v L = £c i cp i 

La décomposition de T dans la base { cpj } entraîne que les coefficients Cj sont uniques. On 
peut montrer que cette décomposition résulte simplement de la projection de ¥ dans la 
base cpj . Pour cela, multiplions le développement de T par cp* et intégrons : 

(cpj l'k) = ^<Pj I^CjCPi^) = Zc^cpj |cp^ 

Si la base est orthonormée, le résultat de la projection T selon la direction cpj de la base n'est 
autre que le coefficient de participation de cpj c'est-à-dire : 

(WW = Çcftj = Cj 
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Représentation matricielle d'un produit scalaire : 

Soit {cpi, 92, . . ., cp n } une base de fonctions d'onde servant à développer une fonction 'P 
du même espace : 

v P = Zc i cp i « |'P> = (cp)(c) 

Dans la notation fonctionnelle ou vectorielle, le produit scalaire de T sur lui-même peut 
être représenté comme suit : 


(^|^> = EEc*c j (cp i |cp j } 


r('Pl'P) =(c) + (cp) + (cp)(c) 


-LL CjCjSjj I I = (c) (S)(c) 

i j J v 

On définit alors une nouvelle matrice (S) = (cp) + (cp) dont ses éléments sont Sÿ = ^cp; |cp J, 
appelée matrice de recouvrement des fonctions de base dont sa dimension est d'ordre n x n : 
(Su S 12 ... S ln 
S„ ... S~ 


(s)- 


V ni °n2 


••• S m 


Dans l'hypothèse d'une base orthonormale, Sij = ôÿ, la matrice de recouvrement est une 
matrice unité (S) = 1 : 

r \ 0 ... 0" 

0 1 ... 0 


(S) = 


0 0 


LcjCj = Lcjf 


et le produit scalaire se réduit à : 

('F| X F) = (c) + (c) = CjCj + C* 2 C 2 + ••• + C*C n 

Dans ce cas, la condition de normalisation de T s'écrit : 

\ N 2 =i 

Le produit bra ® ket n'est pas commutatif ('F | V F) ^ | V F)( V P|; en effet, dans l'hypothèse d'une 
base orthonormale on a : 

<¥]¥) = (c) + (c) et |'P)('P| = (cp)(c)(c) + (cp) + 
le premier produit est un scalaire alors que le second est une matrice carrée. En effet : 


C.C, +C 0 C 9 +... + CC 
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C,C, C,C r. 


c i c n 

C 9 C n 


ce, c c~ 


On retrouve donc que le produit de (c) par (c) + n'est pas commutatif: (c) + (c) * (c)(c) + . 


Propriétés : 

• conjugaison : (w i 1 ^2 ) = (^2 | ^1 ) 

• linéarité: (^ 1^1 +X 2 V 2 ) = À 1 (o|'F 1 ) + X 2 (o|'F 2 ) 

• antilinéarité : +^ 2 |o) = |®) + X.*(Y 2 |o) 

• positivité de la norme : || O || = ^ (O |o) > 0 

La norme d'une fonction d'onde est toujours positif || O || > 0 , elle est identiquement nulle 
seulement si la fonction d'onde est nulle 0 = 0 et si la fonction d'onde est nulle la particule 
n'existe nul par dans l'espace. 


II- Opérateurs 

11-1° Notion d'opérateur 

En mécanique quantique, un opérateur est une opération ou suite d'opérations 
mathématiques transformant une fonction d'onde d'un espace S donné en une autre 
fonction du même espace : 

fe 5 - A - > ÂW = 'P ' e 6 

Exemples : ô/ôx, ô 2 /ôx z , sin, cos, yj~ , .... Les opérateurs observables qui représentent les 
grandeurs physiques ont pour propriété mathématique fondamentale d'être linéaires et 
hermitiques : linéaires pour satisfaire le principe de superposition des états et hermitiques 
pour que leurs valeurs propres soient réelles ( voir plus bas 111-1° et 111-2°). 


Propriétés : 

• multiplication par une constante : 

• somme : 

• produit : 

• associativité : 

• distributivité : 

• action sur un ket : 

• action sur un bra : 


Â(XW) = k(ÂV) 

(Â + B)^ = Â'¥ + (B V P) 

(ÂB) Ÿ = Â(B'E) 

(ÂB) CT = Â(BCT) 

(Â + B)Cf = (ÂC V F) + (BC'P) 

â| v e) = |â v f) 

(Â v f| = ( v P|Â + ou ('P|â=^â + 'f| 
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11-2° Notion de commutateur 


En général, le produit de deux opérateurs Â et B dépend de l'ordre dans lequel les 
opérations produits ÂB et BÂ sont écrites. L'opération ÂBT signifie d'abord action de B 
sur T et ensuite action de A sur le résultat BT . Comme le produit AB n'a aucune raison 
d'être égal au produit BÂ, on définit leur différence par un commutateur noté : 

[Â,B ] = ÂB - BÂ 

Si [Â,B ] * 0, les deux opérateurs ne commutent pas, par contre si [Â,B] = 0 les deux 
opérateurs commutent. 


Exemples : 



[x,y] = 0 


Pour l'évaluation de tout commutateur, il faut agir une fonction d'onde quelconque. Ainsi, 
le premier commutateur se démontre de la manière suivante : 




T + x 


ST 5T 

x 

ôx ôx 


T 


On voit dans la liste des exemples ci-dessus que les commutateurs qui opèrent sur des 
variables différentes commutent toujours. En toute généralité, deux opérateurs à variables 
indépendantes commutent : 

[Â(x),B(y)]=0 


Propriétés : 


[Â,B ] 

- -[B,Â] 

[XÂ,B] 

= [ Â , À,B ] - X[B,Â] 

[Â,B+C] 

= [Â,B]+[Â,C] 

[Â,B C ] 

= [Â,B]C + B [Â,C] 

[Â,Â n ] 

- 0 

[Â,B n ] 

= IBj[Â,B]B n -j- 1 
j=o 

[Â,B]T 

= Â (B T) - B (ÂT) * [ÂT,BT] 


Les commutateurs jouent un rôle important en mécanique quantique puisqu'ils sont en 
relation directe avec le principe d'incertitude de Heisenberg. 


Commutateurs et principe d'incertitude : 

On peut facilement montrer que les opérateurs associés aux variables 
complémentaires position-impulsion (x, p x ) et temps-énergie (t, E) ne commutent pas : 

[x,p x ] = i/z, [t,H]=-i h 

Par ailleurs, le principe d'incertitude de Heisenberg concernant les grandeurs physiques 
correspondant nous enseigne que : 

Ax Àp x >fi At AE >h 

Ces deux inégalités peuvent être reliées aux commutateurs correspondant si l'on considère 
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la mesure des grandeurs dans un état T supposé normé : 

Ax Ap x >||('P|[x,p x ]|'î / )|| = ft 
AtAE > | ('Fl [t. H] l'E) [J = fi 

Pour démontrer, par exemple la première inégalité, il suffit de calculer d'abord la valeur 
moyenne de l'opérateur associé au commutateur [x,p x ] = \h : 

(^1 [x,p x ] l'F) = (Y| m |ÿ) = i^('P|'P) = ih 


ensuite calculer le module de la quantité obtenue : 

| <'F|[x,p x ]|'F)|| = ||i h | = V(iA)*(i*) 


= h 


On peut généraliser la relation d'incertitude de Heisenberg à tout couple de grandeurs 
physiques A et B : 

AA AB > | ('F | [Â,B] l 'F) J « AA AB > || J^*[Â,B] v Pdx || 


Deux propriétés importantes peuvent en découler : 


Propriété 1 : 

La condition pour que deux grandeurs physiques A et B soient 
simultanément mesurables avec précision est que les opérateurs Â et B qui 
les décrivent commutent : 

AA = AB = 0 [â,b]= 0 A et B compatibles 

Dans ce cas, on dit que les deux grandeurs A et B sont compatibles à la mesure. 


Propriété 2 : 

Si deux opérateurs commutent, ils partagent un système commun de 
fonctions propres. 

Si par exemple HL est fonction propre de Â , et que Â commute avec B , alors la même 
fonction HL est fonction propre de B : 

ÂT 0 =aT 0 , BŸ 0 »b>F 0 

Propriété 3 : 

Deux grandeurs physiques A et B sont dits incompatibles à la mesure si 
les opérateurs A et B qui les décrivent ne commutent pas : 

AA.AB>/? [â,b]* 0 A et B incompatibles 

Dans ce cas, les fonctions propres de Â ne les sont pas nécessairement pour B . 


11-3° Valeurs propres et fonctions propres d'un opérateur 

Lorsqu'un opérateur Â agissant sur une fonction T redonne la même fonction T 
multipliée par une constante a, 

Â x P = a x P 

on dit que T est fonction propre de A pour la valeur propre a. Par exemple, la fonction e™ 
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est fonction propre de l'opérateur d/dx pour la valeur propre m. En effet : 


T-fe™) = me™ 
dx 


Propriété 1 : fonction propre 

Une fonction propre d'un opérateur est définie toujours à une constante 
multiplicative près. 

Soit 'F une fonction propre d'un opérateur Â, on peut montrer que la fonction X F' = W est 
également fonction propre de A avec la même valeur propre : 

Â'F = a'F 

Â'P' = Â(w) = A,(âÿ) = À,(a x F) = a(w) = a^' 

On peut donc dire que les deux fonctions T et XT décrivent le même état quantique, elles 
différent seulement par la norme. Pour décrire un état quantique par une fonction d'onde 
unique, on peut toujours la normaliser (valeur de X unique) : 

'F'=W tel que ( V F , | V F') = 1 

c'est-à-dire : 



Propriété 2 : opérateurs puissances 

Si T est fonction d'un opérateur Â linéaire, elle est également fonction propre de 
l'opérateurs puissance Â n quelque soit le nombre entier n : 

ÂW = a^ -A Â nx P = a nv F 

Dans le cas par exemple de n=2, la relation est facilement vérifiable : 

Â 2v F = Â(Â'F) = Â(a'F) = a(ÂT) = a(a'F) = a 2x F 


Propriété 3 : spectre discret 

Lorsqu'un opérateur A possède des valeurs propres toutes différentes quelques 
soient les fonctions propres, on dit que A possède un spectre discret de valeurs propres, 
c'est-à-dire non dégénéré : 


ÂTj = 

=a 2 v î , 2 


A^ = a n Ÿ n 




Â v P i =a i v P i , i = 1, 2, ..., n 


Dans ce cas, à chaque valeur propre ai est associée une seule et unique fonction propre v Fi. 


- 28 - 


Mécanique Quantique S MC (S4)~ Tès 


CHapitre II 


Propriété 4 : dégénérescence 

Lorsqu'à une valeur propre a d'un opérateur Â est associée m fonctions propres 
distinctes (¥1, ¥ 2 , . . f m ), on dit qu'il y a dégénérescence d'ordre m : 


Â¥, = a ¥, 
A¥ 2 = a ¥ 2 


Â¥ v 


= aT m 


avec a 


¥, 

^2 

¥ m 


On rencontre, par exemple, des états doublement dégénérés, triplement dégénérés, etc. Un 
état est doublement dégénéré, si à une valeur propre est associée deux fonctions propres 
distinctes, et un état est dit triplement dégénéré si à une valeur propre est associée trois 
fonctions distinctes : 

états doublement dégénérés: 

=a T, 

AV 2 = a ¥ 2 



états triplement dégénérés: 
Â¥, =a T, 

Â¥ 2 = a ¥ 2 
Â¥ 3 = a ¥ 3 


TT ¥ 2 ¥3 


Propriété 5 : combinaison linéaire 

Toute combinaison linéaire des fonctions propres dégénérées d'un 
opérateur Â est aussi fonction propre de Â pour la même valeur propre. 

Soient TT et ¥2 deux fonctions dégénérées d'un opérateur Â : 

Â¥j = a T, et Â¥ 2 = a ¥ 2 


Pour la combinaison linéaire ¥ = A.i¥i + >^¥2, on a : 

Â(à 1 ¥ 1 +à 2 ¥ 2 ) = X 1 (Â¥ 1 )+à 2 (Â¥ 2 ) = A, 1 (a¥ 1 )+À 2 (a¥ 2 ) = a^^ +À 2 ¥ 2 ) 

On voit que la combinaison linéaire est également fonction propre de l'opérateur pour la 
même valeur propre. 

11-4° Valeur moyenne d'un opérateur 

On doit distinguer deux cas de figure : ¥ est fonction propre d'un opérateur Â et ¥ 
non fonction propre de cet opérateur. Dans le premier cas, le résultat de mesure de 
l'observable Â est la valeur propre associée. Cette valeur propre peut être déterminée de 
deux manières selon que l'on opère dans l'espace fonctionnel ou dans l'espace vectoriel des 
états. La correspondance entre les deux écritures est la suivante : 


- 29 - 


Mécanique Quantique S MC (S 4 )~ Tès 


CHapitre II 


espace fonctionnel 

espace vectoriel 

Â v F = a v F 

« Â|'F) = a|'F) 

x ¥*ÂV = a x ¥*'¥ 


j'P*Â'Pdx = aJ'P*Wx 

« (T Â T) = a(T 'P) 

_ j’ v F*Â v Pdx 

<'F|À|y> 

j'PVdx 

(T V) 


Par contre, si T est une fonction d'onde quelconque, c'est-à-dire non fonction propre de 
l'opérateur, on détermine alors une valeur moyenne. La valeur moyenne d'un opérateur A 
est définie par la quantité : 

£ 

P 


IM 

?l*> 


Pour pouvoir calculer cette valeur moyenne, on doit connaître la fonction d'onde T. Pour 
cela, on peut toujours faire appel à une base connue de fonctions d'onde {cpi, 92, . . ., cp n } du 
même espace selon laquelle on peut développer T sous forme de la combinaison linéaire 
suivante : 




Dans ce cas, la valeur moyenne devient : 


(A> = - 



ziqcj^i Âcpj) 
i j 

ÜcîcjAÿ 
i j 

(tw 

H 

ZScJc j (q) i |<p j ) 
i j 

Elc-CjSÿ 
i j 

avec A i j=(cp i Â cpj 

} et s >j = 

=<«pil«pj>- 




Si le jeu de base {cpi, 92, . . ., cp n } est orthonormé, on a = 8^ , et le terme au 
dénominateur s'écrit : 

n n n n n n 2 

XZqCjSÿ = XXqCjÔij = £ C * C . =X| Ci | 

i j i j i i 

La valeur moyenne devient : 

ZIc-CjAÿ 

W- IJ . , 

Dans le cas où la fonction d'onde T est normée, c'est-à-dire £|cj|“ = 1, la valeur moyenne se 
réduit à : 

(A) = (>F|Â|>F) = iîc’cjA, 
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Propriété de dérivation : 

Si l'opérateur Â dépend par exemple de la variable x, la dérivation de sa valeur 
moyenne est : 

Ô ( A ) * (dA\ 
dx \dx I 


En effet, le membre de gauche a pour expression : 


Sx 


(ff l À W + MflÙ + Hlr) 


alors que le membre de droite est simplement : 


£-<W> 


Dépendance temporelle d'une valeur moyenne : théorème d'Ehrenfest 

Considérons une grandeur physique A représentée par un opérateur Â pouvant 
dépendre explicitement du temps et dont sa mesure dans un état quantique T donne la 
valeur moyenne suivante : 

(A) = <T|Â|>P) 

La dérivée temporelle donne : 

La variation temporelle de la fonction d'onde T peut être déduite de l'équation de 
Schrôdinger généralisée : 

espâpe direct .-ësfëàœ, t ^àl;’ 

Hh) = "''1 * (hH + « 

|ht) = ih\^pj (T|H = H hermitique 

Le remplacement de telle quantité dans (A) donne : 


D'où : 


^ = (Ùffi'Ù - iî(' l '|HÂ|'P> + Î C('F|ÂH|'F> 

= (' ï 'l#l' 1 '} + ïLLÂH-HÂl'î') = + jL('p|[Â,h]| , p) 


d( A ) = M\ 

Ôt \dt/ 


[Â,H]) 


De telle variation temporelle d'une grandeur physique A porte le nom du théorème 
d'Ehrenfest. Les deux dérivations S(A)/St et (SA /St) ne sont égales que si les deux 
opérateurs Â et H commutent [Â,H] = 0. 
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11-5° Représentation matricielle d'un opérateur 

Rappelons que dans une base donnée de fonctions d'onde {cpi, cp 2 , . 
l'F) et son bra ('F] peuvent être représentés matriciellement par : 


, cp n }, tout ket 


| v P> = Ec i |cp i > 

(¥| = £c*(cpi| 

(cp) = (|cpi),|cp 2 ), •••,!%>) 


|'F> = (cp)(c) 
('F| = (c) + (cp) + 


(c) + = 



(M) 


V 

((p) + = 


(C) = 

c 2 


>n l 


v C n, 


La valeur moyenne d'un opérateur Â pris dans un état T supposé normé s'écrit alors : 

(A) = ('F|Â|'F> = (c) + (cp) + Â (cp) (c) = (c) + (A) (c) 
où (A) = (cp) + Â (cp) est la matrice représentative de l'opérateur Â dans la base (cp) dont ses 
éléments sont A- = ^cp ; |Â|cp^ . En effet le produit d'un vecteur colonne par un vecteur ligne: 



>if 

(Âcp 1 ),Âcp 2 ),...,Acp n )) 



(A) = 

(<p 2 | 





.Kl 




donne : 






"((PlIÂ^) (CP! |X| Cp 2 ) ... <CPi |Â| CPh}^ 


^A]] A 12 ... A ln N 

(A) = 

(*P2|Â|cpi) (tP2|Â|cp 2 ) ••• (<P 2 | Â |cPn) 

- 

A 2 i A 22 ... A 2n 


v (<Pn|Â|cPl) (<Pn|Â|cp 2 ) ••• (<Pn|Â|cp n ) y 


v A nl A n2 ... A nny 

On peut de même rechercher une représentation matricielle de l'équation au* 


Â|¥) = a|'F) « Â(cp)(c) = a(cp)(c) 

En multipliant à gauche par l'adjoint (cp) + , 

(cp) + Â (cp) (c) = a(cp) + (cp) (c) 
on obtient alors : (A) (c) = a (S) (c) 


Due à la présence de la matrice de recouvrement (S) , cette équation n'est pas une équation 
aux valeurs propres. Pour qu'elle en soit ainsi, il faut choisir une base orthonormée pour 
ramener (S) à une matrice unité (S) = (1) . Dans ce cas, on obtient une véritable équation 
aux valeurs propres : 

(A) (c) = a (c) 
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L'équation matricielle (A) (c) = a (c) traduit en fait un système de n équations homogènes et 
linéaires : 

AijCi + A12C2 h i-Aj n c n = ac^ 

A21C1 + A22C2 "i 1 A2 n c n = ^2 


AjjjCj + A n 2 c 2 + • • • + A nn c n = ac n 

que l'on peut écrire encore : 

(An -a) c i + Ai2 c 2 -i h A ln c n = 0 
A 2 jCi + (A 2 2 — a)c 2 h i-A2 n c n =0 


A nl Cj +A n 2 c 2 h — + (A nn -a)c n =0 

Ce système admet des solutions non nulles seulement si le déterminant correspondant est 
nul : 

(A n -a)cj + A 12 c 2 h 1- A ln c n 

A 2 1 c 1 +(A 22 -a ) c 2 +'" + A 2n c n _ n 


AjjjCj +A n 2 c 2 H — + (A nn -a)c n 

Après développement de ce déterminant, on obtient une équation de degré n en "a" dont sa 
résolution fournira les n valeurs propres possibles (a 1 ,a 2 ,...,a n ) de la matrice (A). A 
chaque valeur propre ai, on doit associer un vecteur propre (c 1 ,c 2 ,...,c n ) dont les 
coefficients s'obtiennent en résolvant le système d'équations correspondant. 

III- Classes d'opérateurs en mécanique quantique 

111-1° Opérateur linéaire 

Un opérateur Â est dit linéaire si pour tout couple de fonctions Ti et ¥j, il vérifie 
l'égalité suivante : 

ÂfÀi'Fi + X j'Fj ) = (ÂTj ) + X j (Â'Fj ) 

En mécanique quantique, les opérateurs observables qui représentent des grandeurs 
physiques, tels que l'impulsion p, le moment cinétique L, l'hamiltonien H, ..., sont 
linéaires. Cette linéarité est imposée afin que ces opérateurs satisfassent le principe de 
superposition des états (postulat II). 

111-2° Opérateur adjoint 

Un opérateur Â admet son adjoint A + , si pour tout couple de fonctions T; et Tj, il 
vérifie l'égalité suivante : 

( v E i |Â + |T j ) = (T j |Â|'E i )* « j x P i * Â +v FjdT = J 'U Â* 'P* dx 

Cette égalité peut être retrouvée facilement : 

(T, lÂ+l'Fj) = (Â'Fj l'Fj) = (’FjlÂ'Fi)* = (Tj |Â|Tif 
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Comme les deux intégrales définissent les éléments de la matrice adjointe (A) + de la 
matrice (A) et comme l'opération adjoint est l'opération transposée-conjuguée, on a : 



L'opérateur adjoint a une grande importance en mécanique quantique, il sert à définir les 
vecteurs bra de l'espace dual S* . En effet, si | V F) est un ket propre de l'opérateur Â dans 
l'espace direct 5, 

1^) — Â|'F) = |'F') e 6 
alors dans l'espace dual S* , on aura : 

('Fl - A - - > ( V F|Â + =( v L'| g s " 


Propriétés : 

|â v f} = â|'f) 
(â'f| = ('f|â + 
('f|â = ^â +v f| 
(Â + ) + =Â 
(ÀÂ) + = à*â + 


(Â + B) + = Â + +B + 
(Â + iB) + =Â + -iB + 
(ÂB) + = B + Â + 

[â,b Mb\â + ] 


111-3° Opérateur hermitique 

Un opérateur Â est dit hermitique si pour tout couple de fonctions 'Fi et 'Fj, il vérifie 
l'égalité suivante : 

|Â| 'Fj ) = (Tj |Â| 'Fj)* « j 'F* Â'Fjdr = j'F j Â* W*dx 

Autrement dit, un opérateur est hermitique s'il est égal à son propre adjoint (hermitique = 
auto-adjoint) : 

Â = Â + 


Dans la formulation matricielle, une matrice est hermitique si les éléments diagonaux sont 
réels, et les extra-diagonaux sont transposés conjugués l'un de l'autre : 


Exemple : 



(A) = 


La matrice n'est est hermitique, c'est-à-dire (A) = (A) + , que si a = 1 - i . 
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En mécanique quantique, on impose que les opérateurs observables représentant les 
grandeurs physiques (impulsion p, moment cinétique L, hamiltonien H, ...) d'être 
hermitiques afin que leurs valeurs propres soient réelles, c'est-à-dire soient mesurables. 

Propriétés : 

• les valeurs propres d'un opérateur hermitique sont réelles et les vecteurs 
propres sont mutuellement orthogonaux, 

• les fonctions propres d'un opérateur hermitique forment une base 
orthonormée complète, 

• Si deux opérateurs hermitiques commutent, ils admettent un système 
commun de fonctions propres, et inversement si deux opérateurs 
hermitiques admettent un système commun de fonctions propres ils 
commutent, 

• Si deux opérateurs Â et B sont hermitiques, leurs transformés Â + iB et 
A - iB sont adjoints l'un de l'autre. 

Exemple : 

Montrons que l'opérateur p x = -i hdl ôx est hermitique sachant que x e ]-oo,+oo[ . 

La propriété d'hermiticité de p x impose que : 

HT T n j dx 

Intégrons par partie le membre de gauche : 

CT Ÿ n ih J^mCbc = -i^['P*'P m t + y ^n dx 

Comme les fonctions Tn et T m doivent s'annuler à l'infini (finitude des fonctions d'onde), 
alors le premier crochet est nul. D'où : 

[ +00x P n 1 -i/z— Vmdx = i*r> m -<dx = f + ^ m i,AVdx = f +co W |-i/jAw ) dx 
J - œ dxj m J - œ m dx J -“ m dx J - œ dx n J 

et la propriété d'hermiticité est démontrée. 

IV- Propriétés de parité et applications 


IV-1 Notion de parité 

La parité est une propriété de symétrie qui consiste à inverser le signe des 
coordonnées d'espace d'un point M : 


(x.y.z) (— x,— y,— z) 

ou 

r -> -r 

avec r = OM et | r | = r 
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En coordonnées sphériques, l'opération d'inversion est équivalente à une rotation d'angle n 
autour d'un axe suivie d'une réflexion par rapport au plan perpendiculaire à l'axe de 
rotation : 

r — > r 0 — > 71 — 0 cp — I 7i + cp 

En mécanique classique, il découle de cette propriété le principe de conservation de parité 
qui s'énonce : 

"les lois physiques ne se modifient pas lors de l'opération d'inversion 
des coordonnées d'espace d'un point matériel par rapport à l'origine" 

Ainsi, les grandeurs physiques, comme l'énergie totale, moment cinétique, impulsion, ... 
conservent leurs valeurs ou leurs normes après la transformation subie. 

En mécanique quantique, la parité est d'une importance capitale puisqu'elle permet 
d'étudier le comportement des fonctions d'onde lors de l'opération d'inversion et par suite 
elle importe des importantes simplifications des intégrales. Cette opération d'inversion en 
mécanique quantique est décrite par un opérateur P unitaire tel que : 

P 2 = 1 


appelé opérateur parité. II consiste à transformer une fonction d'onde ¥(?) en une autre 
fonction x P(-r) de coordonnées renversées : 

PF(r) = v P(-r) 


Comme P est supposé unitaire, ses valeurs propres sont simplement +1 ou -1 : 

P'P(r) = t'Pfr) => 'P(-r) = +'P(r) 

Les états pour lesquels 'P(-r) = +'P(r) sont appelés états pairs, et ceux pour lesquels 
'P(-r) = -'P(r) sont des états impairs. 

Exemples: fl'fx) = sinx impaire A^x) = xsinx paire 

2 2 
'P(x) = e paire A'(x) = xe impaire 

Pour le calcul des intégrales, cette propriété de parité apporte de grandes simplifications 
puisqu'elle permet à n'étudier les fonctions que sur la moitié de leur intervalle de 
définition, l'autre moitié étant déduite par symétrie. En effet, si [-a,+a] est l'intervalle 
servant pour l'étude d'une fonction f(x) donnée, deux cas se présentent : 

- si f (x) est pair, on a : 

j^*f(x)dx = J° a f(x)dx + j^ a f(x)dx = 2j^ a f(x)dx 

- si f (x) est impair, on a : 

f>) dx = f° a f ( x ) dx + J 0 +a f( x ) dx = ~ j 0 +af W dx + j 0 +af W dx = 0 


Dans le cas d'une intégrale portée sur le produit de deux fonctions fi(x) et f 2 (x), on a : 


j;; fl (x)f 2 (x)dx “ 


fi( x )f 2 ( x ) impair 
fi( x )f 2 ( x ) pair 


Ces résultats peuvent être appliqués au cas des intégrales ou éléments de matrice en 
mécanique quantique selon le corollaire suivant : 
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Corollaire 1 : 

" L'intégrale ('F; |Â| 'Fj ^ = J 'F* Â'Fjdx est non nulle seulement si le produit 
v F*Â v Fj est pair par inversion du signe des coordonnées d'espace" 

Si l'opérateur Â est invariant sous l'effet de l'opération d'inversion (par exemple 
A(x) = A(-x) ), la valeur de l'intégrale suit le corollaire 2 : 

Corollaire 2 : 

" L'intégrale ^ |Â| 'Fj ''j est non nulle seulement si le produit v F* v F j 
est pair par inversion du signe des coordonnées de position" 


IV-2° Applications 


Cas d'une particule oscillante : 


Considérons une particule de masse m réalisant des oscillations harmoniques 
suivant l'axe x. Son hamiltonien est : 


H = T + V(x) 


h 2 d 2 
2m dx 2 


2 


L'origine O (x=0) de l'axe x'Ox constitue le centre d'inversion du système. On voit que les 
deux termes de H, énergie cinétique T et potentiel V, se conservent lors de l'opération 
d'inversion : 


T (x) = T (-x) V(x) = V(-x) 


car elles impliquent des quantités x 2 . Par suite, l'hamiltonien est invariant H(x) = H(-x) 
sous l'opération d'inversion et ses fonctions propres l F i (x) sont soit paires soit impaires. 
Considérons les trois premières fonctions propres de H : 

^ =Nj e - ax2/2 

*F 2 =N 2 xe " ax2/2 
'F 3 =N 3 (2a x 2 -1) e ~ ax2 /2 


On voit que 'Fi et T 3 sont pairs alors que ¥2 est impair. Si nous souhaitons calculer les 
éléments de matrice LL = (w i |h| 'Fj 'j de l'hamiltonien dans la base {'Fi/HL, ¥ 3 }, on écrit : 

f H 11 H 12 H 13 ' 

(H) = H 12 H 22 H 23 

H„ H„ 


Des importantes simplifications peuvent être apportées en faisant appel au corollaire 2 : 
h 12 = (^î |H| *F 2 ) , T1T2 impair ^ H 12 = 0 

H 13 = ('Fi |H| W 3 ), pair H 13 * 0 

H 23 = (*F 2 |H| ^3 ) , ^ 2^3 impair H 23 = 0 
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(H) = 


H n 0 H 

0 H 22 0 

«13 0 H 3: 


13 


On peut aussi utiliser les propriétés de parité pour le calcul explicite des éléments de la 
matrice (H) . Par exemple, pour l'élément Hn, on a : 

H„ = <'P 1 |H|'P 1 ) - + = -£<'î'.l^l' i '3) + |('i' 1 |x 2 |-P 1 ) 

Après dérivation et substitution, l'élément matriciel prend l'expression suivante : 


/T a 




k tC a‘ 


ttl* 2 !*.) 


Ainsi, on peut utiliser la propriété de parité des fonctions pour écrire : 
( V P 1 | V F 1 ) - N?J*V ax2 dx - 2N?J 0 +c V ax2 dx 

('P 1 |x 2 |'i' 1 ) - Nf£^x 2 e _ax2 dx - 2Nf| 0 +c °xV ax2 dx - 

La même démarche peut être appliquée aux autres éléments H 13 , H 22 et H 33 . 


-fJ! 

N? [V 
2 


Cas d'une particule dans un champ de potentiel central : 


Considérons une particule de masse m se mouvant dans un champ de potentiel de 
type central, l'électron de l'atome d'hydrogène en est un exemple. L'hamiltonien de tel 
système est : 


H = - 


h 2 _ 
2m 


A + V(r) 


Comme V(r) est à symétrie sphérique (r 2 = x 2 +y 2 +z 2 ), l'hamiltonien est invariant sous 
l'opération de symétrie. Les fonctions propres v F n ^ m (r,0,cp) de H s'écrivent comme le 
produit des fonctions radiales R ni ?(r) et angulaires Y |m (0,cp) (voir chapitres suivants) : 

'P nto (r,0,cp) = R n |(r)Y to (0,cp) 

où n, i et m sont les nombres quantiques principal, orbital et magnétique des orbitales 
atomiques de l'atome d'hydrogène. 

Comme l'opération d'inversion transforme les coordonnées (r, 0, cp) en (r, tü- 0, Ti+cp), 
la fonction radiale R nJ? (r) est invariante, alors que la fonction angulaire Y to (0,tp) peut 
changer de signe. On peut montrer que l'application de l'opérateur parité P à v P ntol (r,9,cp) 
donne : 

=<-!)' » T n , m (r,e,<p) = (-l)''P n(m (r,n-e,7 t + <p) 

La parité des états propres de l'atome d'hydrogène est donc (-l) f , c'est-à-dire si I pair 
'F nto est pair, et inversement si i impair v F ném est impair. 
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Sachant que i décrit la symétrie des orbitales atomiques : 

i 0 1 2 3 4 5 ... 

symétrie s p d f g h 

alors les états atomiques s, d, g sont pairs et les états p, f, h sont impairs. D'après les 
expressions des premières orbitales atomiques ls, 2s et 2p de l'atome d'hydrogène : 

ls = N ls e~ r/a ° 

2s = N 2s (2a 0 -r) e _r/2a ° 

2p z = N 2p r e _r/2a ° cos0 

on voit que lors de l'opération d'inversion seul le terme cos0 change de signe, alors que 
r = (x 2 + y 2 + z 2 ) 1/2 demeure toujours positive. En effet, 

cos(7i-0) = cos7icos0-sin7rsin0 = -cos0 
et l'orbitale 2p z est impaire. 
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Chapitre III 


Le moment cinétique est une grandeur qui caractérise le mouvement angulaire des 
particules dans l'espace, il intervient principalement dans la résolution de l'équation de 
Schrôdinger de l'atome d'hydrogène et du rotateur diatomique ainsi que dans la 
caractérisation de la symétrie des états électroniques des atomes et molécules diatomiques. 
On distingue dans la physique des particules trois classes de moments cinétiques: moment 
cinétique orbital L, moment cinétique de spin S et moment cinétique total J. 


I- Moment cinétique orbital 
1-1° Définition 

En mécanique classique à un point matériel se mouvant avec une vitesse v autour 
d'un point fixe O est associé un moment cinétique L donné par : 

L = r a mv = r a p 

où r est le rayon vecteur définissant la position de la 
masse m relative au centre O et p le vecteur impulsion. 

Dans le repère trirectangle Oxyz portant les vecteurs 
unitaires i, j, k selon les trois axes, les vecteurs r , p et L 
se définissent par : 

r = xi + yj + zk, p = p x i +p y j + p z k, L = L x i + L y j + L z k 

Compte tenu de la définition du produit vectoriel, le moment cinétique s'écrit sous la forme 
d'un déterminant : 

ï j k 

L = x y z 
Px Py Pz 

Après développement de ce déterminant, on obtient les trois composantes classiques de ce 
moment : 

L x =yPz- z P y 

Ly = zp x -xp 
L z = x Py-yPx 


A 



Notons qu'en partant, par exemple, de l'expression de la première composante L x , les deux 
autres peuvent s'obtenir simplement par permutation circulaire des indices x, y et z. Le 
module du vecteur moment cinétique | L | se définit à partir de son carré : 


ou à partir de : 


L 2 =L 2 X +L 2 y +L 2 z 
| L | = | r | | p | sina 
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a étant l'angle que fait le vecteur position r avec l'impulsion p . Si une particule est 
astreinte à se mouvoir sur un cercle de rayon r fixe, on a a=90° et le moment cinétique 
devient simplement L = rp = mvr. 

Pour passer aux expressions quantiques, on applique le principe de correspondance 
en substituant p q par son l'opérateur associé p q = yj^-, ce qui donne : 


L x 


h 




Ly 

L z 



avec L 2 =L 2 +L 2 +L 2 


En physique, on choisit souvent l'axe Oz comme axe privilégié par rapport auquel les deux 
directions Ox et Oy jouent des rôles équivalents. Ainsi, le jeu des trois 
composantes L x ,L y ,L z peut être remplacé de façon équivalente par le jeu des opérateurs 
L + ,L-,Ê Z tels que : 

L + = L x +iL y 
L_= L x — iL y 
L z = L z 


L + et L_ sont appelés opérateur escalateurs. Etant donné que L x et L y sont hermitiques, 
les deux opérateurs L + et L_ sont adjoint l'un de l'autre : 


(L + ) + =L_ (L_) + =L + 


1-2° Relations de commutation 


En faisant usage aux définitions ci-dessus des composantes du moment cinétique, on 
peut démontrer facilement les relations de commutation suivantes : 


[L x ,L y ] = ïhh z 
[L y ,L z ] = iftL x 
[L Z ,L X ] = ihL 


O 


Ainsi, les trois composantes du moment cinétique ne commutent pas. En conséquence, elles 
ne sont pas mesurables simultanément et ne possèdent pas de système commun de 
fonctions propres. En revanche, chacune de ces composantes commute avec le carré du 
moment cinétique : 

[ L 2 , L q ] = 0 q = x, y, ou z 

Il est donc possible de mesurer simultanément le carré du moment cinétique et l'une de ses 
trois composantes. Cela signifie qu'ils admettent un système commun de fonctions propres. 
En particulier, si T est fonction propre de L 2 , il l'est également pour L z : 


L 2v F = a*? 
L z 'F = b v E 
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Avec les opérateurs escalateurs L + et L_, on peut retrouver facilement les relations de 
commutation suivantes : 

[L + ,L X ] = ~[L_,L X ] = C z 
[L+,L y ] = [L_,L y ] = ifcL z 
[L + ,L Z ] = -ftL + , [L_,L Z ] =+«L_ 

[L + ,L_] = C z/ [L + ,L 2 ] = [L_,L 2 ]=0 

Le carré du moment cinétique L 2 est défini soit à partir des trois composantes L x , L y , L z : 
L 2 = L 2 + L 2 + L 2 

soit à partir du jeu d'opérateurs L + , L_ , L z : 

L 2 = L + L_+L z (l z -/?) 
l 2 =l_l + +l z (l z + h) 

La première relation, par exemple, se démontre de la manière suivante : 

L + L_= (L x +iL y )L x -iL y = L 2 +L 2 -i(L x L y -L y L x ) = L 2 - L 2 -i(i^L z ) = L 2 -L Z (L Z +h) 


1-3° Expressions en coordonnées sphériques 


Un point M de l'espace peut être repéré soit par ses coordonnées cartésiennes 
(x,y,z ) soit par ses coordonnées sphériques (r,0,cp). La relation entre ces deux systèmes de 
coordonnées est : 


x = rsinBcoscp 
y = rsin0sincp avec 
z = rcos0 


0 < r < oo 

0 < 0 < 71 
0 < cp < 271 



L'élément de volume dx = dxdydz entourant un point M a pour expression en coordonnées 
sphériques : 

dx = r 2 sin0 drd0dcp 


où on y est identifie l'expression de l'angle solide dS = sin 0d0dcp , angle sous lequel un solide 
est repéré. 




- 42 - 


Mécanique Quantique S MC (S4)~ Tès 


Chapitre III 


Pour pouvoir exprimer les composantes L x , L y et L z en coordonnées sphériques, on doit 
d'abord établir les dérivées partielles d/dq (avec q=x,y,z) en fonction de ces coordonnées. 
Pour cela, on utilise les relations courantes de changement de variables : 

_d_fd l\JL 


ôq 


dqjôr ^dq J dG ^ dq J dtp 

("M 


dx 


“V 

dz J 


q = x, y, z 
Q = r, 0, q 


D'abord partant de la relation r 2 = x 2 + y 2 + z 2 et calculons sa dérivée : 
2rdr = 2xdx + 2ydy + 2zdz 

dr = (f') dX+ (r) dy + (f) dZ 

ensuite comparant à la différentielle totale : 


■(IMS 


nous en déduisons : 


-I dy+ I&l dz 

= sin G coscp 
= sin G sin cp 


= — = cos G 


De même en partant de la relation cosG = z/r, nous en déduisons : 


dG _ 

xz 

cos G coscp 

dx 

r 3 sin G 

r 

dG _ 

yz 

cos G sin cp 

dy 

r 3 sin G 

r 

dG _ 

z 2 - r 2 

sin G 

t dz 

r 3 sin G 

r 


Enfin à partir de la relation tgcp = y/x, nous obtenons : 
y cos 2 cp 


smcp 

rsinG 


cos cp 


coscp 

rsinG 


En substituant les relations ainsi obtenues dans les expressions d/d q (avec q=x,y,z), on 
obtient : 
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'5 

( sin 0 coscp)— + 

( cos0coscp^| 5 i 

r sin cp ) 

5x 

ôr 

l r J 50 1 

^rsin0 ) 

5 

(sin 0 sincp)— + ( 

^cos0sincp^ 5 t | 

cos^ 

5y 

ôr \ 

, r J 50 + 1 

^rsin0 ) 


— = (cose)— - 

5z V 1 lr J 00 


En substituant ces dérivées partielles dans les expressions des composantes du moment 
cinétique. 


O 

_z A] 


L - H (x 8 

-y-1 

^ ôz 

5yJ r 

y “il 0x dzj' 

X il 5y 

5x J 


on obtient les expressions en coordonnées sphériques que nous recherchons : 

r J . 5 n 5 

L x = m\ sincp — + cotg 0 coscp— | 

50 


ôcp 


L y -«(- 


8 „ . 3 I 

L v =m\ -coscp — + cotg 0 sin(p — 
50 5 cp J 


z i 5 cp 

De même, on peut exprimer les opérateurs L + , L_ et L 2 dans ce même système de 
coordonnées : 


L_ =he 1<p + icotg 0 — 

^ 50 ôcp 

T 2 +2 ( 1 d ■ n d 1 5 2 

sin 0 50 50 sin 2 0 5 cp 2 

Remarquons que l'opérateur L 2 peut s'exprimer en fonction de L z : 

f2 n 2 5 . _ 5 L 2 
sin 0 50 50 sin 2 0 

On voit que la dépendance en cp de L 2 est entièrement contenue dans le terme L 2 . 

Faisons remarquer que les opérateurs du moment cinétique orbital L x , L y , L z , L + , 
L_ et L 2 s'expriment uniquement en fonction des variables angulaires 0 et cp, ils sont 
indépendants de la variable radiale r, et c'est pour cette raison que le moment cinétique 
orbital est appelé aussi "moment angulaire". 

1-4° Relation avec l'hamiltonien de Schrodinger 

Considérons une particule soumise à un potentiel de type central V (r) et se mouvant 
dans l'espace tridimensionnel Oxyz . L'opérateur hamiltonien correspondant est : 

H=-L-A + V(r) 

2m 
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où A = + -^ 5 - + ~~z~ est le Laplacien exprimé en coordonnées cartésiennes. Le fait que le 

Sx 2 5y 2 5z 2 r r n 

potentiel V(r) est de type sphérique, on est amené à exprimer A en coordonnées sphériques : 


, 1 5 2 d 1 

A = — — r — + — 
r 2 5r 5r r 2 


1 5 . _ 5 1 

sin 0 — + — 

sin 0 50 50 sin 2 0 c 


La partie entre parenthèse peut s'exprimer en fonction de l'opérateur L 2 : 


a = J-A r 2A _ 

r 2 5r 5r h z r 


L 

2 2 


Ainsi, l'hamiltonien prend la forme 
H 


h 15 2 5 L 1 , x 

T — r — + + V(r) 

2m r 2 5r 5r 2mr 2 


On voit que la dépendance angulaire 0 et cp dans H est entièrement contenue dans le 
terme L 2 . 

On peut montrer facilement les relations de commutation suivantes : 

[H, L 2 ] = 0 
[H, L z ] = 0 
[L 2 ,L z ] = 0 

Les trois opérateurs H, L 2 et L z commutent deux à deux, ils possèdent donc un système 
commun de fonctions propres. Soit v P(r,0,cp) une fonction commune, on peut écrire : 

H'LCrAcp) =E'F(r,0,(p) 

L 2 v F(r,0,cp) = a x F(r,0,cp) 

L z v L(r,0,cp) = b x P(r,0,(p) 

Comme L 2 n'agit que sur les variables angulaires (0,cp), on peut procéder à la séparation 
des variables en posant : 

'L(r, 0, cp) = R(r)Y(0,cp) 

tel que : 

HRY = ERY 
L 2 Y = aY 
L Z Y = bY 

De même, comme L z n'agit que sur la variable cp, on cherche aussi des solutions de la 
forme: 

Y(0,cp) = F(0)<D(cp) 

tel que : 

L 2 F® = aF® 

L z ® = b® 

Finalement, la fonction d'onde d'une particule dans un potentiel central est formée du 
produit de trois fonctions d'onde indépendantes : 

v F(r, 0, cp) = R(r)Y(0,cp) = R(r)F(0)® (cp) 
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Pour la normalisation de telle fonction, on peut procéder à la séparation des intégrales 
comme suit : 

('P|'P) = J'PVdx = J J J W" Vr 2 sinGdrdGdcp = |R*Rr 2 drf J Y*Ysin0dedcp = jR*Rr 2 drf Y*YdS 
D'où : 

[(R|R) = jR*R r 2 dr 
( V P|'F) = (R|R)(Y|Y) avec J 

|^( Y | Y ) = J J Y * Y sin0d0 d(p 

Ainsi, lorsqu'il s'agit de calculer le produit scalaire, il est nécessaire de tenir compte 
correctement des éléments d'intégration r 2 dr pour la partie radiale R(r) et dS = sin0d0 pour 
la partie angulaire Y (0, cp) . L'intégrale relative à la partie angulaire peut aussi se séparer en 
produit des intégrales relatives aux variables 0 et cp : 

(Y | Y) = J J Y*Y sin0 d0d(p = J F 2 (0) sin0 d0 J<F*d> dtp 


Exemple: 

Normalisons la fonction radiale et la fonction angulaire suivantes : 
R(r) = Ne _r/a ° , Y(0,cp) = N’cos0 

sachant que J”x n e -ax dx = . 


(R|R) = j®R 2 r 2 dr = N 2 {“r V 2r/a °dr = N z | 

(Y | Y) = 0 2 *Y 2 sin 0d0dcp = N' 2 J 27t cos 2 0 sin 0d0dcp = J 271 dtpj^cos 2 0 sin 0d0 


f 2! ; 

= N 2 K) = 1 -> N = /U 3 

l(2/a 0 )% 

t 4 J UoJ 


= N’ 2 [cpïH-(cos 3 0)/3f = N’ 2 (2k) (2/3) = 


-> N'= I — 1 


Rotateur rigide : 

Un rotateur rigide est un système physique dont une 
masse m est libre d'effectuer des mouvements de rotation par 
rapport à un point fixe O de l'espace (rotateur spatial) sans 
aucune déformation de sa distance à l'origine (r fixe). Comme 
la particule est libre de se mouvoir, son énergie potentielle 
V(r) est nulle. Dans ces conditions, l'hamiltonien de 
\ & S L 2 

Schrôdinger H = — — r 2 — H — + Y(r) se réduit à: 

6 2m r 2 dr ôr 2mr 2 




où I = mr 2 est le moment d'inertie. Les fonctions propres communes aux deux opérateurs H 
et L 2 s'obtiennent par la résolution de l'équation aux valeurs propres suivante : 
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HY (0, cp) = EY (0, cp) 

c'est-à-dire : 

L 2 Y(0, cp) = 2IEY(0, cp) avec L 2 Y(0, cp) = aY(0, cp) 

Connaissant la valeur propre a de l'opérateur L 2 , on peut donc déduire les énergies propres 
du rotateur rigide E = . 

1-5° Valeurs propres des opérateurs L 2 et L z 

Il s'agit de déterminer les valeurs a et b associées respectivement à L 2 et L z : 

L 2 Y = aY 
L Z Y = bY 

Etant donné que le moment cinétique a la dimension de la constante h , on peut toujours 
choisir a et b sous les formes suivantes : 

<x = h\ 
b = tira 

tel que k et m sont des nombres sans dimensions et nécessairement réelles car ils 
correspondent à des opérateurs hermétiques. D'où : 

L 2 Y = h 2 k Y 
L Z Y = hm Y 

Pour rechercher les nombres k et m, on parte des deux relations suivantes impliquant les 
opérateurs escalateurs : 

L + L _= L 2 - L z (L z - h) 

L _L + = L 2 - L z (L z + h) 

En supposant que la fonction d'onde Y est normée, les valeurs moyennes correspondant 
aux deux opérateurs produits s'écrivent : 

(Y |L + L _| Y) = h 2 [k - m(m - 1)] 

(' Y |L _L+ 1 Y) = h 2 [k - m(m + 1)] 

Compte tenu de la propriété adjoint d'un ket (^Â = /Â +V F | , on peut écrire : 

(È _ Y |È _ | Y) = fi 2 [k - m(m - 1)] 

(l + Y |l + | Y) = h 2 [k - m(m + 1)] 

Les deux quantités correspondent aux normes des fonctions L + Y et L_Y qui doivent être 
nécessairement des quantités positives ou nulles. Il vient : 
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k - m(m - 1) > 0 

<=> 

k-m(m + l) > 0 


m“ - m-k < 0 
m 2 +m-k < 0 


On obtient alors deux inéquations de second degré en m dont leurs racines doivent obéir 
aux deux inégalités suivantes : 

i-VÂ ^ ^ i+VÂ 

< m < 

2 2 

-î-VÂ , -î+VÂ 

< m < 

2 2 

où A=l+4k est le discriminant commun aux deux inéquations, il doit être un nombre réel et 
positif pour que le nombre m soit réel. Si l'on pose l = (-1 + VÂ) /2, c'est-à-dire k = £(£ + !), 
il vient : 

-l < m < £ + 1 

-I - 1 < m < l 


Vérifiant graphiquement ces deux inégalités : 


-£- 1 -t 


On voit que les deux inégalités ne sont simultanément vérifiées que si le nombre m satisfait 
la condition : 

-I < m < +£ 

Il s'en suit que m prend 21 + 1 valeurs possibles. Finalement les deux équations aux valeurs 
propres s'écrivent : 

L 2 Y = h 2 £(£ + ])Y 
L Z Y = frmY 

Nous verrons plus bas que le nombre m doit être un entier positif, négatif ou nul, et le 
nombre £ doit être strictement positif ou nul. 


1-6° Fonctions propres des opérateurs L 2 et L z 

Pour trouver les fonctions propres Y(0,cp), on est amené à résoudre simultanément 
les deux équations aux valeurs propres relatives aux opérateurs L z et L 2 : 


ndY * v 
= hmY 

i dtp 

t2 1 d . A a L 2 

sin0 50 50 sin 2 0 


Y = a^ + 1)Y 


Après substitution de L 2 Y par h 2 m 2 Y et simplifications, nous obtenons : 


- 48 - 


Mécanique Quantique S MC (S4)~ Tès 


Chapitre III 


ÔY 

dcp 


= imY 


— -f 

sin 0 30 V 


n d Y 


l{l + 1 ) 


sin 2 ©J 


Y = 0 


Etant donné que les nombres I et m sont justes des paramètres dans ces équations, les 
fonctions Y(0,cp) doivent être recherchées pour chaque jeu de ces paramètres, elles 
porteront donc les labelles I et m : 

Y />m (0,<p) - |l,m) 


De telles fonctions sont appelées les harmoniques sphériques. En procédant à la séparation 
des variables en posant : 

Y,, m (0,<p) = F^ m (0) 0> m (cp) 

on obtient : 


dO m 

dcp 


1 d 
sin 0 d0 



£(£ + l) 



= 0 


Solutions de l'équation en cp : 

Les solutions de la première équation sont immédiates, elles sont de la forme : 
O m (cp) = N e im<p 

La normalisation de cette fonction permet de déduire le facteur N : 


= N 2 71 (e — imtp )(e imtp ) dcp = N 2 }^ dcp = 27iN 2 = 1 

1 


d’mfcp) = ~f= e imip 
Q2k 


0 < (D < 271 


Dans l'intervalle [0, 27 t], on peut vérifier que la fonction O m est continue, finie, carré 
sommable et dérivable mais pour qu'elle soit uniforme il faut que : 


c'est-à-dire : 


® m (cp) = ® m (cp + 27t) 


e i27tm _ cos (27 t m) + i sin(27nn) =1 — » 


Jcos(27im) = 1 
[sin(27im) = 0 


Ainsi, cette condition mathématique impose que le nombre m doit être un entier positif, 
négatif ou nul : 

m = 0, ±1, ±2, . .., ±£ 


Par suite, la valeur propre de l'opérateur L z est quantifiée par le nombre m appelé dans le 
langage d'atomistique "nombre quantique magnétique". On peut montrer que les fonctions 
propres ® m (cp) forment un système orthonormé complet : 


( 0 m| 0 m ') = Ô mm , 
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Comme à chaque valeur propre m est associée une seule et unique fonction propre ® m , on 
dit que l'opérateur L z est non dégénéré. Par ailleurs, comme on a m = 0, ±1, ±2, ±£ le 

nombre £ doit être un entier strictement positif ou nul : 

i = | m | = 0, 1, 2, ... 

appelé nombre quantique orbital (ou nombre quantique secondaire). Ainsi, la valeur propre 
de l'opérateur L 2 est quantifiée par le nombre i . Pour une valeur de l donnée, le nombre 
quantique m prend 2f +1 valeurs possibles et la quantité 2£+l est appelée multiplicité 
orbitale. 


Solutions de l'équation en 6 : 

Les fonctions F i m (0) s'obtiennent à partir de la seconde équation ci-dessus : 


1 d f . 

si 

sin 0 d0 ^ 


d0 


£(£ + !) - - 


Effectuons un changement de variable en posant = cos 0 : 


d 9 dF« m . 

— (1 — ^ 2 ) — — — + \ £(£ + l) - 
dÇ s dÇ 1 


l-^ 2 


On obtient alors une équation différentielle de second ordre et à coefficients variables. Dans 
la théorie des équations différentielles, l'intégration d'une telle équation possède des 
solutions de la forme : 

F^m(Q = N,, m P| m| (Q 

où N^ m est le facteur de normalisation et pj™^) le polynôme associé de Legendre ayant 
les expressions suivantes : 

(2£ + 1) £-\m\ 

2 (£ + |m|) ! 

i ( ni+IH M A e+\M 


Finalement, les fonctions Y f . m fonctions propres à la fois de L 2 et L z , s'écrivent : 

Y^ m (0,(p) = i^ p| m l(cos0)e im(p 
QIti 

Les harmoniques sphériques forment un système de fonctions propres orthonormées : 

Les expressions explicites des premières harmoniques sphériques sont données dans le 
tableau suivant : 
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l 

m 

Y^ m (0,<p) 

Orbitale 

atomique 

0 

0 

Y — 1 

Y °’°“V4^ 

s 

1 

0 

y io=Æ cos0 

Pz 

1 

+1 

Y m = +jJsin0e ±i(p 

(Px Py) 

2 

0 

il 

?R 

O 

O 

<X> 

1 

d z 2 

2 

±1 

Y 2)± i = sin 0 cos 0 e ±ltp 

(d xz ,d yz ) 

2 

+2 

Y — / 15 n P ±2icp 

Y 2- ±1_ V32^ 06 

(d X y’ d x 2_ y 2 ) 


Remarquons que les harmoniques sphériques Y e m (0,cp) = F e m (0) cD m (cp) sont réelles 
lorsque m=0 et complexes lorsque m^O; ceci découle de la forme des fonctions O m (cp) . En 
pratique, les harmoniques sphériques Y t m (0, cp) décrivent la partie angulaire des orbitales 
atomiques d'un atome monoélectronique x V n ç m = R n ^(r)Y to (0,cp) de type ns, np ou nd : 

n : nombre quantique principal (n = 1,2, . . ., œ) 

i : nombre quantique orbital (0 < f < n - 1 ) 

m : nombre quantique magnétique ( - 1 < m < +£ ) 

Partons de l'expression générale des harmoniques sphériques : 

Y 4m (0,cp) = N Am pj m l (cos0)e imip 

on peut montrer que l'application de l'opérateur parité P à une harmonique sphérique 
PY^ m (0,cp) = Y^ m (7t-0, 7i + cp) conduit à : 

PY,, m =(-l)% m 

La parité des fonctions Y^ m (0,cp) est (— 1)^ , c'est-à-dire : 

Y f>m (tc - 0, n + cp) = (- 1) ■ 1 Y t%m (0, cp) 

Si I pairY f m est pair, et si I impair Y^ m est impair. Ces fonctions sont en général d'ordre 
I en x/r, y/r et/ou z/r. Par exemple, pour les quatre premières harmoniques sphériques 
on a : 


Y 0 ,o 

il 

V4ït 

_ï : ^ i 

V4n 

= V P 1 

Y l,0 

=Æ cose 

= /X£ 

V 4n r 

= ^2 

Yi,i 


= [T x+i y 

\8n r 

= ^3 

y m 

' = i/ï sinee_i<p 

= /~3~ x ~iy 

V 8n; r 

= ^4 


On voit que si Y 0 0 restent inchangé lors de l'inversion (fonction paire) alors Y 10 , Y,] et 
Y 1 _ 1 changent de signes (fonctions impaires) car x, y, z changent de signe. 
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Exemple d'application : 

Supposons que les quatre premières harmoniques sphériques Y 00 , Y 10 , Yy et Y, forme 
une base dénommée { q>j , q> 2 , q> 3 , q> 4 } pour exprimer une matrice (A) associée à un opérateur 
Â supposé invariant sous toute opération de symétrie. Dans ce cas, les éléments de matrice 
entre fonctions de parité différente sont nuis et la matrice se réduit à: 

'A„ 0 0 0 N 

^ _ 0 A 22 A 23 A 24 

0 a 32 a 33 a 34 

v 0 A 42 A 43 A 44 y 

1-7° Action des opérateurs escalateurs 
sur les harmoniques sphériques 

Il s'agit d'établir les deux relations suivantes qui montrent l'action des opérateurs 
escalateurs L + et L _ sur une harmonique sphérique Y lm : 

L + Y*, m = hj£{£ + l)-m(m + l) Y />m+1 m < l 

L_ Y i m = hyl£{e + l)-m(m-l) Y^ m _, m>-£ 

Ces équations ne sont pas des équations aux valeurs propres, on voit que Y £m se 
transforme en Y^ m+1 sous l'action de L + et en Y i m _ x sous l'action de L_. 

Pour pouvoir démontrer les deux dites équations, on doit d'abord montrer que 
L + Y fm et L _Y £m sont fonctions propres de l'opérateur L z . Pour cela, on part du 
commutateur : 

[ L z , L+ ] = ± hh ± -> L Z L ± = L+L z ± 

Appliquons l'opérateur L Z L + = L + L z + hL + à une fonction Y^ m : 

L Z (L + Y^ m ) = L + (L z Y />m ) + h( t + Y £m ) = L + hmY im + hh + Y lm 

d'où : 

L Z (L + Y^ m ) =h( m + l)(L + Y 4m ) 

La fonction L + Y im est donc fonction propre de l'opérateur L z avec la valeur 
propre h(m + 1) . Par ailleurs, on sait qu'une harmonique sphérique Y t m+1 est aussi fonction 
propre de L z avec la même valeur propre h(m + 1) : 

L z Y/, m+ i =Hm+l)Y £m+l 
En comparant à l'équation aux valeurs propres, 

L z (L+Y />m ) =*( rn+l)(L + Y fm ) 

on voit que les deux fonctions L + Y^ m et Y f m+1 possèdent la même valeur propre h ( m + 1) , 
or comme l'opérateur L z est un opérateur non dégénéré, on a nécessairement : 

L+Y />m =Y £m+l ou L+Y^ m = c + Y^ m+1 


- 52 - 


Mécanique Quantique S MC (S4)~ Tès 


Chapitre III 


où c + est une constante à déterminer. Pour cela, calculons le produit scalaire suivant en 
supposant que les fonctions L + Y />m sont normées : 

(L + Y^ m |£ + Y Am ) = c 2 + (L + Y tm |L + Y^ >m ) = c+ 

On a par ailleurs : 

(L + Y tm |L + Y^ m ) = /z 2 [k-m(m + l)] = h 2 [l(l+Y)-m(m+ï)\ 

Par identification, on obtient : 

c + = hyj MJ. + 1) - m(m + 1) 

et la relation est démontrée : 

L + Y 4tn =^^ + l)-m(m + 1 )Y, ;m+ i 

De la même manière, on démontre la seconde relation concernant l'opérateur L : 

L-Y^, m =^V^ + l)“ m (m-l)Y 4m _ 1 

Les quantités à l'intérieur de la racine doivent être positives ou nulles, par suite la première 
relation n'est applicable que si m < i et la seconde si m > -£ : 

L + Y tm = hyl £(£ + 1) - m(m + l)Y <>m+1 si m < £ 

L_Y lm = fi ^ 1<JL + 1) - m(m - 1^^ si m > -£ 

On remarque que l'action des opérateurs escalateurs L + et L_ sur une harmonique 
sphérique Y^ m font respectivement monter et descendre le nombre quantique m 
d'un pas d'unité. D'où l'appellation de L + par l'opérateur de " montée " et L_ par 
l'opérateur de "descente". Ces deux relations sont très utiles pour la construction d'un jeu 
complet des harmoniques sphériques associées à une valeur de I donnée mais de m 
différents. Connaissant par exemple l'expression de Y f m , on obtient les expressions des 
deux jeux de fonctions suivantes : 



par action respective des opérateurs L + et L_. Notons que pour les harmoniques 
sphériques Y^ +i et Y e _ e qui correspondent aux valeurs extrêmes m = +£ et m = -l,ona: 

L + Y^ £ = 0 et L_Y^ _£ = 0 

mais : 

L_Y £ j = hyfïï Y L £_ x et L + Y^ = fijïi Y( _( + \ 

Afin de distinguer le nombre quantique m de la masse d'une particule m, on notera par la 
suite le nombre magnétique par m f . Dans ce cas, les ket propres communs aux opérateurs 

L 2 et L z se notent simplement par \l,m : 
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Remarquons ici qu'un état décrit par un ket donné renferme 21 + 1 sous-états 

appelés composantes. Dans ce cas, 21 + 1 est dénommée "multiplicité 
orbital", c'est-à-dire le nombre de composantes associées à une orbitale de symétrie I . Par 
exemple, une sous-couche 2p = |l,m^ d'un atome renferme trois sous-couches 
|l,+l), 1 1,0), |l,-l) correspondant aux trois orbitales atomiques 2p x , 2p y et 2p z . 

1-8° Moment cinétique d'un système de plusieurs particules 

Le moment cinétique d'un système à n particules est défini classiquement par un 
vecteur L somme des vecteurs individuels L(v) de chaque particule v : 

L = £l(v) =XL x (v)i+XL y (v)j + ZL z (v)k = L x i+L y j+L z k 

avec : 

L x =ÈLx(v), L y =XL y (v), L z =ÉL z (v) 

Les composantes quantiques s'écrivent : 

L x =ZL x (v), L y =XL y (v), L z =ÎL z (v) 

Les relations de commutations que nous avons établi pour une particule sont valables ici 
pour un système à plusieurs particules : 

[L x ,L y ] = i&L z Q 

En effet : 

[L x , L y ] = [SL x (v),ZL y (p)] =ZZ[L x (v),L y (p)] = £[L x (v),L y (v)] + S £[L x (v),L y (p)] 

vu vu v 

Le second terme où v ^ p est nul, car les opérateurs individuels L x (v) et L y (n) qui portent 
sur des variables v et p différentes commutent. D'où : 

[L x , Ly] = £[L*(v),L y (v)] = i^L z (v) = i»L z 

De même, on peut montrer que L 2 commute avec chacune des composantes du moment 
cinétique : 

[ L 2 , L q ] = 0 q = x, y ou z 

Le carré du moment cinétique se définit comme le produit de l'opérateur L par lui-même : 
L 2 = L.L = EEL(v)L(n) 

V H 

D'où : 

L 2 = ZZ[L x (v)L x (p) + L y (v)L y (p)+L z (v)L z (^)] 

V H 

On peut montrer facilement que l'opérateur L 2 , pour un système à plusieurs particules. 
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peut se mettre sous la forme suivante : 

L 2 = L + L_+L 2 -/zL z 

c'est-à-dire : 


tels que : 


L 2 =ZZL + (v)L-(m) + ZZL z (v)L z(m )-^L z (v) 

vu vu v 

L + =lL + (v), L_=£L_(v), L z =£L z (v) 


Dans le cas d'un système à plusieurs particules, les kets propres communs aux 
opérateurs L 2 et L z se notent plus simplement par : 


où L et Ml sont les nombres quantiques totaux. Dans ce cas, on distingue les équations 
relatives à une seule particule v de celles relatives à plusieurs particules : 


L z (v)|^,m^ = 

L+Cv^.m*) = h^£(M) - m(m+l) 
L _(v)|^,m / ) = h^£(£+\)-m(m-\) 
-l< m £ <+£ 

= 0 , ± 1 , +£ 

£ = \m. t | max 


l 2 |l,m l ) = ^ 2 l(l+i)|l,m l ) 
l z |l,m l ) = /m l |l,m l ) 

L + |L,M l ) = ^L(L+1 )-M l (M l +1) |L,M l + 1) 
L _|L,M l ) = ^L(L+1 )-M l ( m l-1) |L,M l - 1) 
- L < M l <+L 
m l = E m K v ) 

L=ZW(v)-U^0 


Dans la série L = ^f(v),^^(v)-1,...>0, le nombre L varie de pas d'unité à partir de sa 
valeur maximale L(max) = 'Yjl(y) jusqu'au dernier nombre qui doit être nécessairement 
positif ou nul. Cela revient aussi à écrire la série comme suit : 

L(max), L(max)-1, ..., > 0 

Pour L donné, il existe 2L+1 kets | L, M L ), c'est-à-dire | L, L), | L, L-l), ..., |L, - L), formant 
ce qu'on appelle les microétats, c'est-à-dire les composantes possibles de l'état. 

Exemple : 

On considère un système à deux électrons dans une configuration p 2 . La sous-couche p est 
de symétrie 1 = 1 et chaque électron peut occuper l'une des trois orbitales atomiques p x , p y 
et p z représentées par les kets : 

|1, l),|l, 0) et |l, -l) 

Les différents états | L, M L ) du système à deux électrons s'obtiennent en précisant les 
valeurs possibles de L et M L . Deux méthodes existent : 
l?re méthode : 

On part de la définition L = £(i) + 1(2), l(\) + £( 2) - 1, . . . > 0 pour déterminer la valeur 
maximale de L c'est-à-dire L(max) =1+1=2. En variant par pas d'unité on obtient L = 2, 1, 0.: 
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D'où : 

L = 2 -> M l = -2, -1, 0, 1, 2 

L = 1 M l = -1,0,1 

L = 0 ->• M l = 0 

Ainsi, la configuration p 2 engendre 9 microétats qui sont désignés par : 

1 2, -2), 1 2, -1), 1 2, 0), 1 2, +l), 1 2, +2) , 

| 1,4 | 1 , 0 ), | 1 , 4 ), 

| 0 , 0 > 

2ème méthode : 

La deuxième méthode consiste à composer les nombres quantiques individuels selon la 
relation M L =Xm f (v), elle est dite méthode de composition des moments. Dans le cas de 
la configuration p 2 , on compose m e (1) avec rry (2) : 


\ m ^( l ) 

1 

0 

-1 

1 

2 

1 

0 

0 

1 

0 

-1 

-1 

0 

-1 

-2 


Nous identifions ainsi trois séries de nombres : 


M l = -2, -1, 0, 1, 2 

-> 

L = 2 

M L = -1,0,1 

-> 

L = 1 

m l = 0 

-> 

L = 0 


qui correspondent bien aux 9 microétats trouvés par la première méthode. 



(q e > “e) 


1-9° Moment magnétique orbital d'un électron 

Classiquement, une charge q e en mouvement circulaire 
crée un moment magnétique p proportionnel au moment 
cinétique L : 

f 9e 

p = y l y=ir~ 

2m e ^ 

y étant appelé rapport gyromagnétique, m e est la masse de la particule. Dans le cas 
d'un électron (q e =-e), le moment magnétique orbital g crée par le mouvement de cet 
électron autour du noyau est : 

e - e 

g = - L ÙZ=-^ L 7 


2m„ 


2m p 


L'opérateur p z associé est : 


p z = L 

2m e z 

Comme les deux opérateurs commutent, ils admettent un système commun de fonctions 
propres. Soit : 

Pz'P = — Lz'L =.^Æ =0, ±1, ...±l 


2m e 


2m„ 


- 56 - 


Mécanique Quantique S MC (S4)~ Tès 


Chapitre III 


On voit que les valeurs propres du moment magnétique sont quantifiées par le nombre 
quantique magnétique . Le rapport p B = etil 2m e est appelé le magnéton de Bohr. 

L'énergie d'interaction du moment magnétique g avec un champ magnétique B 
appliqué est : 

W = -pB 


Si ce champ magnétique est appliqué dans la direction Oz, cette énergie devient : 
qui lui correspond l'opérateur : 


eB 7 

w=-u 7 b 7 = — m 7 

z 2m z 


W=i^L 2 

2m, 


et les valeurs propres associées sont : 


w ehB z 

W = -YO. e 

2m„ * 


= 0 , ± 1 , ... ±£ 


On voit que l'énergie d'interaction est dépendante du nombre quantique magnétique . 
Ainsi, l'action du champ magnétique sur un atome provoque une levée de dégénérescence 
des états atomiques suivant le nombre quantique magnétique m^. En d'autre terme, la 
dégénérescence d'ordre 2^ + 1 d'une sous-couche I est levée par l'action du champ 
magnétique : 


- cas des orbitales ns : 

£ = 0, m ^=0 => W = 0 

le champ magnétique n'a pas d'effet sur les états ns 

- cas des orbitales np 


E 



m= +1 
= 0 


, i , _ ohB 7 _ eM3 7 

£ = 1, rri/ = -1, 0, +1 => W = —, 0, H — 

2m e 2m e 

le champ magnétique provoque une levée de 


dégénérescence des orbitales np 

B = 0 

0 


Exemples : 

- U : ls 2 2si 

m, = 0 0 +10-1 

Ml = 2(0+0) 


Ml = 0 

- B : ls 2 2s 2 2p T 

mm pi 1 1 

Ml = 2(0+0) +1 

— > 

Ml = 1 

- C : ls 2 2s 2 2p 2 

QDffl 1 1 |t | | 

Ml = 2(0+0) + 1 + 0 


Ml = 1 

- Ne : ls 2 2s 2 2p 6 

lui nul itiluini 

Ml = 2(0+0) + 2(l+0-l) 

-+ 

M l = 0 


Dans cette série d'atomes, seuls B et C possèdent un moment magnétique non nul, donc 
interagissent bien avec un champ magnétique. Notons que les sous couches complètement 
remplies et les sous-couche s ont toujours Ml = 0 et dans le comptage on ne prend en 
compte que les sous-couches incomplètes. Par exemple, pour l'atome de scandium qui 
possède la configuration [Ar]4s 2 3d 1 , le cœur Ar et 4s ont Ml = 0 et seul 3d donne la valeur 
effective du nombre quantique magnétique de Sc. 
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II- Moment cinétique de spin 


11-1° Notion de spin 

L'expérience de Stern et Gerlach (1922) a permis de 
mettre en évidence les propriétés de spin des systèmes 
électroniques. Ainsi, un électron possède non seulement un 
moment cinétique orbital L dû à son mouvement orbital 
autour du noyau, mais aussi un moment cinétique de spin S lié à son mouvement interne. 
Le dispositif expérimental de cette expérience est illustré par le schéma suivant : 


source 

d'argent 


magnétique , P"**"* . 

photographique 

Une source produit des atomes d'argent qui sont collectés par un collimateur afin d'avoir 
un faisceau atomique fin pouvant frapper une plaque photographique situé à une distance 
du collimateur. Lorsque aucun champ magnétique n'est appliqué, on observe une seule 
tache au point O, alors que si ce champ est appliqué dans une direction perpendiculaire au 
trajet du faisceau on observe deux taches aux points A et B symétriques par rapport à O 
(OA = OB). Les deux taches observées indiquent que le moment magnétique (ou moment 
cinétique) de l'atome a interagit avec le champ magnétique. Or comme l'atome d'argent 
dans sa configuration [Kr]4d 10 5s 1 a Ml = 0, ce champ a dû interagir avec un moment autre 
que le moment magnétique orbital. Goudsmit et Uhlenbeck (1925) ont dénommé ce 
nouveau moment par le moment magnétique de spin jî s lié au moment cinétique de spin S 
par la relation : 



où g est le facteur de Landé prenant la valeur 2 dans le cas d'un électron. 

Les deux taches observées sont la preuve évidente de l'existence de deux états 
quantiques de spin et doivent être alors caractérisés par un nouveau nombre quantique de 
spin noté s. Comme la multiplicité d'un état (nombre de composantes possibles) est donnée 
par 2s+l, on aura : 

2s + 1 = 2 s = i 



champ 



Ainsi, cette expérience montre que le spin d'un électron est 1/2. Les deux taches observées 
en A et B doivent correspondre à deux nombres quantiques magnétiques +m s et -m s 
respectivement tel que : 


-s< m s <+s 
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La manifestation de la propriété de spin n'est pas seulement propre à l'électron, 
d'autres particules manifestent aussi cette propriété. Il existe des particules douées de spin 
entier s=l comme le photon, ou bien de spin nulle s=0 comme le méson. Ainsi, en 
mécanique quantique, les particules élémentaires doivent être distinguées entre elles non 
seulement par la masse m et la charge q, mais également par le spin s : 

(m, q, s) 

Le moment cinétique de spin S n'a pas d'équivalent classique comme le moment cinétique 
orbital L; il doit être considéré comme une quantité purement quantique. Son appareil 
mathématique est fondé par analogie formelle au moment cinétique orbital où il doit 
posséder les mêmes propriétés que L : 

- être linéaire et hermitique, 

obéir aux mêmes lois de commutations, 
obéir aux mêmes équations aux valeurs propres. 

Ainsi, on admet l'existence de l'opérateur de spin S défini par ses trois composantes 
(S x ,S y ,S z ) tel que S 2 =S X +S y +S 2 . On définit alors les relations de commutation 
suivantes : 

[S x ,S y ] = i/ïS z [S y ,S z ] = i/zS x [S Z ,S X ] = iftS y 

[S 2 ,S x ] = 0 [S 2 ,S y ]=0 [S 2 ,S z ]=0 

On définit également les opérateurs escalateurs S + = S x + iS y et S_= S x - iS y qui doivent 
vérifier les relations de commutation suivantes, 

[S Z ,S ± ] = ±ihS+ [S + ,S_] = 2hS z [S 2 ,S ± ]=0 

et obéir aux relations produits : 

s + s_=s 2 -s z (s z -h) 
s_s + =s 2 -s z (s z +h) 

La propriété de commutativité des opérateurs S 2 et S z leur confère la propriété d'avoir un 
système commun de fonctions propres. Pour une particule de spin s donné, une telle 
fonction commune X doit vérifiée les deux équations aux valeurs propres : 

S 2 X = /? 2 s(s + l)x 
S Z X = fin i s X 

tel que : -s < m s < +s 

La fonction de spin X est une fonction particulière impliquant une variable discrète notée a 
ne pouvant prendre que certaines valeurs. De telle fonction X(a) décrit l'état quantique de 
spin d'une particule, elle dépend des valeurs des deux nombres quantique s et m s et qu'on 
la note par X s ,m S . Dans l'espace d'états direct, ces fonctions sont représentées par les vecteurs 
kets |s,m s ) et leurs conjuguées dans l'espace dual par les bra (s,m s | associés : 

|s,m s ^ e ë — > (s,m s | e S* 
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Dans le cas d'un seul électron, cet espace est de dimension deux formé par les ket + -0 
et | ^ , - -0 qu'on note communément par les fonctions a et (3 respectivement : 

a = 1 2 _ 5 2 ') symbolisé par | f | 

(3 = J 2 » " 2 } symbolisé par | | | 

Etant donné que les opérateurs S 2 et S z sont supposés hermitiques, les deux fonctions de 
spin a et (3 forment une base orthonormée complète à deux dimensions: 

( a l“> = (p|p) = 1 

(a|p) = (P|a) = 0 

Il s'ensuit que S 2 et S z appliqués à ces fonctions donnent : 

S 2 a = -|ft 2 a S 2 p = |^ 2 p S z a = +^fia S Z P = -^P 

Les applications des relations relatives aux opérateurs escalateurs : 

S + 1 s, m s ) = h^js ( s + 1) - m s (m s + 1) | s, m s + 1) 

S_|s,m s ) = ^^/s(s + l)-m s (m s -l) |s,m s -l) 

permettent de déduire les résultats de l'action de ceux-ci sur les fonctions de spin a et P : 

S + a = 0 S + p = /ux S_a = h$ S_p = 0 

11-2° Relation avec l'hamiltonien de Schrôdinger 

Dans l'approximation non relativiste, l'hamiltonien d'une particule soumise à un 
potentiel de type central : 

fi 2 

H = A + V(r) 

2m e 

est un hamiltonien indépendant de la variable de spin a. Par suite, les opérateurs H , S 2 et 
S z commutent deux à deux : 

[H,S 2 ]=0 [H,S z ]=0 [S 2 ,S z ]=0 

Les trois opérateurs possèdent donc un système commun de fonctions propres. Soit 
T(r,0,cp,a) une telle fonction, on écrit : 

fPF (r , 0, cp, a) = ET (r , 0, cp, a) 

S 2v F(r, 0, cp, a) = h 2 s(s + l)p(r, 0, cp, a) 

S z 'F(r,0,cp, a) = ^m s 'P(r,0,cp,a) 

Comme l'hamiltonien et les opérateurs de spin opèrent sur des variables indépendantes, on 
procède à la séparation des variables en posant : 

¥(^0^,0) = ®(r,0,cp)X(a) 
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Il vient : 

HO(r,0,cp) =EO(r,0,cp) 

S 2 X(q) = /? 2 s(s + 1 )X(g) 

S z X(g) = /?m s X(a) 

Ainsi, la résolution de la première équation se fait indépendamment des deux autres 
équations relatives au spin. 

11-3° Spin d'un système à plusieurs particules 

Pour un système à plusieurs particules, les opérateurs S x , S y et S z sont simplement 
somme des opérateurs individuels : 

S x =ZS x (v) S y =I§ y (v) S z =I§ z (v) 

Le carré du moment cinétique de spin S 2 = S 2 +S 2 +S 2 est donné par : 

s 2 =s + s_+s 2 z -ns z 

avec : 

S + =lS + (v) S_=I§_(v) S z =I§ z (v) 

On peut montrer facilement que les opérateurs S 2 et S z d'un système à plusieurs particules 
commutent : 

[S 2 ,S z ]=0 

Par suite, ils possèdent un système commun de fonctions propres. Soit X(ai, 02 , . . ., a n ) une 
telle fonction, il vient : 

S 2 X(a 1 ,a 2 ,...,a n ) = /i 2 S(S + l)x(G 1 ,G 2 ,...,a n ) 

S z X(a 1 ,cj 2 ,...,cy n ) = ^M s X(a 1 ,o 2 ,...,o n ) 

Les nombres quantiques totaux S et Ms doivent vérifier les conditions suivantes : 

-S < Ms < +S 
M s =Im s (v) 

S = Es(v),Es(v)-l,..„>0 

Un état quantique de spin est caractérisé par les deux nombres quantiques S et Ms, et 
chaque composante ou microétat sera représenté par un ket |S,M s ^. Comme Ms prend 
2S+1 valeurs possibles, l'état quantique du système est caractérisé par la quantité 2S+1 
appelée multiplicité de spin, elle prend les appellations suivantes selon les valeurs du 
nombre quantique S : 


S 

0 — 1 — 

u 2 1 2 

2S+1 

Appellation 

1 2 3 4 

Singulet Doublet Triplet Quadruplet 


-61 - 


Mécanique Quantique S MC (S4)~ Tès 


Chapitre III 


Exemple 1 : 

Soit un système à trois électrons occupant la sous couche p, c'est-à-dire un système dans la 
configuration électronique p 3 . Pour recenser les nombres quantiques S et Ms et les 
microétats |S,M s ^ possibles, on utilise la méthode de composition des moments selon la 
relation Ms= m s (l)+m s (2)+m s (3). A cet effet, on doit construire un tableau à trois entrés 
faisant intervenir les trois nombres quantiques m s (l), m s (2) et m s (3). Pour simplifier, on 
procède d'abord par la composition de m s (l) avec m s (2) et ensuite le résultat obtenu 
Ms(l,2)=m s (l)+m s (2) se compose avec m s (3); cela 
revient à considérer la configuration p 3 comme 
étant le résultat de la composition de deux 
configurations p 2 ®p 3 . Pour la configuration p 2 , 
cette composition donne les résultats dans le 
tableau ci-contre. Ainsi, on obtient quatre nombre 
Ms(l,2) = 1, 0, -1, 0. Ces derniers doivent se composer avec le nombre m s (3) pour donner le 
nombre quantique total Ms(l,2,3) de p 3 : 


10-10 

1 1 1 _i 1 

2 2 2 '2 2 

1 1111 
2 2 ~_2 ~_2 ~ 2 




Nous identifions ainsi trois séries de nombres : 


M s 

_ 3 1 1 3 

2’ 2’ 2 ’ 2 


s = ! 

2S + 1 = 4 

M s 

_ 1 1 
2’ 2 

->■ 

s ' = | 

2S'+Ï = 2 

M s 

_ 1 1 

2’ 2 

-> 

s” = ~ï 

2S"+1 = 2 


La configuration p 3 engendre donc trois états distincts: un quadruplet et deux doublets. Les 
microétats contenus dans chaque état sont : 



Etats (2S+1) 

Microétats |S,M s ) 




Quadruplet 
Doublet (1) 
Doublet (2) 

MM 

MM 

Ib+i)' 

M)- |M>- |f -f) 

2 ,_ l)' 

1_1\ 

2’ 2/' 



Exemple 2 : 

- Li : ls 2 2si 

E B 

M s 

= 4-i)=o -, 

s = o. 

2S + 1 = 1 

- B : ls 2 2s 2 2p 1 

e e iïrr 

M s 


Ç — 1 
2 , 

2S + 1 = 2 

- C : ls 2 2s 2 2p 2 

E E 1 1 1 1 1 

M s 

= 2(ï-ï)+2(t)=1 -> 

S = l, 

2S + 1 = 3 

- Ne : ls 2 2s 2 2p 6 

E E «ittltt 

M s 

= 4-i)+3(l-i)=0 -, 

s = o. 

2S + 1 = 1 


On voit que l'état fondamental de Li et de Ne est un singulet, alors que B est un doublet et 
C un triplet. 

Notons que la première méthode ( exemplel ) donnent tous les états possibles 
fondamental et excités issus d'une configuration électronique donnée, alors que la seconde 
méthode ( exemple 2) ne donne que l'état fondamental de la configuration considérée. 
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III- Moment cinétique total 

Pour une particule quantique, les moments cinétiques 
orbital L et de spin S se composent pour donner une résultante 
J appelé moment cinétique total : 

J = L + S 

Si L et S caractérisent respectivement le mouvement orbital et le mouvement de spin de la 
particule, alors J traduit la manière dont les deux moments L et S sont couplés et orientés 
l'un par rapport à l'autre. 

Les opérateurs composantes J x , J y , J z s'écrivent simplement : 



Ces composantes obéissent aux même lois de commutation que celles établies pour les 
moments cinétiques orbital et de spin : 

[J x ,Jy] = iW» ^ 

[J 2 ,J z ]=0 ^ 

Pour les opérateurs escalateurs J + = J x + i J y et J_ = J x - i J y , on définit aussi : 


J 2 =J + J_+J z (J z -/0 

J 2 = O + +J z (J z +/0 

En calculant le carré du moment cinétique total J 2 , on fait apparaître un terme important 
en physique appelé terme de couplage spin-orbite : 

J 2 = (L + S) 2 = (L + S)(L + S) = L 2 + S 2 + LS + SL 
Comme L et S opèrent sur des variables différentes, ils commutent. Il vient : 

J 2 = L 2 +S 2 +2LS 

D'où l'on tire le terme : 


LS = ^(J 2 -L 2 -S 2 ) 

appelé terme de couplage spin-orbite ou couplage LS. 

La commutativité des opérateurs J 2 et J z entraîne la possibilité de les mesurer dans 
le même état quantique : 

J 2 'F = /* 2 j(j + l>F 
J z v F = /zm j v F 

avec : - j < m j < + j 

m j = xx\ p + m s 

Dans ce cas, la fonction d'onde d'une particule doit dépendre des coordonnées d'espace 
(r,0,cp) et de spin (a) ainsi que du jeu des nombres quantiques j et m j : 

l P = 'P jm (r,e,(p,CT) 

Pour simplifier l'écriture de cette fonction d'onde, on utilise la notation vectoriel | j, m^. 
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Pour un système à n particules, cette fonction d'onde est évidemment fonctions des 
coordonnées des n particules 'P(r 1 ,0 1 ,cp 1 ,a 1 , r 2 ,0 2 ,cp 2 ,a 2 , •••> r n ,0 n ,cp n ,ü n ) et dépendant 
des nombres quantiques totaux J et Mj : 


tels que : 


j 2 t = /? 2 j(j+i) v f 

J z 'F = ^M j v P 
-J < Mj < +J 
Mj = Ml + Ms 


M i=Z m j(v) = Z[ m A v ) +m s( v )] 


On peut montrer que le nombre quantique total J obéit à l'inégalité suivante : 

| L — S | < J <|l + s| 

où J parcourt les valeurs de | L + S | à | L — S | par pas d'unité : 

f L entier positif ou nul 

J = L+S, L+S-l, — , > | L — S | s S entier ou demi-entier, positif ou nul 
L J entier ou demi-entier, positif ou nul 

Pour démontrer l'inégalité |L-S|<J<|L + S|,on précise d'abord les conditions auxquelles 
doivent obéir les nombres quantiques : 


-L < Ml < +L M L (max) = L M L (min) = -L 

-S < Ms < +S Ms(rnax) = S Ms (min) = -S 

-J < Mj < +J Mj(max) = J Mj(min) - -J 

La valeur maximale que peut prendre Mj est : 

Mj(max) = ML(max) + Ms(max) = L + S 

Or, comme Mj(max) = J donc la valeur maximale que peut prendre le nombre J est : 

J(max) = L + S 

De la même façon, on peut démontrer que la valeur minimale que peut prendre le nombre J 
est : 

J(min) = | L S | 

La quantité L-S doit être prise en valeur absolue pour s'assurer que le nombre J soit 
toujours un nombre positif. Le nombre J peut être entier, demi-entier ou nul selon les 
valeurs de L et S. On distingue ici deux cas : 

- Si L > S, on a L - S < J < L + S et le nombre variable dans ce cas est S, par suite J 
prend (2S+1) valeurs possibles. 

- Si L < S, on a S - L < J < L + S et le nombre variable dans ce cas est L, par suite J 
prend (2L+1) valeurs possibles. 
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Exemples : 

- Configuration électronique p 1 : 

Le seul électron a les nombres quantiques 1=1 et s=l. Par suite, le nombre j vérifiant la 
condition | t -s | < j < t +s est : 

|i-i|<j<|ui| - i -H 

En tenant compte du couplage spin-orbite, la configuration p 1 engendre deux états distincts 
l'un a j = ^ et l'autre a j = -| . Ici on a I > s, donc j prend (2s+l)=2 valeurs possibles. 


- Configuration électronique p 2 : 

Pour chacun des électrons, onaf(l)= I (2) = 1 et s(l) = s(2) = j-. Avant de pouvoir recenser 
les valeurs possibles de J, on doit d'abord déterminer les valeurs des nombres L et S selon 
les relations : 

L = I^(v),^(v)-l,...,>0 
S = Es(v),Es(v)-l,..„>0 


on obtient : 


L = 2, 1, 0 et S = 1, 0 


Il existe six manières pour composer les nombres L et S afin d'engendre les nombres J 
possibles : 

(L=2, S=l), (L=2, S=0), (L=l, S=l), 

(L=l, S=0), (L=0, S=l), (L=0, S=0) 


Par exemple, on a : 


(L=2, S=0) ^ 

2 < J 

(L=l, S=l) ^ 

0<J 

(L=0, S=0) -> 

0<J 


<2 


j = 2 

<2 


j = 0, 1, 2 

<0 


j = o 
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I- Introduction 

Un système quantique simple est un système dans le quel une particule est soumise 
à un champ de potentiel bien défini permettant la résolution exacte de son équation de 
Schrôdinger : 

A + V(x, y, z) j W = E'P 

Toute la difficulté pour la résolution de l'équation de Schrôdinger provient du terme 
potentiel V. A chaque champ de potentiel correspond un modèle en mécanique quantique : 

- au potentiel nul V = 0 correspond le modèle de la particule libre ou rotateur libre, 

- au potentiel constant V = Vo correspond le modèle du puit de potentiel ou barrière 
de potentiel, 

- au potentiel central V(x,y,z) = V(r) correspond le modèle de l'oscillateur harmonique 
ou atome d'hydrogène. 

En toute généralité, une particule est apte à effectuer des mouvements dans l'espace de type 
translation, rotation, vibration ou gravitation. 

II- Particule libre 

Une particule est dite libre si elle n'est soumise à aucun champ de force, son énergie 
potentielle est nulle (V=0) et seule l'énergie cinétique caractérise son mouvement : 

h 2 

AT = E V P 

2m 

Dans ce cas, deux types de mouvements peuvent être réalisés par cette particule : 
mouvement de translation ou mouvement de rotation. 


11-1° Mouvement de translation : 


Soit une particule de masse m libre d'effectuer un 

mouvement de translation suivant l'axe x. Un tel mouvement P ^ x 

est caractérisé par la quantité de mouvement (p = mv). m 

L'hamiltonien d'un tel système est : 

~ Pl û h 2 d 2 

2m 2m dx 2 

Comme H et p x commutent, alors une fonction propre de l'un est également fonction 


propre de l'autre : 


P x v E = a p v E « = 

2m 

a p et E sont les valeurs propres associées qui sont nécessairement réelles puisque les 
opérateurs qui les donnent sont hermitiques. Comme p x T = a 2x V , la valeur propre de H 
a p 

est E = — — > 0 . Ainsi, pour une particule effectuant un mouvement de translation libre, 
2m 
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les énergies propres sont positives. Comme E n'est conditionnée par aucun nombre discret, 
le spectre des valeurs propres est continu; l'énergie de translation est donc non quantifiée. 

Les fonctions propres communes aux deux opérateurs H et p x peuvent s'obtenir par 
simple intégration de l'équation différentielle relative à p x : 

T-“j— = ap'E -> — = -> ¥(x) = Aexpfex] 

i dx dx h V n J 

Une telle fonction est particulièrement non finie lorsque x— >+œ et non carré sommable, 
donc elle n'appartient pas à l'espace d'Hilbert décrivant les états liés. En conclusion, le 
mouvement libre de translation est un mouvement non quantifié, il est caractérisé par un 
spectre continu, énergies positives (E > 0) et états non liés ('E non carré sommable). Parmi 
les exemples courants de particules libres, on peut citer l'électron envoyé à l'infini lors d'un 
processus d'ionisation d'un atome : 

A -» A + + e“ 


11-2° Mouvement de rotation : 

a] Rotateur rigide plan : 

Considérons une particule de masse m astreint à se 
mouvoir sur un cercle de rayon r fixe autour de l'axe z avec 
une vitesse angulaire co uniforme : 



Un tel système est dit rotateur rigide plan, la position de la particule est repérée à chaque 
instant par l'angle cp qui parcourt les valeurs de 0 à 2k dans le plan xOy. Comme la 
particule est libre (V=0), son énergie de rotation est purement cinétique : 

E = — mv 2 =— m(cor) 2 = — (mr 2 )co 2 
2 2 2 


En identifiant l'expression du moment d'inertie I = mr 2 , cette énergie devient : 


E = — Ico 2 
2 


Un tel mouvement doit être caractérisé par le moment cinétique ( L = r a p ) suivant l'axe z, 
c'est-à-dire : 

L z = r • mv = mr 2 oo = Ico 

L'expression de l'énergie cinétique du rotateur en mécanique classique est : 



qui lui correspond en mécanique quantique l'opérateur suivant : 

h - a _ 

21 21 dcp 2 
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L'équation de Schrôdinger du système est : 


d 2v F 

dtp 2 


En posant k 2 = 2TE!Ti 2 l'équation devient : 


d 2 T 

dtp 2 


k 2 T = 0 


La solution générale est de la forme : 

T = Ae ikip 


On peut vérifier que cette fonction est une bonne fonction d'onde dans l'intervalle [0, 2n], 
elle est particulièrement de carré sommable ce qui permet de déterminer le facteur de 
normalisation A : 

ff^Wcp = L 27 I A 2 dcp = A 2 271 = 1 A = -L 

J ° J ° V2jt 

mais T n'est uniforme que si ¥( 9 ) = 4 / (cp+27i), c'est-à-dire : 

e lk271 _ cos (] c 27r) +isin(k27i) = 1 

Cette relation n'est vérifiée que si m est un nombre réelle positif, négatif ou nul : 
k = 0, ± 1, ± 2, . . . avec k 2 =2IE /^ 2 


Finalement, les solutions pour le problème d'une particule tournant sur un cercle de rayon 
fixe (rotateur rigide plan) sont : 


E k = 



*k(<p) = 


Æ 




On voit que la quantification de l'énergie est introduite grâce à la condition aux limites 
(finitude de la fonction d'onde dans l'intervalle 0 et 2n). 

Quelque soit le nombre quantique k, c'est-à-dire quelque soit l'état quantique 'Fk, la 
densité de probabilité prend une valeur constante : 


La séparation énergétique entre deux niveaux successif k et k+1 est donnée par : 

AE = — (2k + l) 

21 V ’ 


L'intervalle AE croit avec k. Si le système tournant est un objet macroscopique (masse m 
grande, moment d'inertie I grand), l'intervalle énergétique tend vers zéro (AE— >0) et la 
quantification disparaît; on obtient alors le résultat de la mécanique classique (énergie 
continue). 
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b) Rotateur rigide sphérique : 

Un rotateur sphérique est un système pouvant effectuer des rotations dans toutes les 
directions d'espace. Les molécules diatomiques sont des exemples de rotateurs sphériques. 
Considérons une molécule diatomique AB où les deux atomes de masses iua et mB sont 
séparés par une distance r supposée fixe (rotateur rigide sphérique) : 



Le système pivote autour du centre de masse G et si l'on prend ce centre comme origine on 
peut montrer que le système à deux corps se ramène en un seul corps de masse réduite g 
tournant autour d'un axe £, perpendiculaire à l'axe intemucléaire. La position du centre de 
masse est définie à partir des deux égalités suivantes : 


m A r A = m B r B / 


r = r A+r B 

111 A 


r B =- 


Dans l'hypothèse d'un rotateur où les deux masses ne sont soumises à aucun champ de 
force, l'énergie de rotation est une énergie purement cinétique : 

E = Lm A vi +^-m B v| 


Si la rotation est uniforme, c'est-à-dire que les deux masses tournent avec la même vitesse 
angulaire co , on écrit : 

V a V R 

CO = — = — 

*A r B 

L'énergie cinétique devient : 

E = |m A r|co 2 +^m B r|co 2 = ^-m A rl + m B rg co 2 = ^Ico 2 
où I dénote le moment d'inertie de la molécule : 


I - m A r A + m B r B 


m A m B j.2 
m A +m B 


= ùr 


2 


Nous trouvons donc que le moment d'inertie des deux masses mA et mB se ramène à celui 
d'une simple masse g dite masse réduite. En vertu de ces hypothèses, le seul moment 
d'inertie non nul est le moment par rapport à l'axe E,, appelé axe principal d'inertie. Ainsi, 
la seule rotation que l'on doit considérer dans le cas de molécules diatomiques est celle qui 
s'effectue autour de cet axe. 
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Comme le moment cinétique est une caractéristique du mouvement de rotation, on 
doit exprimer l'énergie E en fonction de celui-ci. Le moment cinétique d'un système à deux 
masses est donné par : 

L = m A (r A AV A ) + m B (r B av b ) 

Comme ij _L Vj, le module de L se réduit à : 

L = m A r A v A +m B r B v B = (m A r A +m B r B )co = Icq 
L' énergie de rotation classique devient : 



Ce résultat montre que le rotateur est en continuel pivotement autour du centre de masse et 
son énergie E est donc continue. 

En mécanique quantique la situation du rotateur est différente, on montrera que son 
énergie est discontinue. L'opérateur hamiltonien H associé à l'énergie classique E est : 



Dans ce cas, les variables qui repèrent la position de la masse 
réduite g dans l'espace sont les angles (0,cp), la distance r à 
l'origine étant supposée fixe. On doit donc utiliser les 
coordonnées sphériques pour exprimer l'hamiltonien : 

x = rsinG coscp 
y = rsinG sincp 
z = r cosG 

En reportant l'expression du moment cinétique en coordonnées sphériques dans H, on 
obtient : 


Comme les deux opérateurs H et L 2 commutent, ils admettent Y(0,cp) comme fonctions 
propres communes : 

HY(0,cp) = E Y(0,cp) 

L 2 Y(0,cp) = ft 2 J(J+l)Y(0,cp) 

D'où, les valeurs propres du rotateur sphérique rigide : 
h 2 

E t = — — J(J +1) J = 0, 1,2, ... 

J 21 

J étant le nombre quantique de rotation moléculaire. Ainsi, on montre que l'énergie du 
rotateur diatomique est quantifiée par le nombre discret J. Les fonctions propres Y (0,cp) ne 
sont autres que les harmoniques sphériques données dans le chapitre III. 

III- Puits de potentiel - particule liée 

En principe, une particule liée est une particule constamment soumise à un champ 
de potentiel qui la piège. Une telle particule est dite confinée ou piégée dans le puits de 
potentiel. En physique, plusieurs formes de puits existent pour modéliser les systèmes 


û » i ô . . a î 

H = — sinG — -\ - 

21 |^sin0 50 50 sin 2 0 5cp 
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quantiques, comme le puits constant, le puits parabolique ou le puits coulombien. 


V(x) = V 0 V(x) =1/2 kx 2 V(r) = - e z /r 



autres cas seront traités dans les paragraphes qui suivent. 

111-1° Solutions du problème 

Les deux caractéristiques importantes d'un puits de potentiel constant sont sa 
largeur (distance entre parois) et sa profondeur (hauteur des parois). Considérons une 
particule de masse m piégée dans un puits de largeur L et de profondeur infiniment grand : 


Les parois du puits supposées infinies veut dire que la particule ne peut pas y échapper, 
une fois piégée elle y demeure indéfiniment. A l'intérieur du puits, le potentiel est constant 
qu'on peut poser nul, alors qu'à l'extérieur le potentiel est infini : 

V(x) = 0 0 < x < L 

V(x) = oo ailleurs 

Ainsi, la probabilité de trouver la particule à l'extérieur du puits est nulle, c'est-à-dire que la 
fonction d'onde associée doit être nulle dans ces régions : 

x < 0 ¥(x) = 0 

x > L ¥(x) = 0 

et également nulle lorsque que sa position x atteint les parois du puits : 


¥(0) = 0 ¥(L) = 0 


Ces deux relations constituent les conditions aux bords que la fonction d'onde doit 
satisfaire. Par conséquent, on ne doit résoudre l'équation de Schrôdinger qu'à l'intérieur du 
puits : 


d 2v P 

dx 2 



avec 0 < x < L 


En posant k 2 = 2mE//z 2 , l'équation devient : 


d 2 ^ 

dx 2 


k 2x ¥ = 0 
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On obtient donc une équation différentielle linéaire à coefficients constants, sa résolution 
peut se faire au moyen de la méthode de l'équation caractéristique (ou équation auxiliaire). 
Dans ce cas, on cherche les solutions de la forme : 

sx d 2v P 2 

'F(x) = e => — — = s 2v P(x) 

dx 2 

s étant en général un nombre complexe. Après substitution dans l'équation différentielle, on 
obtient : 

s 2 + k 2 = 0 -> s 2 = -k 2 = (ik) 2 | 1 

[s 2 = -ik 

Deux solutions particulières sont alors possibles : 

v P 1 (x) = e +ik et v F 2 (x)=e- ik 

La solution générale est une combinaison linéaire des deux solutions particulières, c'est-à- 
dire : 

'F(x) = ae +ik +be -lk = A cos kx + B sin kx 

avec : 

A = a + b et B = i(a - b) 

Les coefficients A et B sont des constantes arbitraires mais qui peuvent se préciser en 
appliquant les conditions aux bords et la condition de normalisation : 

- la fonction s'annule en x = 0, 

'Ffx^O => A = 0 => 'L(x) = Bsinkx 

- la fonction s'annule en x = L, 

¥(01 = 0 => B sin kL = 0 => sinkL = 0 

La relation sin kL = 0 n'est vraie que si le facteur kL est un multiple de k : 
kL = ±njr et n = 1, 2, 3, . . . 

On trouve ainsi que le nombre n doit être un entier positif mais jamais nul; car si n=0 on 

aurait V L=0. De cette relation, on tire les énergies propres et les fonctions propres de la 

particule dans un puit infini : 

E r-n 2 , 'Ljj (x) = B sinf — xl 

8mL 2 IL J 

Le facteur de normalisation B a pour valeur : 

{'r|'P) = B 2 J„ L sm 2 (=]dx=l => B = Jï 

D'où : 

' p “ (x)= Jr sin (f x ) 
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On peut montrer facilement que les fonctions propres du problème sont orthogonales deux 
à deux : 


e?» h»’) =flo Lsin ( n r) sin ( i r^) dx = s " 


111-2° Quantification de l'énergie 

L'énergie d'une particule dans un puits de potentiel est quantifiée. 


E n 


h 2 

L n 

8mL 2 


2 


elle ne peut prendre que certaines valeurs conditionnées par le nombre entier n. Cette 
quantification est une conséquence directe des conditions aux bords (x=0 et x=L). L'écart 
énergétique entre deux niveaux successifs n et n+1. 


AE = 


h 2 

8mL 


(2n + l) 


croit monotonement avec le nombre quantique n. Deux cas peuvent être envisagés : 

- si la masse m croit (objet macroscopique), AE— >0, la quantification disparaît, la 
variation de l'énergie est continue et on obtient le résultat de la mécanique classique. 

- si la largeur L du puits croit, l'espacement entre les niveaux d'énergie décroît, 
lorsque L devient très grand, AE->0, et la quantification disparaît. 


111-3° Densité de probabilité 

La densité de probabilité d'une 
particule piégée dans un puits de potentiel 
infini est : 



Cette densité est une fonction sinusoïdale 
montrant que la particule effectue des 
mouvements d'oscillation entre les deux 
parois du puits. La fréquence d'oscillation 
peut être déterminée à partir de la relation 
de Planck : 

v = AE/h = — (2n + l) 

8mL 

Aussi, le nombre de nœuds dépend du 
rang énergétique de l'état quantique : 



Un nœud est un lieu où la densité de probabilité est nulle 
| | 2 = 0 
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D'après le diagramme, on voit que le nombre de nœuds augmente comme (n-1) lorsque n 
augmente. Le niveau fondamental (n=l) n'a pas de nœuds, le premier niveau excité (n=2) a 
1 nœud, le second a 2 nœuds, . . .etc. Les abscisses x n donnant ces nœuds sont : 



Par exemple, pour le niveau fondamental (n=l), on a xi = kL (k=0,l,..), et à l'intérieur de 
l'intervalle [0, L] il n'existe aucun noeuds. Par contre pour le premier excité (n=2), on 
a x 2 = jL, ce qui donne un nœud à l'intérieur de l'intervalle. 


IV- Oscillateur harmonique - puits parabolique 


IV-1° Energie potentielle 

Sous sa forme à une dimension, l'oscillateur harmonique est un système constitué 
d'une particule de masse m plongée dans un potentiel parabolique dit aussi potentiel 
harmonique : 

V(x) = — kx 2 
2 


k étant la constante de force de rappel (F=-kx). Le point x = 0 correspond à la position 
d'équilibre de la masse oscillante, en ce point l'énergie potentielle est minimum et la 
courbure est positive (parabole positive) : 



= 0 


= k > 0 


Ces deux quantités correspondent en fait aux deux premiers termes du développement en 
série de Taylor de V(x) autour de x=0 : 


V(x) = V Q + 


xUfeï) 


ôxj 0 2 l[ôx 2 J 0 


Vo=V(x=0) étant l'énergie du minimum qu'on peut choisir comme origine. Une telle énergie 
potentielle V(x) = V 2 kx 2 est représentée schématiquement sur la figure (a) : 


V(x) 



V(x) 



Figure (b) 


V(x) 
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Si le l'origine du potentiel est quelconque, on peut toujours développer V(x) en série de 
Taylor autour de la position d'équilibre x 0 # 0 (figure b) : 


av 


v « = v o + bd(*- x o)+d 


1 ( a 2 v 


2! i dx 2 


(x-x 0 ) 2 +•• 


Pour x = x Q , l'énergie est minimale, le gradient est nul et on écrit : 

V(x) = V Q + ~k(x - x Q ) 2 H — 

Si l'on choisit Vo comme origine (Figure c), on obtient finalement l'énergie potentielle d'un 
oscillateur harmonique centré en point x Q : 

V(x) = |k(x-x 0 ) 2 

IV-2° Traitement classique : 

En vertu de la relation fondamentale F=my de la mécanique classique, on peut 
déduire l'équation de mouvement : 


d 2 x , d 2 x k __ d 2 x 2 

- = -kx — —h — x = 0 — — + MX = ü 

dt 2 m 


dt 2 dt 2 m dt 2 

Les solutions de telle équation sont de la forme : 

x(t) = Asin(cot + cp) 

avec (o = Vk/m : pulsation propre de l'oscillation 
cp : phase de l'onde 

A : amplitude 

La fréquence de l'oscillation v = co/ 2 tt est une caractéristique du mouvement vibratoire : 

vmJ-M 

2n V m 

Pour un oscillateur donné (co ou v donné), on montre que l'énergie totale ne dépend que de 
l'amplitude A : 

„ 1 fdxŸ 1 , 2 1 2 . 2 

2 Ut J 2 2 

Comme l'énergie est constante, la mécanique classique prévoit que le spectre d'énergie de 
l'oscillateur harmonique est continu. 

IV-3° Traitement quantique : 

En partant de l'expression de l'énergie potentielle centrée au point x=0, l'équation de 
Schrôdinger s'écrit : 


d 2 ¥ 2m f 1 

dx 2+ ^ 2 l 2 

Sachant que co = V k/m , l'équation devient : 


d 2v F f 2mE m 2 co 2 


dx 2 


h 2 h 2 
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On peut facilement montrer que le facteur h ! mco a la dimension de longueur qu'on notera 


h 


x 0 est appelé communément longueur fondamentale de l'oscillateur. L'introduction de la 
nouvelle variable E, et du paramètre s , tous les deux sans dimensions. 


x Q fit o 

conduisent alors à la forme réduite de l'équation de Schrôdinger : 


d 2x F 

dq 2 


+ (s-^ 2 ) v P = 0 


Cette équation n'admet des solutions physiquement acceptables (finies et carré sommables) 
que si le paramètre s est un nombre entier impair : 

s = 2v+l 

De cette condition, on tire les valeurs des énergies possibles : 

E y = (v + l)/ko <^> E y = (v + ^)hv v=0,l,2, ... 

Elle ressort donc que l'énergie de l'oscillateur harmonique est quantifiée par le nombre 
discret v appelé nombre quantique de vibration. Ainsi, contrairement à la mécanique 
classique, la mécanique quantique montre que l'énergie de l'oscillateur est discrète, elle ne 
peut prendre que certaines valeurs permises. Dans ces conditions, la fonction d'onde 
vérifiant l'équation différentielle : 

d 2v P , 

-J- + (2v + l-£ 2 ) x F = 0 

d^ 2 v 

est une fonction formée de produit d'un polynôme H v (£) appelé polynôme d'Hermite et 
d'une exponentielle e -q /2 : 

t v (q = n v h v (q e -^ 2/2 

Sans rentrer dans les détails de calcul, on peut montrer que le facteur de normalisation N v 
et le polynôme d'Hermite H y (^) prennent les expressions suivantes : 

N v =(2 v v!x 0 Vrc) -1/2 

H (4) = (-l) v e^ 2 ~~( e_ ^ 2 ) 
d£, 

Les expressions des premiers polynômes d'Hermite s'écrivent : 

H 0 © = 1 

Hj© = 2^ 

H 2 (Ç) = 4Ç 2 -2 
H 3 (^) = 8E, 3 -12E, 
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On obtient deux familles de polynômes : polynômes pairs (si v pair) et polynômes impairs 
(si v impair). Les premières fonctions d'onde de l'oscillateur harmonique et leurs énergies 
s'écrivent : 


Les fonctions propres de l'oscillateur harmoniques sont réelles et forment un système 
orthonormé complet : 


En vertu de l'expression générale des fonctions d'onde de l'oscillateur harmonique : 


la parité de celles-ci est (-l) v . Ainsi, HL est pair si v est pair et inversement. La fonction HL 
n'a pas de racine E, donc paire, alors que HL a une racine donc impaire. On retrouve un 
résultat générale pour tous les états liés d'une particule à une dimension : les fonctions 
d'onde sont alternativement paires et impaires (l'état fondamental est toujours pair) et le 
nombre de racines (nombre de nœuds) correspond au degré d'excitation v. Cette propriété 
de parité est très importante surtout pour l'évaluation des intégrales telles que : 


L'intégrale est non nulle si le produit H / v x I1 H / v . est pair. Comme H^HV a la parité (-l) v+v et 
x n a la parité (-l) n , alors la parité totale est (-l) n+v+v . L'intégrale (H^ |x n |T v .) est non 
nulle si n + v + v' est pair, c'est-à-dire si n et v+v' sont de même parité. 

IV-4° Spectre d'énergie 

Les énergies propres de l'oscillateur harmonique sont dépendantes du nombre 
quantique de vibration v, elles peuvent s'exprimer de trois façons différentes : 

V(x) 


E o = 2 hv -> T 0 =(x 0 V^r 1 / 2 e -^ 2/2 

E i = i hv -> %=(^r 1,2 ^ /2 

E 2 =fhv -> T 2 = (2x 0 V^r 1/2 (2^-l)e^ 2/2 
E 3 -Jhv -> H> 3 =(3x 0 V^)- 1/2 (2^ 3 -3^)e-^ /2 


x _ mco 
x 0 n 


(W y |HV) = J^H / V H / V .dx = ô vV 



(T v | x n |T v .) = J^H / v x n H / v ,dx 


E 










avec v = 0, 1, 2, 3, ... 


o 
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a» = Vk/m est la pulsation, v = ta Un la fréquence, k la constante de force et m la masse de 
la particule. Les énergies E v sont strictement positives et jamais milles. En effet, même pour 
le niveau fondamental (v=0) l'énergie est non nulle : 

E 0 =ihv 

Eo est appelé "énergie du point zéro" qui montre que la vibration d'une particule ne cesse 
jamais. Aussi, nous pouvons remarquer que les niveaux d'énergie sont équidistants de la 
quantité hv : 

E v +i - E v = hv 


IV-5° Vibration de molécules diatomiques 


Une application évidente du modèle de l'oscillateur harmonique est l'étude du 
spectre de vibration de molécules diatomiques. Pour une molécule AB, l'énergie potentielle 
V(R) en fonction de la distance intemucléaire R 
peut être modélisée par la fonction de Morse : 

V(R) = -V 0 +D 0 [l-e" a(R_Ro) ] 2 

Ro : distance d'équilibre 
Vo : énergie d'équilibre 
Do : énergie de dissociation 
a : paramètre caractéristique de la liaison A-B 
Nous remarquons qu'à la distance d'équilibre Ro l'énergie est minimale V(R 0 ) = -V 0 . 
L'énergie de dissociation s'obtient par simple différence d'énergie entre produits et 
réactants : 

D 0 = V(A) + V(B)-V(AB) 



Si l'énergie V(A) + V(B) des atomes A et B sont choisies comme origine, c'est-à- 
dire V(A) + V(B) = 0, alors Do n'est autre que Vo (profondeur du puits). L'énergie de 
dissociation est une quantité nécessairement positive car elle correspond à un processus 
exigeant de l'énergie. 


Fa F b 

WWW - 

M a G M b 


Pour traiter la vibration moléculaire en 
mécanique classique, on suppose que la liaison 
chimique est assimilée à un ressort de constante de 
raideur ou constante de force k. Dans ce cas, les deux R A Rb 

atomes peuvent osciller (ou vibrer) autour de leur R 

position d'équilibre le long de l'axe internucléaire (liaison est non rigide). La force de rappel 
est proportionnelle à l'élongation. Cette élongation est la différence R - Ro entre la distance 
intemucléaire d'équilibre Ro et la distance instantanée R. En appliquant la loi fondamentale 
de la mécanique classique à chacun des atomes, on obtient : 



d2R A 


m a 

dt 2 

= -k(R-R 0 ; 


d 2R B 


M b 

“dt 2 " 

= -k(R-R 0 ) 


où Ra et Rb sont les distances instantanées des atomes A et B par rapport au centre de 
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masse G de la molécule. Or, d'après la définition du centre masse (voir II- rotateur rigide), on 
a : 


R a = 2 R 

A M a +M r 


Rr = — R 

B M a +M r 


et les deux équations de mouvement se ramènent en seule équation : 
M a M b d 2 R 


M a +M r dt 2 


= -k(R -R 0 ) 


^ = -k(R-R 0 ) 

p étant la masse réduite. On définit R -R 0 = x comme étant la variable de déplacement ou 
d'oscillation. D'où : 


d 2 x 

dt^ 





d 2 x 

dt^ 


co 2 x = 0 


On obtient une équation identique à celle rencontrée pour une masse m oscillante (voir plus 
haut) dont ses solutions sont de type : 

x(t) = Asin(cot + cp) 


En mécanique quantique, le traitement de la vibration d'une molécule diatomique est 
analogue à celui de l'oscillateur harmonique : 


d 2v P 

dx 2 



ikx= 

2 


W = 0 


à condition de remplacer la masse m par la masse réduite p et de considérer x comme étant 
la variation R - R 0 . 


V- Atome d'hydrogène - puits coulombien 

V-l° Equation de Schrôdinger et ses solutions 

L'atome d'hydrogène est un système coulombien où l'électron est piégé dans un 
champ de potentiel de type central. Le système est composé de deux particules : un électron 
de masse m e (charge -e) et un proton de masse m p (charge +e). Leur interaction 
électrostatique est donnée par le potentiel d'attraction : 

V(r) = 

r 

r étant la distance qui sépare l'électron au noyau. Ce problème à deux corps (électron, 
proton) peut se réduire en un simple corps de masse réduite p si l'on choisit l'origine au 
centre de masse G. 



x 
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Etant donné que la masse de l'électron est beaucoup plus petite que celle du proton 
(m p =1836m e ), on peut négliger m e devant m p : 


Cette approximation est toujours valable quelque soit le système moléculaire, ce qui permet 
de prendre l'origine O du trièdre xyz au noyau plutôt qu'au centre de masse G. Dans telle 
approximation, l'hamiltonien de Schrôdinger s'écrit : 

2rn 


H = 


-A 


Comme le potentiel V(r) est de symétrie sphérique, l'utilisation des coordonnées sphériques 
s'impose : 


~ h 2 f 1 ô , ô 1 ô . A ô 1 ô 2 ^ e 2 

H = — — r 2 — + — sin 0 — + — 

2m e ^r 2 ôr ôr r 2 sin0 50 50 r 2 sin 2 0 5cp 2 J r 


Cet hamiltonien peut encore s'exprimer en fonction de l'opérateur carré du moment 
cinétique L 2 : 

H - _ A + L 2 _ e^_ 

2m g r 2 5r ôr 2m e r 2 r 


Les trois opérateurs H, L 2 et L z commutent deux à deux, ils admettent donc un système 
commun de fonctions propres. Soit T(r,0,cp) cette fonction commune, on écrit : 

HT = ET 
L 2 T = ^ 2 £(^ + 1)t 
L Z T = hm i x ¥ 

Comme L 2 n'agit que par ses variables angulaires (0,cp), on peut faire appel à la méthode de 
séparation des variables en posant : 

T(r,0,cp) = R(r)Y to (0,cp) 


où R(r) est une fonction radiale à déterminer et Y^ m (0,cp) les harmoniques sphériques 
précédemment connues (voir chapitre III). Le remplacement de T dans l'équation de 
Schrôdinger donne : 


_^_J__5 2 5_ h 2 t(l + 1) _ e^_ 
2m e r 2 ôr ôr 2m e r 2 r 


ERY, m 


D'où l'équation radiale de l'atome d'hydrogène : 

__ft^J__ô r 2 8_ + h 2 £(£ + !) 
2m e r 2 ôr * ôr 2m e r 2 



R(r) = ER(r) 


qui s"écrit encore : 


1 d 2 

d 

y+i) , 

fAvljl + 

f 2m e E^l 

r 2 dr 

dr " 

■ ~ 

[* 2 h 

l h 2 )_ 
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Dans cette équation, on identifie l'expression du rayon de la première orbite de Bohr : 

h 2 

a 0 = y = 0.53A 0 

m e e 


Lorsqu'on impose aux fonctions radiales R(r) d'être normalisables et finies dans l'intervalle 
re[0,œ[ l'énergie E ne peut prendre que certaines valeurs. Dans ces conditions, la 
quantification de l'énergie est assurée par un nombre entier n appelé nombre quantique 
principale : 


E n = 



n = 1, 2, 3, . . ., oo 


En remplaçant E n par son expression dans l'équation radiale, nous obtenons : 


J_d_ 2 d_ 

r 2 dr dr 



R^(r) = 0 


On voit que les fonctions radiales dépendent des valeurs des nombres quantiques n et i. 
Ces deux nombres quantiques sont liés par la condition : 

0 < i < n-1 


Pour permettre à l'équation d'être sans dimension faisons le changement de variable 
suivant : 


il vient : 



na o 


1 d 2 d t(£ + l) 2n 

d^ 


R n ^) = o 


Les fonctions R nf (lQ qui satisfont cette équation différentielle sont de la forme : 

R n( W=-N 11( èe- 5,2 L 2 n , + + ( 1 ft) 

N n ^ est le facteur de normalisation ayant pour expression : 


m 3 ' 2 

' (n-^-l)î ' 

v na o J 

_ 2n((n + £)î ) 3 _ 


et dénotent les polynômes associés de Laguerre qui sont d'ordre 21 + 1 et de degré 

n + i en Un polynôme de Laguerre d'ordre k est donné par : 


Son associé est simplement : 


L k fe) = e^(^) 


L ïft)=^r(L k ) 
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Dans le tableau suivant, on donne les premiers polynômes L 2 ^ 1 (^) et les premières 
fonctions radiales R ne (r) de l'atome d'hydrogène avec c, = 2r/na 0 : 


n 

t 

Ctfe) 

Rnf(r) 

Partie 
radiale OA 

1 

0 

Life)=-1 


ls 

2 

0 

L‘fe)=^l+2Ç 


2s 

2 

1 

L 3 3 fé) = -6 


2p 


Les fonctions radiales R nf (r) de l'atome d'hydrogène sont des fonctions réelles et forment 
un système orthonormé complet : 

{r f | R n'f >)*= £ R n< (r)R„ r (r)r 2 dr = 8 nn .8„, 


Finalement, les états de l'atome d'hydrogène sont décrits par des fonctions d'onde 
dépendant des trois nombres quantiques : n nombre quantique principal, t nombre 
quantique orbital et m nombre quantique magnétique : 


T nAn (r,0,cp) = R n ,(r)Y, m (0,(p) 


n = 1, 2, oo 
0<£ <n-l 
-£<m <+£ 


De telles fonctions d'onde décrivent les diverses orbitales atomiques de l'atome: ls, 2s, 2p, 
etc. Les expressions explicites des premières orbitales sont résumées dans le tableau 
suivant: 


n 

l 

m 

Notation 

'En/m 

1 

0 

0 

ls 

JeU) 

3/2 

e 

-r/a 0 



0 

0 



f 2 | 

f j 

' | -r/2a 0 



V8ti v a o 

J 1 

^ 2a 0 J 

2 

1 

0 

2po 

V87r( a o 

n 

L 

e- r/2a ° COS0 

2 

1 

+1 

2p±i 

1 f 1 
4 VE v a o 

n 

& 

| e _r/2a ° sin 0 e ±itp 


V-2° Parité des orbitales atomiques 

Lors de l'opération inversion des coordonnées de position, les orbitales atomiques 
'P n£m (r,0,cp) = R n ^(r)Y to (0,cp) peuvent changer de signe. Etant donné que la coordonnée r 
reste inchangé, la parité des orbitales atomiques est celle des harmoniques sphériques 
Y^ m (0,cp) c'est-à-dire : 


^ntm h, 9, 9) = (-1) ' e 'Enfin (r, 71 - 0, 7T + Cp) 
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Autrement dit, les orbitales de symétrie i - 0, 2, 4, . . . (ns, nd, ng, . . .) sont paires et celles de 
symétrie ^ = 1,3,5,... (np, nf, nh, ...) sont impaires. Cette propriété de parité permet des 
simplifications considérables dans les calculs des intégrales. En effet, pour calculer par 
exemple la matrice (x) associée à l'opérateur position x = x dans la base des orbitales 
atomiques ls, 2s, 2p z : 

f (ls|x|ls) (ls|x|2s) (ls|x|2p z ) ^ 

(x)= (2s|x|ls) (2s|x|2s) (2s|x|2p z ) 

v (2p z |x|ls) (2p z |x|2s) (2p z |x|2p z ) y 


l'utilisation de la parité permet de préciser les intégrales non milles. Comme x est impaire, 
les intégrales non milles sont celles qui font intervenir les orbitales de parité différentes. 
D'où : 


(x) = 


0 

0 


l( 2 PzH ls ) 


0 

0 

(2p z |x|2s) 


( 1 s H 2 p z )" 

(2s|x|2p z ) 


0 


Comme la position x et les orbitales atomiques ls, 2s et 2p z sont réelles, la matrice est 
symétrique, c'est-à-dire (ls|x|2p z ) = (2p z |x|ls) et (2s|x|2p z ) = (2p z |x|2s). Il demeure enfin deux 
intégrales à calculer au lieu de neufs au départ. 


V-3° Spectre énergétique 


Les niveaux d'énergies des orbitales atomiques 4 / n/m de l'atome d'hydrogène sont 


donnés par : 



(eV) 


L'indépendance de E n des deux nombres quantiques t et m entraîne à une dégénérescence 
des états. Pour un niveau n donné, comme m prend [21 + \) valeurs possibles et l est 
comprise entre 0 et n-1, le degré de dégénérescence est donné par la suite arithmétique : 

n-l 

2](2^ + l) = l + 3 + --- + (2n-l) = n 2 

£=0 

Ainsi, à un niveau n donné correspond n 2 orbitales atomiques dégénérées, comme illustré 
par le diagramme d'énergie ci-après. Tous les 
niveaux d'énergie (l<n<oo) de l'atome 
possèdent des énergies négatives E n <0, ils 
correspondent donc à des états liés. Si l'atome 
est isolé, l'unique électron occupe le niveau 
de plus basse énergie (Ei=-13.60eV) qui 
correspond à l'état fondamental 4L S . Sous 
l'effet d'une radiation extérieure d'énergie hv, 
l'électron peut occuper les niveaux d'énergie 
supérieure (E 2 =-3.4 eY, E 3 =-1.5 eY, ...) qui 
correspondent aux états excités de l'atome. 

La dépendance de l'énergie en 1/n 2 implique que les niveaux excités se rapprochent 


4p 

'“2p" 
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davantage à raison que n augmente. A la l im ite n— »oo, l'énergie s'annule (Eoo=0), l'électron 
est envoyé à l'infini et l'atome se trouve dans un état ionisé (H — » H + + e). Cette ionisation 
nécessite une énergie | E œ - Ei | = 13.60eV appelée énergie d'ionisation ou potentiel 
d'ionisation de H. L'électron éjecté à l'infini (n— »oo) qui se trouve dans un état non lié (état 
de continuum) est une particule libre où elle n'est soumise à aucun champ de potentiel, son 
énergie est continue. 


V-4° Système des unités atomiques : 

Trois systèmes d'unités communément employées en mécanique quantique : 
systèmes MKS, CGS et u.a. Le système des unités atomiques (u.a) est un système qui a été 
introduit afin d'apporter des simplifications des valeurs des constantes physiques 
fondamentales telles que fi (constante de Planck), m g (masse de l'électron) et e (charge de 
l'électron) : 

fi = 1 u.a m e = 1 u.a e = 1 u.a 


De ce fait l'énergie de l'atome d'hydrogène E n et le rayon de la première orbite de Bohr a 0 
prennent des valeurs simplifiées. Ces deux quantités ont pour expressions dans les trois 
systèmes des unités : 


Système MKS 
1 


\{A%& Q Ÿ fi 1 )2n 
fi ^ 

a 0=( 4 7T8 0 ) y ( m ) 

m e z 


— j (joule) 


Système CGS 



Système u.a 

E n = ~~~p ( hartree ) 

a 0 = 1 ( bohr) 


Dans le système u.a, les deux nouvelles unités "le hartree" pour l'énergie et "le bohr" pour la 
longueur s'obtiennent de la manière suivante : 

1 hartree = lu.a = m e e 4 /ft 2 = 27.21 eY 
lbohr = lu.a = fi 2 ! m e e 2 = 0.529 A° 


VI- Méthode des variations 

VI-1° Principe variationnel 

Jusqu'ici nous n'avons considéré que le cas des systèmes simples où l'on a pu 
déterminer les solutions de l'équation de Schrôdinger de façon exacte. La majorité des 
problèmes en mécanique quantique ne peuvent être traité que dans une façon approchée au 
moyen des méthodes d'approximation. La méthode des variations est une méthode 
largement utilisée pour résoudre l'équation de Schrôdinger HT 7 = ET où il est impossible 
de connaître ses solutions de façon exacte. Pour tout problème, l'inconnu est toujours la 
fonction d'onde T. Si l'on part d'une fonction d'onde quelconque (appelée fonction 
d'essaie T essai ), on obtient pas l'énergie exacte Eo du système, mais une énergie moyenne 
que nous notons W : 
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La méthode des variations consiste à faire varier 'F essai de telle manière que l'énergie W soit 
la plus proche possible de Eo. Pour résoudre ce problème, on doit faire appel au théorème 
variationnel suivant. 


Théorème : 


"L'énergie moyenne W est toujours une limite supérieure de 
l'énergie exacte Eo de l'état fondamental" : 

(Y M ^|H|Y es aai ) 


En - W 


E n S 






ai) 


Pour démontrer ce théorème, on part de l'hypothèse que l'opérateur H admet une famille 
de fonctions propres {HL, 'Pi, , . ., 'En} et valeurs propres {Eo, Ei, . . ., E n }. Supposons que ces 
valeurs propres sont classées par ordre croissant ▲ 

Eo < Ei <,..., < E n . Dans ce cas, on a : E n — j 

H v E i =E i v E i _J 


Comme H est hermitique, les fonctions propres 
{ 'Ej } forment une base orthonormée complète, 

(Vj l 'EjJ) = Sÿ, on peut donc développer la fonction 
d'essai en série des 'Ej : 


El 


, W 

E 0 '' 

solutions solution 

exactes approchée 


Cessai = Zc^ 
i=0 

Pour simplifier, supposons que E essai est normalisée pour pouvoir écrire : 

W = fessai |H| Cessai) = |h|] Z c j'Ej ^ 

= ZZ c *c j ( v E i |h| 'E- \ = ZZc*c j ( v E i |E j | v E.\ 
i j 1 i j 1 

= EIc-cE^I*) = ZZc'cES 
i j i j 

Etant donné que les fonctions de base sont orthonormées, les seules termes non nuis sont 
ceux qui correspondent à i = j . D'où : 

W = ZcjCjEj = Z|ci| 2 E. 

Soustrayons membre à membre l'énergie exacte Eo, l'équation devient : 
w-E 0 = z|c,| 2 (e, -e 0 ) 

Par hypothèse Eo est l'énergie de l'état fondamental, donc Ej-Eq^ 0, et le module des 
coefficients est toujours positif |cj| 2 > 0. Ainsi, le membre de droit est toujours positif ou nul 
et le théorème est démontré : 

E„<W o E 0 < (ŸejHl*^) 

La valeur moyenne associée à la fonction d'essaie est toujours supérieure à l'énergie 


- 85 - 


Mécanique Quantique S MC (S 4 )~ Tès 


CHapitre tV 


exacte W > E 0 , l'égalité W = E 0 n'étant réalisée que si et seulement si la fonction d'essai est 
identique à la fonction d'onde exacte du système ( v ï / essai = Ÿ Q ). Le théorème qu'on vient 
d'énoncer traduit donc le principe variationnel où il est à la base du procédure de calcul des 
fonctions d'onde inconnues. 

En pratique, on choisit une famille de fonction d'essai T a dépendant d'un paramètre 
a ou d'un certain nombre de paramètres {ai, ai, ...} appelés paramètres variationnels. On 
calcule l'énergie moyenne W associée : 


W(a) = 


('i'.lfth.) 

bal '•'«) 


La variation de W(a) en fonction de a peut se 
traduire schématiquement par le graphique ci- 
contre. Pour que W(a) soit la plus basse, c'est-à- 
dire la plus proche de Eo, on dérive l'énergie par 
rapport à a : 


W(a op t) 



ÔW 

da 


Cette équation constitue la condition de minimisation de l'énergie, elle fournit la valeur 
optimale du paramètre variationnel (a op t) et l'énergie minimale W(a op t) recherchée. On 
espère alors que v P(a opt ) sera une bonne représentation de l'état fondamental 'Lo, et que 

W(a op t) sera une bonne approximation de son énergie Eo. Le choix de fonctions d'essai est 
important pour accéder à des solutions acceptables. Ce choix se base sur la nature du 
système physique à traiter et aussi sur les conditions de bonne fonction d'onde. 


VI-2° Application à l'oscillateur harmonique : 

On se propose de rechercher par la méthode des variations la fonction d'onde de 
l'état fondamental de l'oscillateur harmonique. Supposons que nous ne pouvons pas 
résoudre cet exemple par une méthode exacte. Pour l'oscillateur harmonique, on sait que la 
fonction d'onde 'F(x) de l'état fondamental doit être une fonction symétrique par rapport au 
centre x=0 : 

¥(x) = T(-x) 

et qui s'annule aux bornes de x qui peuvent être finies ou infinies : 
xe [-a,+a] -> ¥(+ 3 ) = 'E(-a) = 0 

xe ]-oo,+oo[ — » ¥(-00) = ¥(+00) = 0 

Nous pouvons alors choisir deux fonctions d'essai possibles : 

T^X) = COSGCX -^ < X < 

X ¥ 2 (x) =e ~ ax2/2 -oo<x<+oo 


où a est le paramètre variationnel (on peut vérifier que les deux fonctions s'annulent aux 
bornes correspondant). L'énergie moyenne associée à une fonction d'essai est : 


( v e|h| v e) 

( X F| V P) 


s n 2 d 2 i, 2 

H = - + -kx 2 

2m dx 2 2 
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En remplaçant chaque fonction d'essai dans W, on obtient les deux résultats suivants : 

0 1 


T.(x) -> W,(a) = a 2 +k — 

1 1 ^ I 24 4 J a 2 


T 2 (x) 

Les dérivations de W(a) donnent 
SW, 


W 2 (a) = a+ — k a 2 — 


— " a k 1 ‘ 0 

12 2 a 3 


Sa 

SW k_ 

Sa 4m 4a 2 


= 0 


_ Vkrn 
h 

_ Vkm 
h 


La substitution de a t dans W donne l'énergie minimale dans chaque cas : 


W, (min) = | — - - ?iJ— = 0.57 h J— = 0.57 hv 


W 2 (min) 


= *,11 = 0.5 

2 Vm 


La comparaison à la valeur exacte (obtenue par la résolution directe de l'équation de 
l'oscillateur harmonique, voir plus haut ) correspondant à l'état fondamental : 

E q = 0.5 hv 

révèle que la fonction d'essai ^ (x) = cos ax donne une énergie Wi > Eo avec une erreur 
relative de 14%. Par contre, avec la fonction ^(x) = e _ax “ /2 nous obtenons la valeur exacte 
de l'énergie Wi = Eo, car cette fonction a exactement la même forme que la solution exacte 
de l'état fondamental ( voir paragraphe III). Avec ces deux exemples, on voit que le principe 
variationnel est bien vérifié (Eo < Wl) : 

"Plus la fonction d'essaie est proche de la vraie fonction d'onde plus 
l'estimation de l'énergie est meilleure" 


VI-3° Application à un atome hydrogénoïde : 

Un atome hydrogénoïde est un système constitué d'un seul électron gravitant autour 
d'un noyau de charge +Ze (H, He + , Li 2+ , . . .). L'hamiltonien d'un tel système est : 

H - %2 1 8 r 2 8 \ L 2 Ze 2 

2m e r 2 dr ôr 2m e r 2 r 


qui dans le système u.a prend la forme simplifiée suivante : 

- 110,3 L 2 Z 

H = r 2 — + — 

2 r 2 dr dr 2r 2 r 

Pour l'état fondamental ls, les solutions exactes sont : 




r y / 

W 


c ls - 
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Pour l'application du principe variationnel, on doit choisir une fonction d'onde 
physiquement acceptable surtout d'être finie dans l'intervalle [0,+œ[. Deux fonctions 
d'essaie sont alors possibles : 


1 F|M = (Ÿ) 3 ' 


^ 2 (X) = 



où a est le paramètre variationnel à optimiser dans chaque cas. Comme les deux fonctions 
sont normalisées, l'énergie s'écrit : 


W= ( v F|h| v F) avec H = — l —— r 2 — - - 

' 11 ' 2r 2 dr dr r 


Remarquons que nous avons réécrit l'opérateur H sans la partie angulaire (contenue 
dans L 2 ), ceci est valable lorsque les fonctions d'onde à appliquer sont indépendantes des 
angles (0,cp). La substitution de 4^ (r) et ensuite v P 2 (r) dans W conduit aux résultats 


suivants : 

W 1 (a) = |a- A /|zV^ 

Les dérivations donnent : 


oW, 

s(V^) 

6W 2 _ 
da 


a-Z = 0 


Z = 0 


W 2 (a) = ^a 2 -Za 


a op t=^Z“=0.2829Z z 
ci t = Z 


La substitution de a opt dans W donne l'énergie minimale dans chaque cas : 

Wj (min) = -^Z 2 = -0.4244Z 2 
W 2 (min) = -iz 2 =-0.5000 Z 2 

La comparaison à la valeur exacte correspondant à l'état fondamental, E 0 = -0.5 Z 2 , révèle 
que la fonction d'essai Tjfr) donne une énergie Wi > Eo avec une erreur relative de 15%. Par 
contre, avec la fonction H^fr) nous obtenons l'énergie exacte Wi = Eo. La fonction d'essai 
T 2 (r) est alors le meilleur choix pour représenter l'état fondamental ls de l'hydrogénoïde. 
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Université Sidi Mohamed Ben Abdellah 
Faculté des Sciences D. Mehraz Fès 
Département de Chimie 


SMC-S4 

Complément de Physique 
TD de Mécanique Quantique 

Exercice 1° : On donne les fonctions d'onde suivantes définies dans l'intervalle x e [0,oo[ : 
'F 1 ( x ) = ax 2 ^(x^e ” 2 x F 3 (x) = e“ ax T 4 (x) = e _ax2 

a étant un paramètre défini réel et positif. 

(a) Lesquelles vérifient les conditions de bonnes fonctions d'onde? Pour celles qui ne les 
sont pas, indiquer la raison. 

(b) Parmi les bonnes fonctions d'onde, indiquer celles qui sont fonctions propres de 

l'opérateur Â = d 2 /dx 2 . A quelle condition (x ) = e -ax est-elle fonction propre de 

l'opérateur B = -d 2 /dx 2 + |3x 2 ? 

(c) Parmi les opérateurs suivants, indiquer ceux qui sont linéaires : 

Â = sin B = d/dx C = /~ D = d 2 /dx 2 

Exercice 2° : En utilisant le concept de fonctions paire et impaire, assurez-vous par simple 

inspection lesquelles des intégrales suivantes sont identiquement milles : 

(a) £siiupcos cp dcp (b) j^sincpcos (p (c) j^xcosxdx 

(d) 2 y dx dy (e) xe _x2 dx (f) J+“x 2 e _x2 dx 

Exercice 3° : On donne les deux matrices A = j et B = ^ ^ j . 

(a) La matrice A est-elle hermitique, inversible et/ou diagonalisable ? 

(b) A quelle condition la matrice B est-elle hermitique? Trouver ses valeurs propres et ses 
vecteurs propres. Les valeurs propres sont-elles réelles et les vecteurs propres sont-ils 
orthogonaux? 

(c) Calculer AB et BA . Le produit des deux matrices est-il commutatif? 

(d) Calculer l'adjoint de AB sachant que a est réel. Vérifier que (AB) + = B + A + . 


Exercice 4°: On donne les deux fonctions d'onde v P 0 (x) = N 0 e _ax /2 et v t'i(x) = Nixe -ax /2 
définies dans l'intervalle xe]-oo,+oo[. Normaliser les deux fonctions à l'unité et montrer 
qu'elles sont orthogonales. On donne : 


F 


2n-ax 2 , _ lx3x....x(2n — l) 


Exercice 5° : Montrer que le produit scalaire de deux fonctions d'onde vérifie la relation 
({D i |q'j) = ( v Fj|<D i ) . Soit |O i ) = |Â'F i ^ un ket de l'espace direct S des états. En utilisant son 
dual (OileS*, démontrer la relation ('f' i |Â + | x Fj} = ('Fj|Â| v î'i}* entre l'opérateur Â et son 
adjoint Â + . 

Exercice 6° : Dans l'hypothèse d'une base orthonormée formée des trois kets jcpj), |cp 2 ) et 
|cp 3 ) de l'espace direct, on exprime un ket |t) du même espace par le développement 
linéaire |’F}=i|q> 1 ) + (1 - i)|q> 2 ) + V2|cp 3 ) . 


- 89 - 


Mécanique Quantique SMC (S4)~ Tès 


Travatoc r Diriqés 


Sachant qu'à chaque ket |cpj) de la base est associé son bra (cpj | : 

- quel est le dual de | v k) ? 

- quelle est la représentation matricielle de l'k) et (W | dans la base considérée? 

- calculer le produit scalaire (H / |4 / ) et en déduire le facteur de normalisation de ¥ . 

Exercice 7° : Soit Â un opérateur hermitique défini dans une base {9 J. A quelle 

condition on a l'égalité (cp i |Â| cp j ) = ^cp j |Â| cp j ) ? Quelles sont les conditions requises pour 
qu'un élément de matrice (cpi|Â|cpj^ soit nul? 

Exercice 8° : Montrer les identités mathématiques suivantes : 

[â,b]=-[b,â] [â,b+c]=[â,b]+[â,c] 

[aâ,b]=a.[â,b] [â,bc]=[â,b]c+b[â,c] 

Développer les deux commutateurs [ Â , B 2 ] et [ Â 2 , B ] . 

Exercice 9° : Pour une particule de masse m se mouvant dans un potentiel scalaire V(x) , 
calculer les commutateurs suivants : 

[x,p x ] [x,p 2 ] [x, V(x)] [p x ,V(x)] [x,h] [H,xp x ] 

Dans le cas d'un potentiel dépendant du temps V(x,t) que vaut [t, h]. 

Exercice 10° : Soient Â et B deux opérateurs hermitiques non dégénérés. 

(a) Montrer que les valeurs propres d'un opérateur hermitique sont réelles et que ses 
fonctions propres sont orthogonales. 

(b) Montrer que si Â et B commutent, ils admettent un système commun de fonctions 
propres. Réciproquement, montrer que s’ils possèdent un système commun de fonctions 
propres, ils commutent. 

(c) A quelle condition on a l'égalité (Â + B) 2 = Â 2 + 2ÂB + B 2 ? 

Exercice 11° : Considérons une particule de masse m mobile dans un référentiel tri-direct 

Oxyz . 

(a) Etablir les expressions quantiques des composantes L x , L y et L z du moment cinétique 
en coordonnées cartésiennes. Utiliser le système u.a pour exprimer ces composantes. 

(b) Calculer les commutateurs [l p ,lJ et [l 2 ,lJ avec p,q = x,y ou z. Commenter. 

(c) Pourquoi les opérateurs L 2 , L x , L y et L z doivent être hermitiques? 

(d) Pourquoi les opérateurs escalateurs L + =L x +iL y et L_=L x -iL y ne sont pas hermitiques 
alors que leurs produits L + L_ et L L + les sont ? 

Exercice 12° : En mécanique quantique l'opérateur carré L 2 du moment cinétique orbital 
appliqué sur une fonction angulaire Y(0,cp) vérifie l'équation aux valeurs propres 
suivante : 

L 2 Y(0,cp) = -h 2 J-A sin0 A + — L — ^_| Y (0,<p) = h 2 £(£ + l)Y(0,cp) 

50 sin 2 0 5cp 2 J 

Montrer que si on pose Y(0,cp) = F(0) <D(cp) l'équation précédente se sépare en deux 
équations différentielles indépendantes. Déterminer l'expression normalisée de <D(cp) et 
vérifier que l'opérateur L z est quantifié et non dégénéré. 
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Exercice 13° : Les harmoniques sphériques sont définies par : 
N 


Y lm ( 0 , 9 ) = Pp (cos 0) e im(p 
Q2n 


\ (2/ + l)(/-:m) il : 

,m "L 2(*+M>l J 


En posant £ = cos 0 , les polynômes associés de Legendre sont : 


pH© = <=i> 


,«H ,a d *t 


T <x-vr 


(a) Vérifier que les fonctions Y ^ 0 sont réelles. En déduire leur parité et celle des 
harmoniques sphériques Y l m en général. 

(b) Trouver les expressions des premières harmoniques sphériques Y 00/ Y 10 et Y 1±1 . 
Donner leurs expressions en coordonnées cartésiennes. Sous cette forme, justifier la parité 
de chacune de ces harmoniques sphériques. 

(c) Représenter les surfaces de densités associées à Y 0;0 et Y 10 . Donner des exemples où 
Y 0 ,o et Yj 0 peuvent intervenir. 


Exercice 14° : On donne les harmoniques sphériques Y(0,cp) = N et Y'(0,cp) = N'cos0 

fonctions propres de L 2 et L z : 

(a) Normaliser les deux fonctions et montrer qu'elles sont orthogonales. 

(b) Attribuer les nombres quantiques l et m associés à chacune des fonctions. 

(c) Partant de Y' (0, 9 ) , trouver les expressions des harmoniques sphériques ayant même l 
mais m différents. 

Exercice 15° : Etablir l'équation de Schrôdinger pour chacun des systèmes stationnaires 
suivants : 

(a) rotateur rigide sphérique (c) atome He et son isoélectronique H 

(b) oscillateur harmonique (d) ion moléculaire Yp 

Dans chaque cas, préciser les variables impliquées dans la fonction d'onde et les 
éléments d'intégration impliquées dans la condition de normalisation. 

Exercice 16° : On se propose de déterminer les états stationnaires d'une particule de 
masse m plongée dan un puits de potentiel défini par : 

V(x) = 0 pour 0 < x < L 

V(x) = oo pour x < 0 et x > L 

(a) Quelles sont les parties d'espace interdites à la particule? En déduire les conditions aux 
limites que doit satisfaire la fonction d'onde x F(x) du système. 

(b) Ecrire l'équation de Schrôdinger de tel système en introduisant le nombre d'onde 
k = (2mE//z 2 ) 1/2 . Trouver les solutions de l'équation de Schrôdinger T n (x) et E n et montrer 
qu'elles sont sous la dépendance d'un nombre quantique n à définir. Pourquoi le cas n = () 
est-il impossible? Montrer que deux fonctions d'onde T n et Y n . différentes sont 
orthogonales. 

Exercice 17° : On donne la fonction d'onde T(x) =(l/ 47 i) l/ 4 ( 2 x 2 - 1 ) e -x /2 représentant un 
état quantique de l'oscillateur harmonique. 

(a) Vérifier que T(x) est normée. 

(b) T possède-t-il des points nodaux? Trouver la distance de maximum de densité x max et 
la distance moyenne <x>. Commenter. 
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(c) Tracer la densité de probabilité D(x) = |T(x)| 2 en fonction de x . 


Exercice 18° : Considérons un atome d'hydrogène exposé à un champ magnétique B 
statique et uniforme dirigé suivant la direction Oz . Au moment cinétique L de l'électron 
est associé le moment magnétique fl propre qui lui est colinéaire et proportionnel 
jI = -(e/2m e )L. L'énergie d'interaction entre ce moment et le champ magnétique B est 
W = — pLB . Cette interaction donne lieu au phénomène connu sous le nom d'effet Zeeman. 

(a) En utilisant le système des unités CGS , rappeler l'expression de l'hamiltonien H° de 
l'atome d'hydrogène isolé et ses énergies propres E° associées aux fonctions propres 

'W 

(b) Ecrire l'hamiltonien H de l'atome soumis au champ magnétique B. Pourquoi 'P nAn 
sont-elles également fonctions propres de H. De quels nombres quantiques dépend 
l'énergie des états stationnaires de l'atome soumis au champ magnétique? 

(c) Montrer qu'il n'y a pas d'effet Zeeman sur l'état ls , alors que sur l'état 2 p cet effet 
provoque une levée de dégénérescence. Représenter les résultats sous forme d'un 
diagramme illustratif. Interpréter et commenter le spectre de raies qui en résulte. 

Exercice 19° : Soit Â un opérateur hermitique dont sa valeur moyenne observable dans 
un état 'T peut dépendre explicitement du temps : 


(a) Montrer que la variation temporelle de (A) obéit à la relation d'Ehrenfest suivante: 


(b) Il s'agit d'appliquer cette relation au cas des opérateurs position x = x et impulsion p x . 
Dans le cas d'une particule de masse m soumise à un champ de potentiel V(x), montrer 
que : 


(A) = ( x i'|Â| v ï') = J 'F*(x,t)Â't'(x,t)dx 





et 



dt 


m 


On utilisera les relations [p x , x| = 2— p x et [p x ,V(x)] = — — sans les démontrer. 

i i dx 

(c) Etablir la relation de Heisenberg dans le cas de l'observable produit xp x : 
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Problèmes de synthèse 


Problème I - Puits de potentiel 

On donne la fonction d'onde de l'état fondamental d'une particule piégée dans un puits 
de potentiel V(x) = 0 de largeur L : 



1°- Calculer les valeurs moyennes (x), (p x ) et (T x ). Pourquoi (T x ) donne exactement la 
valeur propre E n ? 

2°- En statistique, la variance a A d'une grandeur A est définie par a 2 A =< A 2 > - < A > 2 , 
elle permet la détermination de la déviation standard : 

AA = -\J< A 2 > - < A > 2 

qui est une mesure de la dispersion d'une distribution statistique autour de la valeur 
moyenne. Cette déviation sera considérée ici comme une incertitude sur la mesure de la 
grandeur A . A l'aide de la fonction (x) on peut montrer que : 


>= L2 h^ 

V 3 2n 2 n 7 


et 


h 2 n 2 n 2 

~ir~ 


Déterminer Ax et Ap x ainsi que le produit AxAp x . La relation que vous obtenez dépend-t- 
elle de la largeur L du puits? Interpréter alors la relation d'incertitude de Heisenberg et 
la discuter pour les trois plus bas états stationnaires. 

3°- Supposons qu'au temps t = 0 la particule était préparée dans un état décrit par : 


m 1 f . nx . 27tx N ) 
x =— ;= sin — + sm 

Æl L L J 


Écrivez la fonction d'onde de l'atome à un temps t > 0 . Déterminez la probabilité de 
trouver la particule au temps t avec : 


i) E = - 


ii) E = - 


iii) E = - 


Calculez la valeur moyenne de E . 


Problème II - Oscillateurs diatomiques 

Dans ce problème, on s'intéresse à l'étude du 

mouvement interne d'une molécule diatomique. Une q x 

molécule diatomique AB peut être assimilée à un >4 

système mécanique où les deux masses m A et m B r A r B 

vibrent harmoniquement autour de leur position r 

d'équilibre. La liaison chimique de tel système est assimilée à un ressort de constante de 

raideur ou constante de force k. 
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Soit r 0 la distance d'équilibre de la molécule, r = r A + r B la distance instantanée, r A et r B 
les séparations des deux masses relatives au centre de masse G . Le déplacement effectif 
des deux masses (élongation ou rétrécissement de la liaison) est mesurée par la quantité 
x = r-r 0 . 

1°- Dans un premier temps, on se place dans un référentiel d'origine O quelconque. 

(a) Démontrer que l'énergie cinétique du système est : 

T = ^M G V£ + i^v 2 

où M g =m A +m B est la masse totale, V G la vitesse du centre de masse G, v la vitesse 
relative des deux masses et n la masse réduite. 

(b) Que devient cette énergie cinétique si on se place dans le référentiel barycentrique 
supposé fixe? 

(c) Déterminer l'énergie potentielle du système. 


2°- On se propose de trouver une forme simplifiée de l'équation de Schrôdinger. 

(a) Ecrire l'équation de Schrôdinger du système. 

(b) Réécrire cette équation en posant a = 2\ Æ/h 2 et (3 2 = \ik/h 2 . 

(c) Etablir les dimensions des deux paramètres a et p . 

(d) Proposer un changement de variable x en E, qui permet l'obtention de l'équation de 
Schrôdinger sans dimension. 

(e) Sachant que les énergies propres E u de l'oscillateur harmonique (u nombre quantique 
de vibration et v fréquence harmonique) sont : 


E U = (u + E )hv, v = /re- 


montrer que l'hamiltonien sans dimension prend la forme simplifiée suivante : 


H = + E 

de 


Quelles sont les valeurs propres s u associées? 

3°- Les fonctions propres de l'hamiltonien H sont de la forme : 

= H u ©e -^ /2 avec H u Œ) = f> u Ç u 


où HyiÇ) est un polynôme généralement de degré infini en E, . 

(a) Trouver l'équation différentielle que doivent vérifier ces polynômes. 

(b) On peut montrer que la condition pour que X F U © soit une solution acceptable, c'est-à- 
dire 11,^© soit de degré u fini, amène à la relation suivante : 

H u =2oH u _ 1 

où H,j est la dérivé première de H,, . Montrer alors que les polynômes d'Hermite 
vérifient la relation de récurrence suivante : 

Hy = 2^H U _ 1 - 2 uH u _ 2 

(c) Connaissant les premiers polynômes d'Hermite H 0 =l et H ( =2^, déterminer les 
expressions des quatre premières fonctions d'onde correspondant à l'état fondamental 
v = 0 et aux états excités u = 1, 2 et 3. Discuter la parité et la forme des solutions ainsi 
obtenues. 
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4°- On donne la fonction d'onde suivante décrivant un état quantique de l'oscillateur : 

¥© = N(4^ 4 - 12^ 2 + 3) e - ^ 2 1 2 

Par simple inspection attribuer le nombre quantique de vibration u associé. Confirmer 
votre réponse par un calcul approprié. 

5°- Applications : 

(a) Sachant que la molécule d'acide chlorhydrique HCl possède une constante de force 
k = 469N/m, calculer l'énergie des niveaux de vibration en joule. Quelle sera la 
fréquence du photon absorbée au cours d'une transition du niveau u = 0 vers le 
niveau o = 1? Quelle est le nombre d'onde (en cm" 1 ) et la longueur d'onde ( en A puis 
en nm) associé au même photon absorbé? De quelle région du spectre 
électromagnétique s'agit-il? 

(b) Comparer les constantes de force k HC1 et k DC1 des deux isotopes HCl et DC1 sachant 
que le nombre d'onde respectif est 2883 cm -1 et 2068 cm -1 . Commenter. 

(c) Comparer les constantes de force k CN et k CN +des systèmes diatomiques CN et CN + 
sachant que le nombre d'onde respectif est 2069 et 1580 cm" 1 . Commenter. 

On donne : 

masses atomiques : mn = 1.008, mD = 2.013, mci = 35.453, me = 12.011, wn = 14.007 

constante de Planck : h = 6.6262 10" 34 J.s 

nombre d'Avogadro : N a = 6.0231 10 23 

célérité de la lumière : c = 2.9979 10 8 m/s 

lnm = lO^m, 1Â = 10 _ï0 m, ldyn = 10~ 5 N, lerg = 10~ 7 } 


Problème III - Représentation matricielle de l'opérateur position 

On donne les deux fonctions d'onde unidimensionnelles décrivant trois états quantiques 
de l'oscillateur harmonique : 


^aW = 


JlT 


AJ 


e-P* /2 


'I'b(x) = 


4(3 J 


-p x 2/2 




2 _ ne -Px 2 / 2 


(2(3x -l)e 


1°- Par simple inspection, identifier le nombre quantique de vibration u attaché à chaque 
fonction d'onde. Quelles sont les énergies propres associées? 

2°- On s'intéresse ici à étudier l'opérateur position x = x dans la base ( x f' a , x f' b , x f' c ). 

(a) Quelle est la représentation matricielle (x) de cet opérateur dans la base considérée. 
Montrer que cette matrice est symétrique. 

(b) Par simple inspection des éléments matriciels de (x), préciser ceux qui sont nuis et 
ceux qui ne les sont pas. Pour ces derniers, déterminer leurs valeurs. 

(c) Calculer les valeurs propres et les fonctions propres de la matrice (x). En déduire la 
valeur de la position moyenne réellement observable? Dans quel état peut être 
mesurée cette position moyenne? Déterminer le poids de participation des fonctions 
de base ( 'f'a , ^ ^ ) dans cet état. 

3°- On suppose qu'une particule se trouve strictement dans l'état x F b (x) . 

(a) Dans quelle position moyenne se trouve la particule? 

(b) Quelle est sa position nodale et sa position la plus probable? 


On donne : 


x( 2n-l ) 


’ ^ [ x 2n+1 e dx = — 

2n+i JO 9 ..«+/ 
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Problème IV - Théorème du viriel 

On se propose ici d'étudier l'équilibre entre l'énergie cinétique et l'énergie potentielle de 
quelques systèmes quantiques dans le sens du théorème du viriel. Soit H l'hamiltonien 
d'un système à une dimension : 


L'énergie potentielle V(x) pour la plupart des problèmes quantiques est une fonction de 
puissances de la variable x , c'est-à-dire : 

V(x) = Àx n 

où n est un entier positif ou négatif et X une constante réelle. 

r. % 

1°- (a) Calculer le commutateur [H,xp x ] . Pour ce calcul utiliser les relations [p x , x] = 2— p x 

i 

et [p x ,V(x)] = — sans les démontrer, 
i dx 

(b) Que vaut la valeur moyenne ^ [H, xp x ] 'j du commutateur [H,xp x ] prise dans un état 

propre 'P de l'opérateur hermitique H (état stationnaire)? Déduire la relation du 
viriel pour un champ de potentiel V(x) quelconque : 

*>-(■£) 

(c) Que vaut cette relation dans le cas d'un champ de potentiel V(x) = Âx n . 

(d) Appliquer le théorème du viriel au cas des deux systèmes quantiques: oscillateur 
harmonique et atome d'hydrogène. 

2°- On se propose d'appliquer le théorème du viriel à l'état fondamental de l'atome 
d'hydrogène décrit par la fonction gaussienne approchée suivante : 

T ls =(2ot/7i) 3/4 e -ar2 

(a) Rappeler la forme de la fonction d'onde exacte de l'état fondamental ls de l'atome 
d'hydrogène. 

(b) Dans le cas de la fonction gaussienne et en utilisant le système des unités atomiques, 
montrer que : 


(c) Trouver la valeur de a qui devrait satisfaire le théorème du viriel. Exprimer l'énergie 
totale E de l'atome en fonction du paramètre a et en déduire sa valeur. Comparer à 
la valeur exacte E = -0.5 u.a . 

(d) Quelle est la valeur de a qui devrait être obtenue par la méthode de variation qui 
optimise l'énergie E? Commenter. 


Expression du Laplacien radiale : A = 1 d r 2 — ; intégrale utile : \' r 4 e z 
r dr dr J0 


8 {(2a) 5 
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Corrections 


Problème I - Puits de potentiel : 

On donne la fonction d'onde d'une particule piégée dans un puits V(x) =0 de largeur L : 


- , . 2 . f nnx à 

""HH-J 


n- 1,2,3, ... 


1°- La valeur moyenne d'une grandeur physique A qu'on peut mesurer dans un état 
® n (x) est donnée par : 

(® n |Â|® n ) J®*Â® n dx 
(®n|®n) J®*® n dx 

Le problème considéré ici est un problème à une dimension, donc l'élément d'intégration 
est dx. Par ailleurs, comme la fonction d'onde ® n (x) est donnée sous la forme normée, 
cette valeur moyenne devient : 

< A > = (® n |Â|®n)= J^nÂ^ndx 
En remplaçant ® n (x) par son expression, on obtient : 

2 rL . f n 


. 2 (-L . f nnx I r . f nnx ) 

< A > = — sin Asm dx 

L J o ( L J v I. 


Valeur moyenne de la position : 


2 r L 

< x > = — x sin 

L Jo 1 




Comme on a sin 2 0 = ^(1- cos 20), il vient : 

1 rL I" f 2nnx Y| 1 f L , 1 A f 2nnx V L 1 f L f 2nnx V 

L J o [ l L JJ L J o { h J 2 L J o l L J 

L'intégration par partie de la fonction x cos kx donne : 

rL ( 2nnx V [ L . f 2nx^| f L f f 2nnxV| 

xcos dx = xsin + cos 

J ° l L J 2nn ^ L J {2nnJ { L )\ 


- -ALi„f2ïîiV— 4— It - aM-t — - « 

2nn|_ ^ L J 2nn Ç L 2nnV2nn 2nn J 


D'où : 


L 

cx> = — 
2 


La valeur moyenne de la position de la particule dépend de la largeur L du puits, elle est 
indépendante du nombre quantique n c'est-à-dire indépendante de l'état quantique 
$ n W. 

Valeur moyenne de l'impulsion < p x > : 

2 fL f nnx à „ . ( nnx V h fL . f nnx V 9 , i nnx , 

<Px> = rJo sm [irJ PxSm u J dx = 


tïrm fL . f nnx f nnx à , 

= — — sin cos dx 

iL 2 0 l L J l L J 
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Comme sin 0cos 0 = ^sin 20 , il vient : 


hnn [ L f2imxY| L h . . .. 

<Px> *' — V cos - (-1 + 1) =0 

2iL 2 [_ 21171 l L JJ Q 2iL V 


D'où : 


<p x > = 0 


La valeur moyenne de la quantité de mouvement de la particule est nulle quelque soit 
l'état stationnaire fondamental (n=l) ou excité (n=2,3,...). 

Energie cinétique moyenne < T x > : 

< T x > =(o„|T x |o n ) avec T x O n =- - 3 ° n 
2m 5x 2 


ô 2 CD n [Y ô 2 f . nn ) [2 ô f nn ütï Ï2 

dit 2 ” l Sm_ L X J * VL &[l7 COS 17 X J ~ VL 

n 2 n 2 f[2.nn) nV 

' ~ IVr 1 ” L *J - 


D'où T x cD„ = ^ n n ® n et la valeur moyenne s'écrit : 

2mL 2 

m /» t H K H | \ h 71 n j | \ 

< T x > = (<t n I — |<D n ) = — (^n^n) 


2mL z 


2mL z 


Comme ® n (x) est normée, on aura : 


h 2 n 2 


On obtient exactement la valeur propre E n de l'hamiltonien H puisque d'une part cet 
hamiltonien est purement cinétique (car V(x) = 0) et d'autre part les fonctions ® n (x) sont 
fonctions propres du même hamiltonien : 


2°- L'erreur 4 A commise lors de la mesure d'une grandeur physique A est donnée par 
l'incertitude quadratique moyenne : 

AA = V< A 2 >-< A> 2 


Incertitudes Ax et Ap x : 
A x = V<x 2 >-<x > 2 


2 T 2 fl 1 

< x > = L“ — - 

U 2 nV, 


Ap x =V<Px >-<Px > 2 

D'où : 


<Px > = 0, 


L fl 2 

: = J-n : 

2mi v 3 


AxAp x = fi 


. h"n~ finn 

Ap ‘ = U T = T 
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On voit que le produit AxAp x ainsi obtenu est indépendant de la largeur du puits, par 
contre il dépend du nombre quantique n. 

On peut montrer facilement que ™ -1 est une quantité toujours supérieure à 1 
quelque soit le nombre quantique n. En effet : 

. n 2 n 2 


n n 
\~ 6 ~ 


> 1 


Vn 


Si n = 1, on a — = 1.64 > 1 . On retrouve donc la relation d'incertitude de Heisenberg : 
AxAp x > h 

L'incertitude est d'autant plus grande que n est grand, c'est-à-dire l'erreur augmente vers 
les états excités. En effet : 



n 

AxAp x 

état fondamental 

1 

0.8 h 

1 er état excité 

2 

2 A h 

2ème état excité 

3 

3.7 h 


3°- Supposons qu'au temps t = 0 la particule était préparée dans un état décrit par : 


v E(x,0) 


1 ( . 7tx . 2roO 

= Ær^ +sm_l 


L L 

Cet état est une combinaison linéaire normée de deux états : état fondamental 
O! = ^sin^- et premier état excité 0 2 = ^sin^^- : 

v P(x,O)=-^(0 1 +0 2 ) 

V2 

L'évolution temporelle d'une fonction d'onde est donnée par le produit : 

iE 

0 n (x,t) = O n (x,O)e n 

D'où l'expression de la fonction à l'instant t : 

( _iK,t ili 2 t à 

0! e n +0 2 e n 

J 

Pour pouvoir déterminez la probabilité de trouver la particule au temps t avec les 
énergies : 


i)E = 




ii) E = 


2/zV 


iii) E s 


9h n 


2mLr mL 2mL 2 

on doit connaître à quel état appartient chacune des énergies sachant que : 

h 2 1 Z 2 2 


E n = 


2m L 


On voit immédiatement q'on a : 


i)Ei - 


h 2 n 2 


ii) E 2 = - 


iü) E 3 = - 
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Pour une combinaison linéaire de fonctions d'onde 'P = Xci'Pj / la probabilité de trouver la 
particule dans l'état ^ d'énergie Ej est proportionnelle au carré du coefficient de 
participation : 

Ry =|cî | 2 

Pour notre fonction d'onde : 

iE.t 

1 ' 


on aura : 


'P(x,t)=-^r O ie n +$ 2 e 


i) E = Ei 


ii) E = E 2 — ^ Pvp — — 

iii) E = Ej — ^ P\p = 0 

Pour le calcule de la valeur moyenne de E, on utilise la relation : 

(3511) 

- m 


j / iEit iEît iEit iE2t 

( V F| V F) =—($> l e h +0 2 6 h 0i e h~ +0 2 e‘T 

= I[( <D ,| <Dl ) + ( <D2 | 1 . 2 ) + 2<<l. 1 |<l, 2 )] = lft + | + 0 ) = 1 


/ iEit iEît I _ I iEit iE2t\ 

(T) H l'P) = /0ie~lT +0 2 e ï| H 0ie~ + 0 2 e“T} 

/ i-Èil J I 1ÎEyt : V / ijjft j _ I iEît 

= / 0ie h H 0i ê~T J + (0 2 e iT H 0 2 e h~ 

i Et t I _ I IBat \ 

( 0ie n~ H 0 2 e n~ ) + (0 2 e a~|H|0ie’ 


iEit I iE2t \ / iE2t I iEit 

( , F|H|'F) = (0i | H |0i) + (0 2 | H |0 2 ) + (0i | H |0 2 ) le td e id \+ (0 2 | H |0i) ( e hT e /T 


(0i|H|0i) = (0i|E 1 |0 1 ) = Ei(0i|0i> = E 1 

(0 2 |h|0 2 ) = (0 2 |e 2 |0 2 ) = e 2 (0 2 |0 2 ) = e 2 
(0 1 |h|0 2 ) = (0 1 |e 2 |0 2 ) = e 2 (0 1 |0 2 ) = o 


D'où la valeur moyenne de l'énergie : 

^ -, tlW 2^27^2 5 % 2 % 2 

E = Ei+E 2 - + s 

2mL2 mL2 2 mL2 
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Problème II - Oscillateurs diatomiques : 

1°- (a) Dans un premier temps, on se place dans un référentiel d'origine O quelconque. 
Dans ce référentiel, le mouvement de déplacement des deux masses m A et m B est repéré 
instantanément par les coordonnées suivants : 
m G 

m . - m B 



r A et r B étant les vecteurs positions des deux masses m A et m B , R celui du centre de 
masse G (M G = m A + m B ) par rapport au point fixe O . La position R du centre de masse et 
sa vitesse R sont définies par : 


m A + m B 


m A +m B 


On pose : 

f A = v A / f B =V B , R=V g 

Si la vitesse de m A par rapport m B (dite vitesse relative) est v , alors on aura : 
v = v A -v B 

En résolvant simultanément les deux équations : 


Vr, =- 


et 


L'énergie cinétique du système est : 


T - - m A v A+- m B v B = ~ n 


1 . \ it 2 1 m A m B 2 

= -(m A +m B )VG+- 

z z m A + m B 

1 2 1 2 
= ~ m G V G "*""2 


avec n = mAmB étant la masse réduite et M G = m A + m B la masse totale au barycentre. 
m A + m B 

(b) Si on se place maintenant dans le référentiel barycentrique supposé fixe, les deux 
masses se repèrent par les vecteurs positions suivantes : 

_ JL — «—JL 

G 

m A m B 

\ f A r B * 

r 

en particulier on a R = 0 et R = 0 , le terme d'énergie cinétique du centre de masse est nul. 
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L'énergie cinétique devient simplement : 


L'étude du système à deux masses m A et m B se ramène à celle d'un système effectif 
constitué d'une seule masse p séparée du centre G par la quantité x = r - r 0 : 


m A < 





G< 


F 


► P 


(c) Comme le mouvement de déplacement s'effectue de façon harmonique, les deux forces 
de rappel F A et F B sont égales et proportionnelles au déplacement x = r - r 0 : 


F a =F b =F F = -kx 

Comme cette force dérive de l'énergie potentielle, on aura : 


F = 

~ V(x) = 

dx 

kf X xdx = -kx 2 
Jo 2 

Le potentiel de l'oscillateur harmonique est donc de type parabolique. 

2°- (a) L'hamiltonien de Schrôdinger correspond à l'énergie totale classique du système : 

E = — + — kx 2 

2p 2 

- 

H = C + I kx 2 = _Lii + I kx 2 

2 n 2 2fi dx 2 2 

L'équation de Schrôdinger s'écrit : 


r_L4 + ikx' 

^ 2p dx 2 2 

J'F(x) = E'F(x) 

dx 2 h 2 y 2 ) 

En posant a = 2pE//î 2 et p 2 = 

= pk / h 2 , l'équation devient 



d 2v F 7 

— y- + (a -p 2 x) 'F = 
dx“ 

0 

(b) Les dimensions de base pour exprimer toute grandeur physique sont : 

[M] : masse, [L] : longueur. 

[T] : temps 

La dimension de E est celle de l'énergie cinétique (T=^mv 2 ), de l'énergie potentielle 
( V = -jkx 2 ) ou bien de l'énergie de rayonnement ( E = hv ). Ainsi, on obtient : 

B 1 2 1 ( 
E = — mv =— m 
2 2 \ 

}) ^ 

E = [M] [L] 2 [T]- 2 

h = — = Et 

=> 

h = [M] [L] 2 [T]- 1 

k -l-A 

x 2 x 2 

=> 

k = [M] [T]- 2 

Les dimensions des paramètres a et P s'écrivent donc 


2pE 

=> 

a = [L] 2 

h 

=> 

P = [L]' 2 
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Les deux paramètres possèdent la même dimension [L] 2 . Sachant que la fonction d'onde 
¥(x) n'a aucune dimension, l'équation de Schrôdinger possède donc une dimension en 
[L]- 2 , inverse de longueur au carré. 

(c) Pour que cette équation soit sans dimension, on doit procéder à un changement de 
variable en écrivant : 

Ç = Ax 

Pour cela, il suffit de rechercher l'expression du facteur A tel que la nouvelle variable E, soit 
sans dimension : 


dx dx d^ 


4-xA A Ui-Ül 

dx 2 dx^dçj d £ 2 


Avec ce changement de variable, l'équation différentielle devient : 


a-p 2 J-|*F = 0 


dû 2 


_tL_Ëlp 
A, 2 A 4 ^ 


Ainsi par identification, on choisit le rapport P 2 / A 4 = % c'est-à-dire A = ^/P . D'où : 




6t A 2 P 


et l'équation de Schrôdinger sans dimension s'écrit 
d 2v P 

d^ 2 

(d) Sachant que les énergies propres sont : 

E, 


- Ç 2 ¥ = 0 


: u = (u+L)hv, v = 
le rapport a/p prend la valeur simplifiée suivante : 

“ = *IËe = )n lï = 2u+ 

p AV k AVk V 2> \\x 

La substitution dans l'équation différentielle donne : 

d 2 v F n 2 

+ (2u + l - Ç 2 )'F V =Q 

d^ 2 

L'équation aux valeurs propres correspondant s'écrit : 


- + * r v = (2u + p'î'v 


On voit que l'hamiltonien de l'oscillateur harmonique exprimé avec la variable E, sans 
dimension prend une expression très simplifiée : 
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Les valeurs propres associées à H sont simplement : 


£ u =2u + l 


3°- Les solutions de l'équation de Schrôdinger à rechercher sont de la forme : 


'Lu© * H 0 ©e- Ç,/Î , 


H„(9 = £a u Ç° 


(a) Les deux premières dérivées de V F Ü (£) sont : 




Après substitution dans l'équation sans dimension, on obtient l'équation différentielle que 
doivent vérifier les polynômes d'Hermite : 

Hj, -2^H U + 2oH u =0 

(b) En partant de la dérivée première de H u : 

H w =2uH u _j 

la dérivée seconde s'écrit : 

H’ = 2oH ü _ 1 = 2o(2oH u _ 2 ) = 4 u 2 H u _ 2 

La substitution dans l'équation différentielle donne la relation de récurrence recherchée : 
4o 2 H ü _ 2 -2^(2oH u _ 1 ) + 2 oH u =0 => = 2£H U _, -2oH u _ 2 

(c) Connaissant H 0 = 1 et H, = 2£, on détermine les polynômes associés à o = 1, 2 et 3 à 
partir de la relation de récurrence : 


Les fonctions d'onde des quatre premiers états quantiques de l'oscillateur harmonique 
s'écrivent : 


Notons que ces fonctions d'onde sont données sans le facteur de normalisation. On voit 
que ces fonctions sont alternativement paires et impaires : 


H 2 =2^H 1 -2H 0 => H 2 = 4E, 2 - 2 

H 3 =2^H 2 -4H 1 => H 3 = 8^ 3 - 122, 


u = 0 H 0 = 1 

o = l H 1= 2^ 

o = 2 H 2 =4Ç 2 -2 

o = 3 H 3 = 8^ 3 - 122, 


'Lq© « e ^ /2 
V L 1 © = T^~ en 
Y 2 © = (4^ 2 -2)e^ /2 
^3© = (8^ 3 -120e^ 2/2 


On doit ramener ces expressions à la variable d'origine x = © 2 £, : 


^o(x) = e- px2/2 
'Li(x) = 2Vpxe- px2/2 
Y 2 (x) = (4px 2 -2)e“ px2/2 
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'ï'o(x) pair 

'î'i(x) impair 

T 2 (x) pair 

T 3 (x) impair 

La parité suit donc le nombre quantique de vibration, c'est-à-dire (-l) u . 

Du point de vue forme, on voit que les fonctions d'onde forment deux séries : série paire 
(^ 211 ) et série impaire ('IVn) possédant respectivement des polynômes du même degré (2n 
et 2n+l) : 

v L 2n ( x ) = ( a 2 nX 2n +a 2 n- 2 x2n_2 + ”' + a o) e_PX ^ n = 0, 1, 2, ... 

^2n+i(x) = (a 2n+ iX 2n+1 +a 2n _iX 2n_1 +--- + a 1 x)e“ px2/2 n = 0,1,2, ... 
et le plus haut degré du polynôme correspond au nombre quantique u de l'état. 

4°- Soit la fonction d'onde décrivant un état quantique de l'oscillateur harmonique: 
ÿ© = (4^ 4 -12^ 2 +3) e“^ /2 

C'est une fonction paire de plus haut degré quatre, donc u = 4. Pour confirmer ce résultat, 
on agit l'hamiltonien sur cette fonction pour résoudre l'équation de Schrôdinger : 

^ = (_4^ 5 + 28^ 3 - 27 Ç)e“^ /2 ^ = (4^ 6 - 48Ç 4 + 1 1 1^ 2 - 27)e^ 2/2 

En substituant dans le premier membre de l'équation, on obtient : 

-■^-+^ 2 'P = (- 4£, 6 +48^ 4 -lll^ 2 +27)e _ ^ 2/2 +^ 2 (4^ 4 -12^ 2 +3) e - ^ 2/2 

« (36^ 4 -108^ 2 +27)e“ ç2/2 
= 9(4£, 4 -12^ 2 +3)e _ ^ 2/2 
« 9^ 

La valeur propre associée est : 

2u + 1 = 9 =4> u = 4 


5°- Applications : 

(a) L'énergie des niveaux de vibration de la molécule HCl s'obtient à partir : 

E u = (u + 2-) hv = (u + l)— — , k = 469N/m 
2n\n 

La masse réduite en uma (unité de masse atomique) est : 


P = 


m H m Cl 

m H +m a 


1.008x35.453 
1.008 + 35.453 


= 0.9801 uma 


Comme dans un atome gramme il y N a atomes réels, on doit diviser p par le nombre 
d'Avogadro N a . On convertit ensuite la masse obtenue en Kg : 


0.9801 
6.0231 10 23 


io - 3 


= 1.6272 10- 27 


kg 
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480 


| H 6.6262 1Q-- 
2 J 2x3.1416 il 1.6272 10~ 27 


|Ë 6 


soit : 


E u = 2.8639 10" /U (2u + 1) en J 
Les valeurs des énergies des niveaux u = 0 et u = 1 sont : 

E 0 = 2.8639 10 _20 J, E l = 8.5916 10“ 20 J 

La fréquence du photon qui doit être absorbée au cours de la transition u = 0 -> u = 1 est : 
v = El ~ E ° = 5-7277 10 = 8.6440 10 13 s -1 (ou Hz) 

h 6.6262 10“ 34 

Remarquons que cette même fréquence peut aussi être calculée à partir de la relation de 
Hook : 


v.J-JÏ 

2n\n 


j— r 

3.1416^1.6 


= 8.6440 10 13 : 


Le nombre d'onde est : 

» - i - 86440 10 ‘L 288335 m- 
C 2.9979 10 8 

La longueur d'onde (1Â = 1(T 10 m , lnm = lCT 9 m) est : 


2.9979 10 8 
8.6440 10 13 


= 0.3468 10" 3 m -1 


v = 2883 cm" 1 


X = 34680 A° 
X = 3468 nm 


Le nombre d'onde (v = 2883cm 1 ) et la longueur d'onde (X = 34680 A°) appartient à la 
région infrarouge du spectre électromagnétique. Rappelons que le spectre 
électromagnétique se subdivise en six importantes régions comme suit : 


0.1 

1 10 

10 2 10 3 

10 4 10 5 

10 6 10 7 

10 8 

! i ' 

V (Ctrl) 

^ i i 




X(A°) 

10 9 

10 8 10 7 

10 6 10 5 

10 4 10 3 

10 2 10 

1 

RY 

RX 

UV + v 

IR 

MW 

FR 








R y : rayons gama RX ■' rayons X UV + V : ultraviolet + visible 

IR : infrarouge MW ■ microonde FR • fréquences radio 

Le domaine infrarouge se situe entre 2x1 0 4 À et 5 x 10 5 A en longueur d'onde ou entre 200 
cm" 1 et 5000 cm" 1 en nombre d'onde. 


(b) Calcul des constantes de force des deux isotopes HCl et DCL : 


v = 



k = 4 ji 2 c 2 (iv 2 
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Pour cela on peut utiliser le système CGS : 

|v H ci = 2883 cm -1 ^ [v DC i = 2068 cm -1 

tact =1-6272 10 -24 g 6 I I^DCl = 3.1626 10“ 24 g 


D'où : 

k HCi = 479871 dyn/cm = 479.9N/m 
k DCl = 479887 dyn/cm = 479.9N/m 

On voit que k HC1 = k DC1/ c'est un résultat général. Une constante de force qui décrit la force 
da la liaison ne dépend que de la nature de la liaison chimique, elle est indépendante des 
masses. Une constante de force décrit la force da la liaison, plus k est grand plus la liaison 
est forte. Ainsi, deux isotopes possèdent des fréquences différentes à causes des masses 
différentes. Le rapport isotopique donne : 

h I k 

v hci =r — J 

2nc \ |a HC1 

h I k 

V DC1 =“ — J 

2 tic \ |LX Da 



H-hci _ m H m p + m ci _ l + ( m d /ai D ) 
M'DCi m D m H +m ci l + (m cl /m H ) 


Comme g DC1 > |i HC1 , la substitution isotopique provoque un déplacement du nombre 
d'onde vers une valeur plus petite v DC1 < v HC1 . 

(c) La molécule CN et son ion CN+ possèdent la même masse (toute la masse d'une espèce 
moléculaire est concentrée dans son noyau) : 


D'où : 



12.011x14.007 
12.011 + 14.007 


= 6.4662 uma = 1.07357 10“ 23 g 


k = 47r 2 c 2 |av 2 = 0.3808 v 2 


J k CN =16301 14 dyn/cm = 1630 N/m 
jk cN+ =950629 dyn/cm = 950 N/m 


Comme k CN > k CN+ , cela veut dire que la liaison chimique de CN est plus forte que celle 
de CN + car la structure électronique des deux systèmes est différente. Le calcul de l'ordre 
de liaison confirme ce résultat : 

OL = — (n-n*) 

2 


En effet, les configurations électroniques des deux systèmes donnent : 

CN : (la) 2 (la*) 2 (2a) 2 (2a*) 2 (ht) 4 (5a) 1 OL = |(9-4) = 2.5 

CN + : (la) 2 (la* ) 2 (2a) 2 (2a*) 2 (Iti) 4 OL = ^(8-4) = 2 

Ainsi OL(CN) > OL(CN + ) est en accord avec k CN >k cN+ : plus la liaison est forte plus sa 
constante de force est grande. 
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Problème III - Représentation matricielle de l'opérateur position : 

Soit T a (x) , y b (x) et ^ c (x) trois états quantiques de l'oscillateur harmonique : 



■l)e-P x2/2 


1 °- Une fonction d'onde de l'oscillateur harmonique est généralement de la forme : 


'Pu (x) = N d H u (Vpx)e _|3x2 ' 2 

où Hu(Vpx) est le polynôme d'Hermite de degré u en x. Ainsi, pour ^(x) ce polynôme 
est une constante donc u = 0, pour ^(x) il est de degré 1 en x donc o = 1 et pour x F c (x) il 
est de degré 3 en x donc u = 2 . Les énergies propres Eu = (u + ^jhv associées sont : 

'Fa(x) o = 0 -> E 0 = |hv 

'PbOO o = l -> E l = |hv 

'Pc(x) o = 2 E 2 = |hv 

2 °- (a) L'opérateur x = x dans la base ( v P a ,'P b , v P c ) est représenté par la matrice suivante : 

^ x aa x ab x ac" 


(x) = 


(b) La matrice est symétrique Xÿ =Xjj, car l'opérateur x = x est hermitique et les fonctions 
d'onde ¥3,^ et 'P c sont réels. En effet : 


x ÿ = ('P i |x|'P j ) = J'P*x'Pjdx = J'Lx'Pjdx = J'PjX'Pidx = (^ j |x|'ï r i ) = Xji 


Il suffit donc de ne calculer que six éléments au lieu de neuf : 

Xa x ab x ac" 


(x) = 


x bb x bc 

x cc , 


On utilise la propriété de parité pour prévoir les termes nuis et non nuis. En effet, on a : 


D'où : 


x = x : impair 


X : 

pair 

X: 

impair 

X : 

pair 


x aa 

= X X 'Pa) 

— > 

<P a x Ya 

impair 

— > 

x aa — ® 

x ab 

= Xl x X> 

-» 

*P a x Ÿb 

pair 

-> 

x ab * 0 

x ac 

= (TaXTc) 

-> 

'P a x 'P c 

impair 


x ac = 0 

x bb 

= Xl x X| 

-> 

'î'bX'î'b 

impair 

-> 

x bb = 0 

x bc 

= Oi'bM'i'.) 

— » 

'LbX 'î'c 

pair 

-> 

x bc * 0 

x cc 

-i'î'cbi'Pc) 


V c x 

impair 

-> 

x cc = 0 


Ces considérations préalables de symétrie permettent donc de simplifier grandement le 
calcul matriciel : 
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(x) = 


U x ab U 
s ab 0 x bc 
0 x bc 0 


Au lieu de calculer au départ 9 éléments, on ne calculera que 2 éléments x ab et x bc : 

^b = (^aM'l'b) = C x2e - Px2dx = 2|3 (f) J 0 + °° x2e “ Px2dx = ^ 

x b c = (^blx^c) = -^[2(3j_ + “x 4 e-P x2 dx-j_ + ”x 2 e-P x2 dx] 




-C 


4P 2 3 Jr 5 


# 


D'où la matrice représentative de l'opérateur position dans la base des trois premiers états 
propres de l'oscillateur harmonique : 


00 = -f= 

# 


0 

1/V2 

0 


1/V2 o N 
0 1 
1 0 


(c) Pour déterminer les valeurs propres de (x) , on calcule le déterminant caractéristique 
associé : 


lx-a.il = -!■ 

-x 1/V2 0 
1/V2 -X 1 

= 0 

1 1 # 

0 1 -X 


Les fonctions propres s'obtiennent : 


' -x 1/V2 0 N 



1/V2 -X 1 

c b 


0 1 -X 

V c cy 


= 0 


fx-o 

=> 

\x = -j3Ï2 

U J 


[X = +yl3/2 


-À,c a +-Uc b =0 

V 2 

-^=c a -À,c b +c c = 0 

V 2 

c b -À,c c = 0 


D'où les solutions normées : 

A, = 0 

X = ±yfv2 



¥ = 0 


On voit que pour les deux valeurs propres X = +V3/2, la fonction d'onde est nulle, ce qui 
veut dire que la probabilité de leur mesure est nulle. Ainsi, la seule valeur de la position 
moyenne réellement observable correspond à X = 0 c'est-à-dire (x) = 0 . L'état quantique 
correspondant est : 

[2 1 

Y = Ca'L.+Cc'L. = J-'F a — 

a a c c h a V3 


qui est en fait un mélange des deux états pairs : ^(x) fondamental ( u = 0 ) et T c (x) excité 
( u = 2 ). Le poids de participation de chaque composante est : 

|c a | 2 = ^xlOO = 67% : état fondamental v P a (x) participe avec un poids de 67% 

|c b | 2 |= 100- |c b | 2 = 33% : état excité v F c (x) participe avec un poids de 33% 
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3°- On suppose qu'une particule se trouve strictement dans l'état 4 y b (x) . 

(a) La position moyenne (x) de la particule dans cet état n'est autre que l'élément 
diagonale x bb de la matrice (x) : 

( x ) = x bb = 0 

(b) Pour déterminer les points nodaux (densités milles) et les positions les plus probables 
(densités maximales) associés à l'état y b (x), on doit calculer la dérivée première de la 
densité de probabilité : 

D(x) = | v F b (x)| 2 = x2 e^ x 


ÔD 

Sx 



[2x-2px 3 ]e~ px2 


= 0 -> 


x = 0 

x = -l/VP 
x = +i/Vp 


On obtient alors trois extrema (minima, maxima) aux quels correspondent les densités de 
probabilité suivantes : 

[D(0) = 0 

D(-l/ ^/p) = (p/7i)*2e -1 =0.2lVp 

[ü(+ 1/VP) = (P/7 i)’2e _1 = 0.21^/p 


Ainsi, on trouve un seul point nodal à la position x = 0 et deux points de maximum de 
densité aux positions x = -1/^P et x = +1/ ^/P dites postions les plus probables. 


Problème IV - Théorème du viriel : 

Soit H l'hamiltonien d'un système à une dimension H = ^ 2L + V(x) avec V(x) = A,x n . 

2m 

1°- (a) Calcul du commutateur [h, xp x ] : 

[H, xp x ] = ^-+V(x), xp x ] ~ x Px] + [V(x), xp x ] 

» x*-[]Px. Px] x + Hl x ] + [V(x), p x ]x + p x [V(x), x] 

2m L J 2m L J 

= I^IpL x ] + [V(x), P x ]x 
2m L J 

- id 2 *g , | _ 


2m L 


(b) Valeur moyenne de [h, xp x ] : 

( [H, xp x ] ) = ('P| [H, xp x ] |T) = |('P|2T-x^|'P) = ^2(^1 T |T)-('P|x^|'P)J 


([H,xp x] )^|2(T)- x ^ 
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Par ailleurs, la relation de Heisenberg nous enseigne que ( voir TD) : 




D'où : 


d ( X Px’ 

dt 


= 2 W- 


Pour un état stationnaire, ou encore pour n'importe quel état propre de H, la valeur 
moyenne (xp x ) est indépendante du temps. D'où l'expression générale du théorème du 
viriel : 



O y n-1 

(c) Dans le cas d'un champ de potentiel V(x) = A,x n , on a — = nÂx et la relation du viriel 
devient : 


2<T) = n(V) 


Une telle relation décrit l'équilibre entre l'énergie cinétique et l'énergie potentielle d'un 
système dans un état stationnaire T (fonction propre de l'hamiltonien H ). Autrement dit, 
la condition 2(T) = n(v) n'est vérifiée que si *P est solution exacte du problème. 

(d) Pour l'oscillateur harmonique, l'énergie potentielle est une fonction homogène de 
degré n = 2 en x ( x étant le déplacement de la masse m ) : 

V(x) = Ikx 2 => (T) = (V) 

Ainsi, le théorème du viriel pour un champ de potentiel harmonique indique que les 
énergies cinétique et potentielle moyennes son égales. Dans ce cas, l'énergie totale 
(E) = (T) + (v) n'est autre que le double de l'énergie cinétique ou de l'énergie potentielle : 

<E) = 2(T> ou <E}=2(V) 

Pour l'atome d'hydrogène, l'énergie potentielle est de degré n = -1 en r (r la distance 
séparant l'électron du noyau) : 

V(x) = -îî => 2{T) = -(V) 

r 

Dans ce cas l'énergie totale moyenne est : 

<E) = -2(T) ou (E) = i<V) 

4°- (a) La forme de la fonction d'onde exacte de l'état fondamental ls de l'atome 
d'hydrogène (exprimée en u.a) est : 

'¥ Is =(l/nf' 2 e- 

L'application du théorème du viriel à cette fonction d'onde stationnaire est rigoureuse, 
elle donne 2 (t) = -(v) . 

(b) Dans le cas de l'utilisation d'une fonction approchée, telle que la fonction gaussienne : 

=(2a/7i) 3/4 e“ cu ' 2 


on doit rechercher séparément les valeurs moyennes cinétique (T) et potentielle (V) : 
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(T) = (T ls |T|V ls ) = J ^T^dx, (V) = |V(r)|H> ls ) = J 'F 1 * s V(r)'F ls dx 

1 d 2 d _ d 2 2 d 


V(r) = , A = 


■ 2 dr dr dr 2 r dr 


dx = r sinOdrdOdcp 


En utilisant le système des unités atomiques (h = 1, m e = 1, e = l), on aura : 


Calcul de l'énergie cinétique moyenne : 
2 21 dr 2 r dr 


d^ls 
I dr 
1 d 2v E] 
[“d? 


= ~2ar v F ls 

- = -2a(l - 2otr 2 ) v F ls 


T'Pjs = (3a-2aV)'ï r ls 


il vient : 

{T) = (Y ls |3a-2a 2 r 2 |Ÿ ls ) = 3a(Y ls |Y ls )-2a 2 (Y ls !r 2 |Y ls ) = 3a-2a 2 (¥ ls |r 2 |Ÿ ls ) 

avec : 

(^ls |r 2 l^ls) = J ^Xdx = sin 0drd0dcp = r dr sin 0d0j ( ^drp 

* '**%'**& - Kt ) 1 Jo Ve " 2at ’ * = Kf] ’ l(^ôy) = ia 


D'où : 


< T > = 3a-2a 2 (A] 


<T) = fa 


Calcul de l'énergie potentielle moyenne : 

1 vt/ 2 \3/2 -2ar 2 

(v) = ('Fis |V(r)|'F ls ) = -(W ls |i|^) = -J-^dx = J JJJ — ^ — r 2 sin0drd0dcp 


*-(v] H“ 


re _zar dr = - — — 4rc- 


2(2a) 




Le théorème du viriel nous enseigne, dans le cas d'un champ coulombien, que 2(T) = -(v). 
Le paramètre a qui doit satisfaire cette égalité est : 

-a = f— 1 => a = — = 0.2829 

2 [n J 9 n 

L'énergie de l'atome d'hydrogène E = (T) + (v) pour cette valeur de a est : 


. n J 
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On obtient une énergie (-0.4244 hartree) plus haute que la valeur exacte (-0.5 hartree) de 
0.0756 u.a, c'est-à-dire de 2.06 eV (1 hartree = 27.211 eV). Ce résultat est attendu puisque la 
fonction d'onde Y ls utilisée n'est pas une fonction d'onde exacte pour l'état fondamental 
de l'atome d'hydrogène : 

'F ls (exact) = 0.564 e -r , Y ls (appro.) = 0.276 e -a283r 


(d) La méthode des variations consiste à minimiser l'énergie par rapport au paramètre a 
inconnu : 


il vient : 


avec 



SE _ 3 
da 2 



a 


_ 8 _ 

9tt 


On obtient une valeur quasiment identique à celle prévue par le théorème du viriel. Les 
deux démarchent sont valables pour rechercher les solutions approchées. 
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